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GETfERAL  E/TRODUCTIOI? 

IMa  report  contains  a  series  or  technical  discussions  and  papers  con¬ 
cerned  vith  the  theory  of  vave  propagation  in  solids  vith  special  applica¬ 
tion  to  ground  shock  phenomena,  It  is  a  summary  vhich  provides  the  theo¬ 
retical  background  to  scene  of  our  investigations  on  these  topics,  initiated 
in  the  early  part  of  1957,  and  vhich  are  still  actively  pursued  at  the 
present. 

The  report  deals  vith  the  tvo  major  subjects  of  interest  namely,  free 
field  effects  and  diffraction  phenomena  both,  in  regions  vhere  neither 
themodynanic  effects  nor  very  large  displacements  have  significant  influ¬ 
ence  on  the  overall  physical  picture.  This  restriction  eliminates  the 
crater  and  its  iiraediate  vicinity  from  the  province  of  our  investigations. 
This  means  that  effects  outside  of  the  crater  region,  but  caused  by  phen¬ 
omena  vithin  this  region  are  only  crudely  represented  by  our  analytical 
vork.  Numerical  evaluations  of  those  effects  should  preferably  be  started 
vith  iLput  data  obtained  from  these  close-in  calculations,  The30  resear¬ 
ches  (conducted  by  3 rode  at  The  RAND  Corporation)  vere  not  conclude!  dur¬ 
ing  the  present  project,  and  for  this  reason  the  validity  of  some  of  our 
conclusions  must  be  bracketed  vithin  suitable  space  and  tine  intervals. 

Attention  mot  be  paid  to  the  relation  of  the  mathematical  model  to 
the  expected  physical  behavior  of  the  material.  Our  researches  show  that 
many  of  the  significant  quantities  due  to  fundamental  consideration  are 
not  very  sensitive  to  variations  of  physical  models.  For  example,  the  P 
vave,  generated  by  a  surface  burst  in  an  elastic  solid,  appears  in  an 
only  slightly  changed  fom  and  intensity  in  an  elastic  fluid,  and  it  is 
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found  again  vlth  Gmail  notifications  in  a  viococlootic  medium.  Similarly 
in  granular  solids,  tho  essential  attenuation  of  porticlo  velocities  and 
otrcoocs  is  not  too  ocnoitivc  to  quite  large  variations  cf  physical  para¬ 
meters  . 

In  the  case  of  Rayleigh  vaveo  vc  vcrc  able  to  ectablish  the  region 
beyond  vhich  they  do  not  contribute  cignificontly  to  tho  freo  field  ctrecoes, 
even  though  th$  onplitude  of  these  vaveo  obtained  by  clastic  theory  vithin 
this  region  i'-ccif  ic  larger  than  vhat  the  physical  properties  of  the 
notorial  con  Justify.  While  these  cccputations  furnish  only  upper  bounds 
of  the  quantities  required,  their  implications  are  physically  meaningful 
and  of  practical  significance. 

The  diffraction  problem  has  so  far  only  been  considered  in  tho  elastic 
medium,  and  it  has  led  to  important  conclusions  regarding  the  anplificr- 
tion  of  free  field  stresses  around  underground  openings.  Ihcoo  go me  effects 
brought  out  ooae  nev  results  regarding  the  applicability  of  frequency- 
amplitude  (or  shock)  spectra:  it  is  found  that  at  high  frequencies  the 
diffracted  spectrum  is  significantly  different  from  that  of  the  free 
field  spectrum. 

Ihe  duplications  and  the  applicability  of  theoretical  findings  are 
discussed  in  thi**  report  and  more  detail  is  contained  in  some  earlier 
studies  ( [  17 1 ,  [18],  and  [19])*  We  ccroe  to  the  conclusion  that  vithin 
the  limitations  outlined,  the  results  of  our  theoretical  work  can  be 
applied  vlth  confidence  to  practical  problems  dealing  vith  ground  shock 
effects. 

lb©  above  work  vaa  undertaken  to  provide  urgently  nooded  answers  on 
problems  related  to  the  vulnerability  of  underground  structures:  florae  of 
the  earlier  classified  feasibility  studies  ([l8]  and  [19] ), prepared  under 
pressure  of  deadlines,  were  issued  containing  only  numerical  results,  or 
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educated  guesses,  based,  however,  on  rather  cxtenoive  theoretical  researches * 
The  purpose  of  the  present  report  therefore  is  to  supply  in  sane 
detail  the  background  and  data  on  which  our  previous  recccsncndaticnc  were 
based.  It  attempts  to  gather  oexae  of  the  loose  ends  and  to  ourzoarize 
our  current  knowledge  gleaned  from  our  own  investigations  and  from  the  work 
or  others.  The  theoretical  results  are  discussed  in  the  report  itself, 
while  some  of  the  detailed  derivations  arc  in  the  Appendices. 

A  considerable  port  of  the  information  which  is  summarized  here  is 
baaed  on  our  earlier  researches  which  li&ve  been  published  or  arc  tovbe 
published  in  the  near  future  in  technical  Journals'  ,  but  recent  theore-  , 
tical  results  which,  at  the  present  are  not  available  in  other  fora  ore 
contained  in  Appendices  C,  D  and  E. 

8otae  of  our  current  investigations  which  could  not  be  included  at  the 
present  time,  but  which  will  be  reported  in  the  near  future  are  a3  follows: 
Static  and  quasi- dynamic  high  pressure  exp^iontal  votk 
on  granular  media; 

Wave  propagation  in  vioco-eloctic  media. 

In  o~lte  of  the  rather  extensive  coverage  of  the  various  topics 
which  ha\e  been  considered  by  us  and  many  other  researcher'?,  there  arc  a 
great  number  of  gaps  in  our  knowledge,  which  have  important  implications 
cm  the  general  problem  of  vulnerability  of  underground  structures.  The 
most  significant  which  are  in  need  of  urgent  solution,  are  as  follows: 

Failure  mechanism  and  failure  criteria  of  cavities  in 
elastic  media. 

Diffraction  of  stress  waves  in  non-circular  openings. 


(#)  These  are  denoted  by  (#)  in  the  list  of  References. 
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Evaluation  of  the  effects  of  tunnel  linings. 

Three-dimensional  wave  propagation  in  granular  media. 

Diffraction  of  stress  waves  in  granular  media. 

Rayleigh  wave  attenuation  due  to  non-linear  and/or 
plastic  effects. 

Experimental  verifications  and  research. 

While  the  above  list  is  far  frcw  exhaustive,  it  is  uncomfortably 
long;  past  experience  with  similar  topics  leads  us  to  suspect  that  studies 
on  these  subjects  may  bring  about  surprising  and  possibly  aggravating 
discoveries.  While  ve  are  working  at  the  present  on  some  of  these  topics, 
we  are  not  as  yet  in  a  position  to  report  significant  progress  in  these 


areas. 
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FREE  FIELD  EFFECTS 
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FRKE-FIEL-D  EFFECTS 


1.  INTRODUCTION 

The  following  contains  a  summary  of  available  theoretical  results  on 
free -field  effects  due  to  (hypothetical)  Megaton  explosions  on  the  surface 
of  semi- infinite  media  of  various  idealized  properties. 

The  ultimate  purposes  of  the  presentation  arc  conclusions  for  actual 
media  like  rock  or  soil  with  complex  properties  difficult  to  analyze  or 
even  to  describe.  The  problem  of  rock- like  media  is  approached  by  consider¬ 
ing  a  succession  of  materials  having  gradually  more  complex  properties.  In 
this  manner  certain  conclusions  drawn  for  cases  of  simple  properties  can  be 
extrapolated  for  more  complex  ones  by  qualitative  reasoning.  The  writers 
believe  that  estimates  applicable  in  judging  realistic  situations  in  rock 
have  been  obtained  in  this  manner  by  considering  the  succession:  acoustic 
invlscid  fluid,  linear  elastic  solid,  non-lincar  elastic  colid* 

As  contribution  to  the  problem  of  soils  additional  studies  of  wave 
propagation  in  locking  media  have  been  made.  Most  of  the  work  concerns 
one- dimensional  wave  propagation,  but  the  case  of  spherical  waves  as  now 
also  been  treated. 

To  have  an  understanding  of  the  meunlng  of  the  analytical  results  it 
is  necessary  to  use  a  realistic  pressure  time  history  in  computing  examples. 
Such  a  history  is  indicated  in  Fig. ( 1-1),  which  represents  ''imputed  values 
of  surface  pressures  for  a  20MT  surface  burst  [l] -  It  is  pertinent  that 
the  pressure  at  any  time  acts  over  a  cir<nilar  area  of  radius  rQ  with  constant 
Intensity  except  at  the  periphery,  where  a  pressure  spike  appears*  It  will 
be  found  very  important  that  the  velocity  of  increase  of  the  radius  of  the 
circle,  i.e*  the  velocity  of  the  shock  in  air,  at  early  times  is  very  large, 
much  larger  than  the  seismic  velocities  in  the  ground.  The  shock  velocity 
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in  air  at  19  ms,  e.g.  is  still  V  »  34,000  feet  per  second,  and  therefore 
larger  than  the  seismic  velocities  of  any  likely  medium;  however,  at  290  ms 
the  velocity  has  decreased  to  V  «  7*200  feet  per  second  which  is  less  then 
the  seismic  velocity  of  some  rocks. 

For  qualitative  considerations  it  is  of  interest  to  have  information 
on  the  impulse  I  « J p  dt  of  the  surface  pressure.  Fig. (l-2) shows  X  at 
the  center  of  the  explosion  as  a  function  of  the  time  tQ. 

The  following  recurring  symbols  are  adhered  to  in  the  body  of  the 

* 

report 

r,  R  and  a  coordinates  as  in  Fig. (2-1). 

c  velocity  of  sound  in  fluid,  velocity  of  shear 

waves  in  linear  elastic  solid* 

c  ,  cR  velocity  of  dilatational  and  Rayleigh  waves,  respee* 

*  tively,  in  linear  elastic  solid. 

0,  \ ,  \i  modulus  of  rigidity,  Loa^s  constants,  respectively 


Different  symbols  may  be  used  in  the  Appendices,  as  they  are  taken 
from  various  sources. 
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HALF  SFAC&  Or  \00U3TIC  ILVICCID  FLUID 


*fhe  effect,  of  I  pressure  p(r,  t)  on  the  our'u^  oi  a  liaif  space  can 

(*) 

be  obtained  for  the  case  of  an  acoustic'  '  fluid  from  the  solution  for  a 
concentrated  point  load.  The  load  ia  a  unit  stop  load  in  :ime,  of  intensity 
H(t).  ISic  solutions  for  the  pressure  can  be  written  in  closed  form  [2] 


CQ3  0 
2  n  Be 


6  (t  -  |)  +  522®  H(t  -  |)  (2-1) 

c  2  K  R  c 


where  r,  R,  z  and  0  are  defined  in  Fig.(2-1). 
p  is  the  pressure  in  the  fluid 


H(x)  ^ 
fc  (r)  « 


0  if  x  <  0 


is  a  symbol  defining  a  unit  step  function. 


^(x) 

lx 


if  r  >  0 

i3  Disc’d  Delta  Function 


Using  the  applied  pressur.  p(r,  t)  defined  in  Fig.(l)and  c  ~  5000  feet 
per  second,  Figc.  (2-2)to(2-h)give  computed  pressures  for  points  r  *  0,  and 
r  «*  2000  feet  at  various  depths. 


a)  Effects  under  Ground  Zero,  r  0. 

Figs.  (2- 2i  (2-  3) consider  points  at  r  ~  0,  directly  under  the  '.enter  of 

the  explosion*  For  the  first  2.2  me  (Fig.  2)  the  given  pressure.,  at 

z  =  0  and  the  computed  ones  rt  5 00  feet  and  2000  feet  depth  are  undistinguish- 

able.  For  later  times,  Fig.<%  differences  occur,  the  deviations  occurring 

earlier  at  greater  depth.  It  2  s  significant  that  for  a  few  milliseconds 

(except  for  time  delay)  the  p rc  5ire  signal  Just  under  ground  zero  is 

independent  of  th;  depth.  Th<  lo  due  to  the  fact  that  the  pressure  p(r,  t) 

spreads  on  the  surface  with  a  large  velocity  that  one  could  aoeume 

simultaneous  application  of  a  uniform  pressure  over  seme  radius  (Fig.  1-1) 

*  Note  that’  in  acoustic  theory  only  first  order  terms  are  retained, 
such  that  the  equations  becceno  linear.  This  linearity  permits  the 
use  of  superposition. 
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ohovs  that  this  pressure  is  at  early  tinea  nearly  unifora,  the  spikes 
containing  only  a  email  portion  of  the  total  load). 

If  a  uniform  pressure  p(t)  is  applied  on  the  surface  of  a  semi- infinite 
fluid,  the  problem  becomes  one  of  one- dimensional  wave  propagation,  and  in 
the  acoustic  case  considered  here,  the  pressure  signal  propagates  then 
without  change.  If  a  uniform  pressure  is  applied  only  over  a  circular 
area,  a  point  directly  under  the  load  receives  a  signal  equal  to  the  applied 
pressure  for  a  limited  time,  later  the  pressure  drops  as  illustrated  in 
Fig.  (2-4)  for  r  «  0,  z  »  500.  The  point  where  the  pressure  drops  is  defined 
by  the  arrival  time  of  a  signal  from  the  periphery  of  the  loaded  circle. 

Hie  situation  in  Figs . (2) and  (3) is  similar,  but  duo  to  the  expansion  of 
the  loaded  area  one  cannot  find  an  exact  time  t  up  to  which  one- dimensional 
wave  propagation  applies;  instead  one  con  estimate  a  time  t  up  to  which 
one-dimensional  wave  propagation  should  be  a  good  approximation.  Hie 
time  in  which  the  approximation  is  valid  will  of  course  be  a  function  of 
the  depth,  z,  (and  of  the  seismic  velocity  and  weapon  size). 

Hiis  rather  elementary  explanation  was  presented  in  ouch  detail  because 
the  above  qualitative  reasoning  is  valid  not  only  for  the  case  of  a  linear 
elastic  solid,  but  also  for  complicated  non-linear  and/or  non- elastic  media. 
Even  in  such  cases  pressure  waves  due  to  the  blast  pressures  on  the  surface 
must  bt,  for  a  certain  time,  one- dimensional  in  nature,  permitting  manage¬ 
able  numerical  solutions.  These  solutions  should  be  good  approximations 
for  a  definable  range  of  time. 

b)  Effects  in  the  high-pressure  (supers^lsmlc)  range  when  z  «  r 

As  a  typical  example  Equations  (2-1)  were  used  to  obtain  the  pressure 
at  a  typical  point  r  ■  2000  feet,  z  *  500  feet  in  the  range  of  high  pres¬ 
sures.  The  result  is  shown  in  Fig.  (2- $) where  the  computed  points  are  compared 


vith  the  surface  pressure  at  the  some  radius,  r  »  2000  feat.  It  is  seen 
that  the  pressures  at  a  depth  of  500  foot  differ  only  little  fraa  those 
on  the  eurxaco. 

The  scull  differences  betveen  the  pressures  on  the  surface  and  at  a 

depth  z  «  x'  are  due  to  the  fact  that  the  velocity  V  of  the  pressure  wave 

in  air  passing  over  the  point  is  very  much  larger  than  the  seismic  velocity 

c  in  the  medium.  The  first  portion  of  this  pressure  signal,  oay  up  to  $0  ns 

after  arrival,  is  due  to  the  airprccourc  .passing  over  the  point,  the  very 

first  signal  originating  from  a  point  F,  see  Fig.(5&)»  A®  first  approx! - 

nation,  the  actual  pressure  on  the  ourfec*?  might  be  replaced  by  a  plane 

pressure  vave  progressing  vith  constant  velocity  V  >  c  and  without  chango 

in  shape  of  the  nrossui  ;  vave.  This  problem  has  an  elementary  solution; 

a  vave  on  the  surface  p  *■  f(t  -  x/v)  produces  an  inclined  plane  vave  in 

the  fluid,  p(t  -  j/c),  Fig.  (2-6).  The  angle  a  depends  on  the  ratio 

V/c,  and  the  pressure  signal  p  at  iuf>y  point  hao  the  same  shape  as  the 

z 

surface  pressure,  p^  ■■  f(t  -  t^);  where  t&  is  a  time  delay, 

In  the  high  pressure  range  where  the  velocity  of  the  shock  :'ront  in 
air  V»  c,  the  pressure  bolow  the  surface  in  tho  fluid  can  theroforo  bo 
approximated  by  the  surface  pressure.  By  purely  qualitative  reasoning  one 
can  conclude  that  this  must  be  00  for  a  very  short  time  after  arrival,  and 
if  the  depth  z  is  very  small.  It  is  therefore  important  to  note  that 
Fig, (2-5)  shove  that  a  reasonably  good  approximation  is  obtained  even  at 
a  depth  g  *  r/4  and  for  the  entire  time  range  where  the  pressures  arc 
substantial.  It  appears  that  z  «  r/4  is  sufficiently  small  to  satisfy  the 
original  assumption  z  «  r. 

T.t  is  interesting  to  discuss  the  difference  betveen  this  approximation 
and  the  computed  points.  At  first  arrival  the  actual  pressure  must  be  equal 


-lt>- 


lo  ih*  air  pr*»r«uro  at  th**  point  ?  (Fig.  >a)  cn  tae  surface,  which 

pr*«**ur«  lrt  a  XitO*  larger  than  the  peak  for  1  «  2v.aj  foot.  (One  could 
there  f one  improve  tho  approximation  for  tha  x>ca!t  proa o are  by  uaing  the 
pressure  hlfttory  at  F. )  Subsequently,  for  a  limited  *  irx- ,  the  computed 
oresauroo  du«  to  tho  loadn  according  to  Fig.(l-l)are  iewer  them  the  ourfaco 
P~ooi3vro»  similar  to  the  situation  for  r  **  0  m  Fig. (2-3)*  However,  at 
later  tinea,  t  >  1^0  js a,  the  computed  pressures  exceed  the  approximate 
ones  somewhat .  This  is  caused  by  the  hign  applied  pressures  vithin  th«< 
central  circle  r  ~  1000  feet,  shortly  after  the  explosion.  Waves  canerated 
from  this  high  pressure  area  arrive  at  this  time;  the  effect  at  the  point 
considered  is  however  only  small,  because  these  pressures  are  remote  and 
act  only  ovux  a  relatively  small  area.  A  computation  confirming  that 
the  effect  of  the  early  high  pressure  in  a  typical  case  is  not  significant 
was  ho rv*v#*r  required,  because  the  relative  magnitude  of  *he  effects  compared 
depends  on  the  rapidity  of  the  decay  of  the  applied  air  pressure.  For  a 
very  much  larger,  or  very  much  smaller  weapon,  the  situation  might  be  different. 

The  simplifying  approximation  of  assuming  that  tho  shape  of  the  sur¬ 
face  pressure  does  not  cnange,  and  that  this  pressure  progresses 
stont  velocity,  V,  has  been  tested  above  at  a  point  where  the  shock  velocity 

V  is  several  times  larger  than  the  seismic  velocity  c.  The  same  assumption, 
ve  will  call  it  the  "steady  state  solution",  can  be  made  with  good  results 
if  the  velocity  V  is  appreciably  smaller  than  c,  but  the  pressure  in  the 
fluid  will  not  have  the  character  of  Fig. (2-6).  A  theoretical  solution  for 
this  case  can  be  derived  from  £3] .  Without  pursuing  the  matter  further,  it 
should  be  stressed  that  the  approximation  must  necessarily  become  poor  if 

V  and  c  differ  little ,  regardless  whether  V  >  c  or  Y  c,  In  the  superceismic 
case  caution  ia  therefore  indicated  unlocc  V  £  1.5e*  Further  numerical  work 
to  find  kov  dose  to  V  »  c  the  approximation  may  be  used  is  planned. 
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Rcturning  to  the  ouperocismic  range,  it  io  important  that  the  analysis 
for  the  fluid,  baaed  on  the  approximation,  indicates  no  attenuation  of  the 
pressure  vith  depth,  (Hilo  conclusion  can  of  course  only  apply  for  the 
range  z«r  for  which  the  approximation  can  be  valid.)  If  the  approximation 
vc re  not  used  small  differences  vith  depth  vo\*J.d  bo  found,  but  there  vould 
be  no  appreciable  attenuation,  until  the  depth  becomes  comparable  to  the 
radial  distance  r.  (in  the  ouboeiomlc  range,  on  the  other  hand,  on  analysis 
using  the  some  approximation  indicates  substantial  attenuation. ) 

Thfe. situation  for  the  fluid,  vhich  in  itself  io  not  pertinent,  liac  been 
discussed  at  length  because  equivalent  approximation  a  can  he  made  for  the 
clastic  solid  or  other  media.  Math coat i colly  speaking  the  steady  state  sol¬ 
ution  in  the  auperseionic  range  is  on  asymptotic  approximation,  for  short 
tine  after  arrival;  because  the  air  pressure  signal  used  is  the  one  for  a 
point  on  the  surface  at  the  distance  r  fron  ground  zero  (not  for  point  F, 
Fig.  Jo,  where  the  first  signal  originates)  the  approximation  as  proposed 
here  is  also  restricted  to  sufficiently  small  depth  z ,  such  that  the  pres- 
auras  at  xh «  two  points  do  not  differ  noticeably. '  f 


It  vould  be  an  unwarranted  refinement  to  use  the  pressure  at  point  F 
as  basis,  because  the  location  of  this  point  vould  require  a  cumbersome 
computation  for  each  depth  z.  Yet,  if  the  actual  air  pressure  varies  so 
rapidly  that  the  pressure  at  F  differs  appreciably  from  the  pressure 
exactly  above  the  point  considered,  the  signal  vould  be  affected  by  these 
variations  within  a  short  time  after  first  arrival.  In  such  a  case  the 
steady  state  approximation  vould  be  poor  and  useless  anyvay. 
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3.  HALF  SPACE  OF  LIME AR-EIASTIC  feOLID. 

The  effect  of  on  applied  surface  pressure  according  to  Fig.(l-l)can 
ho  obtained  again  by  superposition  from  the  effect  of  a  concentrated 
step  load*  However,  the  effect  of  the  latter  is  much  more  complicated 
than  in  the  case  of  an  acoustic  fluid* 

a)  Effect  of  a  concentrated  step  load* 

The  literature  {4],  C5 ]  contains  formal  solutions  for  the  displacements 
of  the  surface  due  to  buried  horizontal  and  vertical  step  loads.  By  virtue 
of  a  dynamic  reciprocal  theorem  (6),  (7 1  it  is  concluded  that  the  vertical 
and  horizontal  components  of  the  displacements,  due  to  a  step  load  on  the 
surface,  are,  respectively,  equal  to  the  vertical  displacements  of  the  sur¬ 
face  due  to  vertical  and  horizontal  step  loads  acting ‘below  the  surface, 
see  Fig.  (3-l)»  The  above-named  papers  furnish  formal  expressions  for  the 
displacements  for  the  present  problem;  expressions  for  the  strains  and 
stresses  can  be  obtained  by  routine  manipulations.  However,  the  expressions 
ore  not  in  closed  form,  they  require  the  very  cumbersome  numerical  evalua¬ 
tion  of  definite  integrals.  A  series  of  displacement  histories  w  >  computed 
and  graphically  presented  in  14].  While  extremely  helpful  for  purposes  of 
the  general  discussion  which  follows,  these  graphs  could  not  be  utilized 
for  the  numerical  determination  of  stresses  or  displacements  due  to  the 
Surface  pressure,  Flg.(l-l).  Numerical  results  reported  below  had  to  be 
developed  from  the  available  theoretical  expressions. 

General  information  on  vrave  propagation  phenomena  iu  elastic  solids 
can  be  found  in  standard  texts,  defining  P- waves  of  dilatation  (velocity  Cp) 
and  S-vaves  of  distortion  (velocity  c<ep).  In  addition  the  existence  of 
surface  waves,  called  Rayleigh  waves  (velocity  Cp<c),  is  established. 

Knowledge  of  theso  vavo  types  permits  a  purely  qualitative  description  of 
the  rftate  of  stress  in  a  iialf  space  due  to  a  concentrated  step  load,  Fig.  (3-2)  • 


At  a  given  tine  t  after  application  of  the  load,  the  discontinuities 
associated  vith  the  fronts  of  the  P  and  S  waves  will  have  reached  the 
surface  of  spheres  of  diameters  tcp  and  tc,  respectively.  Further,  the 
region  between  the  S  and  P  fronts  is  divided  by  a  conical  surface  SP, 
locating  the  Van  Schmidt  front.  The  Rayleigh  waves  will  be  of  importance 
on  the  surface  at  or  near  R,  somewhat  behind  the  S- front. 

Considering  a  point  r,  z,  Fig. (3-^>  the  wave  fronts  will  progress 
vith  time;  the  P  and  S  fronts  will  pass  through  any  point,  while  the  SP 
front  will  only  pass  through  points  for  which  the  ratio  z/r  is  less  than  a 
critical  value  which  is  a  function  of  the  ratio  cs/cp.  Rayleigh  waves 
will  be  only  noticeable  near  the  surface,  i.e.  if  z/r  is  very  small.  This 
qualitative  description  is  borne  out  by  Fig.  (3-^  giving  the  vertical 
deflection  v  «  vp  according  to  Reference  (h)  for  three  ratios  z/r. 

For  r/z  «  0.5  only  the  P  and  S  waves  pass  the  point;  for  r/z  =  5  the  SP 
front  appears,  and  there  is  a  faint  trace  of  the  Rayleigh  wave;  for  r/z  «=  40 
the  Rayleigh  wave  is  fully  developed  and  is  a  major  effect. 

It  will  be  noticed  that  the  response  curves  in  Fig.(3-3)&t  the  P  and  S 
fronts  show  peculiarities,  i.e.  discontinuities  or  infinite  values.  These 
are  due  to  the  singularity  in  the  mathematical  fiction  of  the  assumed 
"concentrated  point  load. 11  If  the  effect  of  a  physically  possible  dis¬ 
tributed  load  is  computed  by  integration  over  the  response  from  the  con¬ 
centrated  point  load,  it  is  found  that  undue  peculiarities  disappear. 

Response  curves  for  other  displacements  or  stresses  will  obviously 
bo  of  a  similar  complicated  nature,  having  various  ranges  as  functions  of 
r/z  and  of  the  non  dime  no  ional  time  t  =  ct/R.  To  find  the  response  due 
to  the  pressure  defined  in  Fig.(l-l)by  a  purely  numerical  integration  process, 
covering  the  entire  response  range  in  r/z  and  t  ,  would  become  extremely 
involved  and  time  consuming  even  using  modern  computer  equipment.  For  this 
reason  an  attempt  will  be  made  to  recognize  the  important  effects  in  each 
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location  of  interest  and  to  find  suitable  approximations,  similar  to,  or 
more  general  than  those  developed  in  Section (2)  for  the  liquid  half  space. 

b)  Effecto  under  ground  zero,  r  ~  0 

Consider  a  point  A  at  r  ^  0  at  some  depth  z  »  r,  Fig. (3-4),  at  a 
time  t  quite  soon  after  arrival  of  the  first  signal  from  a  surface  pres¬ 
sure  p(r,  t)  of  the  type  shown  in  Fig. (l-l).  At  this  time  only  the  surface 
pressure  acting  within  a  circle  of  radius  p  ,  (sec  Fig.  3-4)  con  con¬ 
tribute  to  the  response  of  point  A.  The  integration  over  the  effects  of 
a  concentrated  load  on  the  surface  vill  then  involve  only  curves  which 
have  the  nature  of  the  curve  for  r/z  =  0.5  in  Fig.  (3-3).  To  define  the 
response  in  this  range  (r/z  <  l/v^S  )  the  Jumps  at  the  P  and  S  fronts 
were  obtained  (in  closed  form)  from  the  theory.  The  result  for  the  ver¬ 
tical  displacement  v  io  shown  in  Fig.  (3-5)*  The  portion  of  the  curves 
between  the  P  and  8  fronts,  and  the  portion  following  the  8  front  are 
graphically  given  in  [4].  They  are  30  close  to  straight  lines  that  they 
were  approximated  by  straight  lines  to  obtain  a  single  Idealized  response. 

(The  level  of  the  horizontal  lino  for  — ■  l  is  defined  by  the  s?  stic  deflec¬ 

tion,  known  from  Doussineoq's  solution.) 

Using  thi6  approach  the  vertical  velocity  v  of  a  point  at  z  =  1000  foot 
was  computed  for  several  values  of  t.  The  results  are  shown  in  Fig.  (3-6) 
and  compared  with  the  vortical  velocity  vq  on  the  surfs.ee.  It  is  seen  that 
in  the  time  range  plotteu  the  elocity  histories  on  the  surface  and  at  the 
depth  of  IOOO  feet  eannot  be  distinguished.  The  situation  corresponds  to 
the  one  shown  in  Fig. (2-2) for  the  fluid:  and  we  draw  again  the  conclusion 
that  in  this  location,  r  *  0,  the  assumption  of  one- dimono lonal  iravc  propa¬ 
gation  will  give  a  good  approximation  for  the  important,  early,  high-pressure 
part  of  the  signal. 

For  a  linear-elastic  material  one- dimensional  wave  prop^ig&tion  implies 
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no  attenuation  of  pressure  with  depth;  this  conclusion;  however,  must  not 
be  applied  to  actual  materials  which  will  behave  non-lincnr  at  high  pres¬ 
sure  levels  near  ground  zero.  The  non-linear  case  is  discussed  later > 
Section (5) * 

c)  Effects  in  the  high-pressure  ( supersoismic )  range  at  moderate 
depth,  z  «  r 

In  the  corresponding  case  of  the  fluid,  treated  in  2  (b),  it  was 
found  that  the  response  for  a  considerable  length  of  time  (several  hundred 
ms)  could  be  approximately  determined  from  a  steady  state  solution*  This 
solution  was  the  response  due  to  an  applied  plane  pressure  p  =  f(t  -  x/V), 
wherfe  f(t)  is  the  pressure  history  on  the  surface  above  the  point  consid¬ 
ered,  and  V  is  the  velocity  of  the  3hock  in  air  at  this  point. 

Using  the  same  general  reasoning  as  in  the  case  of  the  fluid,  one 
comes  again  to  the  conclusion  that  at  moderate  depth,  z  «  r,  the  earty  por¬ 
tion  of  the  response  in  the  olaotic  solid  can  bo  obtained  with  good  approxi¬ 
mation  by  computing  the  response  to  the  simplified  steady  pressure  pulse 

p  ®  f(t  -  x/V).  The  solution  for  this  problem  is  available  [8]  and  enclosed 
(*) 

as  Appendix  Ax  .  However,  it  must  not  bo  concluded  without  further  inves¬ 
tigation  that  this  approximation  will  be  sufficient  for  a  considerable  time, 
as  in  the  previous  example  for  the  fluid.  When  discussing  the  latter  case, 
Fig* (2-5^  it  was  pointed  out  that  the  computed  pressures  for  t  >  200  ms 
were  larger  than  the  approximate  values,  which  vac  ascribed  to  effects 
originating  in  the  region  of  very  high  pressures  near  r  «  0  shortly  after 
the  explosion. 

(#)  Thero  is  a  printing  error  in  [8]  for  the  cuporseicmic  case,  which 
is  corrected  in  the  Appendix. 
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To  have  a  basis  for  Judging  the  accuracy  of  the  response  of  the 
elastic  solid  at  later  times,  it  is  necessary  to  consider  effects  originat¬ 
ing  in  the  region  of  high  pressure.  This  could  be  done  by  a  complete 
numerical  determination  of  the  response  to  a  step  load  from  {4]  and  t5] 
and  subsequent  integration  for  the  prescribed  load  p(r,  t).  The  very  ex¬ 
tensive  numerical  work  required  is  however  not  required.  By  qualitative 
consideration  of  the  energy  available,  one  can  reason  that  the  effects  of 
the  P  and  S  waves  from  the  region  of  high  pressures  will  be  of  the  come 
order  of  magnitude  as  the  pressure  wave  effects  in  case  of  the  fluid, 
and  as  far  as  the  effects  of  P  and  S  waves  go  the  steady  state  approxima¬ 
tion  in  the  solid  should  be  as  good,  or  as  poor,  as  in  the  fluid.  This 
reasoning  does  not  include,  however,  the  energy  transmitted  by  surface 
waves.  It  is  obvious  that  the  energy  of  surface  waves  x^osing  through  a 
point  at  some  distance  r  will  be  a  function  of  l/r,  while  the  energy  in  P 
and  S  waves  decreases  as  l/r2.  Rayleigh  wave  effects  decay  therefore 
slower  than  the  other  effects  and  cannot  be  Judged  on  the  basis  of  an 
example  for  the  fluid.  Stresses  or  displacements  caused  by  Fayleigh  waves 
originating  in  the  high  pressure  region  must  therefore  be  computed,  and  if 
substantial  they  nruot  be  added  to  the  "steady  state  solutions *"  As  Rayleigh 
wave  effects  decrease  rapidly  with  depth,  their  relative  importance  compared 
with  the  effects  described  by  the  steady  state  solution  will  depend  on  the  • 
depths  of  the  point  considered,  For  sufficiently  deep  points  the  Rayleigh 
effects  will  become  negligible,  but  the  depths  where  this  occurs  are  not 
known  beforehand.  It  was  therefore  necessary  to  develop  an  analytical  pro¬ 
cedure  for  the  determination  of  Rayleich  (surface  wave)  effects. 

This  suggests  that  the  response  in  this  range  may  be  obtained  approxi¬ 
mately  by  using  the  "steady  ocatc  solution"  and  adding,  if  necessary,  the 
Rayleigh  wave  effects.  It  is  now  appropriate  to  consider  the  steady  Btate 
solution  and  its  limitations. 
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Discussion  of  Steady-State  Solution 

If  a  proscuro  vavG  p(t  -  x/v)  travels  with  constant  supcrceismic 
velocity  cp  >  c  on  the  surface  of  a  (linear)  elaatlc  solid,  it  pro¬ 
duces  on  inclined  plane  P-vave  and  an  inclined  piano  S-vavc  in  the  solid, 
as  shown  in  Fig.  (3-7).  The  intensity  of  these  pressure  ana  shear  waves  is 
proportional  to  p(t),  the  factors  and  angles  of  inclination  depend  on  the 
ratio  V/cp  and  on  Poisson's  ratio,  sec  Appendix  A.  Fig.  (3-7) chows  two 
cases  Y/cp  3  1.5  and  3i  both  indicating  that  the  P  wave  has  nearly  the 
same  intensity  as  the  applied  surface  pressure  while  the  S  wave  is  less 
important.  The  theoretical  solution  boo  the  character  of  Fig.  (3-7) os  long 
as  V  >  Cp*  However,  aa  the  velocity  V  approaches  Cp  the  strengths  of  the 
waves  vary  rapidly  aa  function  of  V/cp,  as  seen  from  the  following  table  v 
(for  Poisson's  ratio  l/4). 


V 

P-Wave 

S-Wave 

cp 

o/p 

f/p 

3.0 

0.99 

0.22 

1*5 

0.96 

0.45 

1.4 

0.96 

0.49 

1*3 

0.97 

0.52 

1.2 

0.99 

O.56 

1.1 

1.05 

0.59 

1.05 

1.19 

O.58 

1.01 

1.71 

0.46 

1.001 

2.4C 

0.21 

1 

3 

0 

{#)  Bie  situation  is  also  presented  graphically  in  the  Appendix,  Fig.  A-2. 
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The  sensitivity  of  the  solution  for  ratios  V/cp  <  1.1  indicates 
clearly  that  tho  steady- state  solution  for  nuch  ratios  cannot  be  used  as 
approximation.  One  can  easily  give  an  alternative  reason  t/hy  the  steady- 
state  solution  be  cane  a  poor  and  invalid  as  V/c^  -♦  1.  When  this  ratio 
approaches  unity  the  inclination  of  tho  P  front  goes  toward  a  90°; 
if  one  considers  the  target  point  A  in  Fig.  (3-8);  the  first  signal  will 
originate  from  a  point  F  whose  horizontal  distance  (because  of  90°) 
is  quite  large  compared  to  the  depth.  Signals  from  other  intermediate 
points,  say  F',  will  arrive  a  very  short  tine  later.  Obviously,  if  the 
horizontal  distance  AF  is  large,  the  pressures  at  F  and  at  intermediate 
points  T*  will  differ  radically,  and  a  v  toady-  state  approximation  becomes 
invalid  except  for  very  small  values  of  z,  and  for  a  very  short  time  after 
arrival  of  the  first  signal . 

One  can  use  the  above  consideration  to  Judgo  the  applicability  in  cases 
of  small  ratios,  V/cp  <  1*5#  by  comparing  pressure  and  velocity  above  the 
point  considered  with  these  at  point  F.  In  the  ease  V/cp  -  1.5  shown  in 
FlS*(3-U  the  angle  of  tho  P-vave  is  a  -  Ul°50*;  considering  a  target  A 
at  a  depth  z  »■  500  feet,  the  horizontal  component  of  the  distance  is 
AF  »  z  sin  a  »  333  feet.  Variations  of  pressure  and  velocity  V  in  such  a 
distance  are  such  that  the  result  for  tills  particular  depth  may  still  bo 
considered  a  crude  approximation. 

Whenever  the  steady- state  solution  becomes  unsatisfactory  the  early 
part  of  the  signal  could  be  determined  by  using  idealized  curves  for  the 
response  due  to  a  concentrated  load,  similar  to  the  one  shown  in  Fig.  (3-5)* 
Further  work  in  this  direction  appears  desirable  because  tho  range  where 
th$  steady  state  approximation  is  unsuitable  is  still  of  practical  interest. 
However,  even  by  using  idealized  response  curves  the  numerical  integrations 
required  are  likely  to  be  cumbersgmo. 
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De torn! nation  of  Rayleigh  Wave  Effecto. 

The  analysis  contained  in  [4]  and  [5]  va*  carried  further*  [9]j  and 
expressions  for  otrcooco  and  displacements  duo  to  the  ourfaco  waves  gen- 
crated  by  a  concentrated  load  aro  given  in  Appendix  3.  The  expressions 
are  in  closed  form,  and  represent  the  major  portion  of  the  total  response 
for  largo  ratios  r/x  and  in  the  vicinity  of  t  -  -  cR/c.  E*B** 

C/vjir^  oDiCr»s  found  will  give  a  very  good  approximation  of  the  portion  near 
T  -  0.9  cf  the  curve  for  r/x  »  40  in  Fig.  (3-3)* 

As  example  the  stress  a zz  due  to  surface  waves  is  given  in  Fig.  (3-9) * 

It  is  seen  that  effects  occur  only  during  a  ahort  period  of  time*  of  the 
order  two  to  three  times  x/c,  where  x  is  the  depth  of  the  point  considered. 

To  obtain  the  Rayleigh  effects  from  the  surface  forces  due  to  the  pres¬ 
sure  history  in  Fig.  (l-l),  ft  formal  double  integral  is  derived  in  Appendix  B. 
Its  evaluation  for  typical  locations  required  the  use  of  an  IbW  704.  The 

analysis  indicates  that  the  stresses  a  are  larger  than  the  other  stresses. 

rr 

9  °xx  9  °rx  Buc^  importance  of  Rayleigh  effects  generally  may 

be  Judged  by  considering  orr#  Figs.  (3-10)  and  (3-H)  show  the  computed 
time  history  of  the  stresses  o ^  at  a  number  of  depths  In  two  iocutions. 

Fig*  (3*10)  preoonts  the  stresses  at  r  -  2000  ft.  for  a  medium  having 
V  -  1/4,  c  ..  10,000  ft/eec  (Cp  ^  17,000  ft/oec).  Fig. (3-ll)  gives  similar 
results  for  r  «  3200  ft.  for  a  medium  with  v  -  l/4,  but  having  c  -  6000  ft/sec. 
(cp  ~  10,000  ft/oec).  The  decrease  of  Rayleigh  stresses  with  depth  is  quite 
evident. 

Comparisons  of  stresses  due  to  Rayleigh  effects  and  due  to  the  effects 
of  P-  and  8 -wave a  are  shown  in  Figs.  (3-12)  and  (3-13)*  The  former  shows 

°nd  £or  r  »  2000  ft.  at  a  depth  of  100  ft.  for  the  same  medium  for  which 
PiR*  (3-10)  applies.  In  this  location  the  velocity  V  of  the  shock  in  air  is 
larger  than  17,000  ft/uec.,  and  the  "steady  state  solution,"  Appendix  A, 
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could  be  used  to  compute  the  pressure  waves  arriving  ahead  of  the  curfacc 
vavco.  Tho  time  t  uccd  in  the  graphs  ic  counted  with  respect  to  the  in¬ 
stant  of  the  explosion  at  r  o  0,  The  two  spikes  in  ozz  represent  the 
arrival  of  the  P-  and  S-vaves,  respectively.  Sinilnr  spikes  occur  in  o  * 
except  that  the  one  due  to  the  S-vavo  ic  not  in  the  cane  direction  a a  the 
one  due  to  the  P-vnve,  ouch  that  the  stresses  remin  high  only  for  a 
3hort  tine.  Fig.  (3-12)  shews  that  in  this  location  the  pea];  values  of  the 
signal  from  P-  and  S-vaves  arc  larger  than  Fayleigh  stresses,  but  the  dura¬ 
tion  of  high  Hayloigh  o tree boo  in  o^r  is  longer.  Depending  on  the  response 
tine  of  the  target  Rayleigh  otrcoocc  nay  or  nay  not  be  ignorablc.  At  greater 
depth  of  course  the  effects  of  the  P-  and  S-vaves  do  not  change  materially , 
vhile  Rayleigh  effects  decrease  rapidly,  such  that  the  latter  vlll  become 
immaterial. 

Fig. (3-13) ohovu  a  similar  comparison  for  r  «  3200  ft.  at  100-ft.  depth 
for  the  medium  uoed  in  connection  with  Fig.  (3-U). 

Additional  numerical  results  ore  presented  in  Appendix  B.  Attention  is 
drawn  to  the  fact  that  the  formulae  developed  in  Appendix  B  canno-  be  applied 
for  points  on  the  surface,  z  =  0,  because  negative  powc  “  of  z  occur.  As  a 
consequence  computational  difficult ies  would  arise  for  very  small  depth, 
and  only  results  for  points  at  a  depth  of  25  feet  or  more  were  obtained. 

Ifce  trend  of  these  recults  as  a  function  of  the  depth  is  such  that 

(*) 

rather  large  stresses  o ^  and  Vqq  can  be  expected  near  the  surface. 

The  existence  of  such  high  stresses  in  the  elastic  analysis  make  it  in¬ 
applicable  to  actual  situations  near  the  surface  z  -►  0  3uch  that  the  addi¬ 
tional  effort  required  to  derive  an  appropriate  analysis  for  shallow  depth 
was  not  expended. 

{*)  The  vertical  stress  s  and  the  shear  stress  c  or  the  free  sur- 

zz  rz 

face  oust  vanish. 
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The  occurrence  of  high  Rayleigh  stresses  near  the  surface  in  itself  is 
a  limitation  to  the  applicability  of  the  linear  elastic  theory  to  actual 
solutions.  Referring  os  example  to  the  situation  in  Fig.  (3-12),  at  a  dis¬ 
tance  of  r  ■  2000  ft.,  the  air  pressures  acting  at  the  surface  arc  such 
that  the  use  of  linear  clootie  theory  is  reasonable  because  the  otresoeo 
obtained  froa  the  steady-state  solutions  are  sufficiently  lov.  The  stresses 
due  to  Rayleigh  vavco  (vhich  occur  at  a  later  tine)  at  thio  horizontal 
diotonce  at  the  depth  z~100  ft.  will  olco  be  snail  enough  to  be  in  the 
elastic  range,  but,  ao  the  sane  con  not  be  noid  near  the  ourfaco,  the 
Rayleigh  ctrcooeo  everywhere  can  not  be  considered  reliable;  the  values 
computed  by  the  theory  arc  presumably  much  too  large.  Combining  thio  con¬ 
sideration  with  the  effect  of  the  crater,  discussed  in  Section  6b.,  the 
Rayleigh  stresses  in  Figs*  (3-10) to  (3-1^  should  be  considered  as  upper  bounds 
only. 

d)  Effects  outside  the  supcrscismic  range,  at  depth  z  «  r* 

It  is  noted  that  even  outside  tho  superscicnic  range  the  effects  at 
shallow  depth  z  «r  can  be  divided  into  a  signal  due  to  the  air  pressure 
passing  over  the  target  point,  and  a  Rayleigh  signal.  The  latter  can  be 
dot eroined  by  the  method  previously  outlined,  which  does  not  depend  on  a 
particular  location  of  the  target  point.  Appendix  B  contains  also  results 
for  oubsoioaic  locations.’ 

To  obtain  the  effect  of  the  air  pressure  passing  over  the  target  point 

one  can  try  to  use  the  steady  state  approximation.  The  solution  to  the 

steady  state  problem  is  available  [0]  but  it  is  somewhat  more  complex  than 

(*) 

in  the  oupersciomic  range.  There  are  actually  two  ranges,  depending  on 
whethor  V  >  c  or  V  <  c.  Further,  a  steady  state  solution  vill  again  be 

useless  if  the  velocity  is  close  to  Cp.  Another  trouble  spot  has  been 

w  Potentially  useful  expressions  for  the  steady- state  solutions  due 
to  exponentially  decaying  pressures  are  contained  in  [?6].  They 
might  be  utilized  to  tabulate  the  solutions. 
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polnted  out  by  Miles  [10].  Tho  steady-state  oolutlon  indicatea  infinite 
otrooooo ;  find  bocotjco  inapplicable  i.C  V  -  vhoro  io  tho  Rayleigh 
vave  speed.  A  formulation  for  a  surface  preooure  vith  varying  speed  vac 
obtained  in  Reference  [10]  to  overcome  thic  difficulty.  No  numerical 
rooultc  oro  available  at  preoont,  5ho  coabination  of  tho  steady  otato 
eolution  and  of  the  Rayleigh  vave  contributions  gives  therefore  only  e,  par¬ 
tial  answer  in  thic  range,  ^  It  should  be  stressed  that  the  cteody-otato 
solutions  outside  the  supersoiomic  range  indicat  attenuation  of  tho  response 
with  depth;  this  vould  indicate— at  least  in  the  theoretical  clastic  medium— 
that  the  surface  vave  effects  outside  the  superoeismic  range  may  remain 
significant  at  larger  depth  than  super  seismic  range. 


(*)  No  caaputationo  to  check  the  accuracy  of  the  steady-state  approxima¬ 
tions  v similar  to  Fig.  2-5)  appear  to  have  been  made  by  anyone. 

Such  a  check  may  be  desirable  because  qualitative  comiderations 
lead  to  the  suspicion  that  this  approximation  outside  the  supers eis- 
rdc  range  might  be  poor. 
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APPLICATION  OF  DYNAMIC  RECIPROCAL  THEOREM 

wp  VERTICAL  SURFACE  DISPLACEMENT  DUE  TO  VERTICAL  FORCE  H(t)  (PEKERIS) 

wc  VERTICAL  SURFACE  DISPLACEMENT  DUE  TO  HORIZONTAL  FORCE  (CHAO) 

w  VERTICAL  DISPLACEMENT') 

f  DUE  TO  VERTICAL  FORCE  ON  SURFACE 
v  HORIZONTAL  DISPLACEMENT  J 


H(t) 


wp(t) 


V(t)  =  Wc(t) 
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G.3-2 


RAYLEIGH  SURFACE  WAVE 
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DISPLACEMENTS  OBTAINED  BY  PEKERtS  AND  '.1FSON 
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IDEAUZED  w(h,t)  FOR 
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F  I  G.  3-10  RALEIGH  STRESSES  AT  R  3  2000  FT 
(C  *  I  0000  FT/SEC) 
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4.  LAYERED  MEDIA. 

In  npito  of  the  fact  that  layering,  or  gradual  changes  of  the  nedium 
vlth  dopth  will  have  pronounced  effects,  only  a  few  applicable  solutions 
are  available.  This  io  due  to  the  complexity  of  the  problems  Of  layered 
media,  The  steady  state  coco  of  a  pressure  pulse  progressing  with  constant 
speed  V  on  the  surface  of  a  fluid  layer  above  a  fluid  half  space,  Fig*  (4*1) 
has  been  considered  [3]  without  restriction  as  to  the  value  of  the  velocity  V. 

The  equivalent  problem  for  a  layered  elastic  solid  has  been  solved  (11] 
in  an  elementary  manner,  but  only  in  the  superseiomic  range,  i.e.  when  V  is 
larger  than  tho  sound  velocities  in  either  medium.  Fig. (4-2) shews  stresses 
o  due  to  a  stop  wave  of  constant  pressure,  when  the  velocities  of  c  (2) 

LL  p 

(2) 

and  c  in  the  lower  medium  are  twice  those  in  the  upper  one.  It  is  seen 
that  the  vertical  stresses  ozz  both  at  pointo  in  the  upper  and  at  points 
in  the  lover  medium  are  increased  by  the  reflections  from  the  interface.  If 
the  applied  pressure  on  the  surface  has  a  decaying  history,  a  case  which  is 
treated  in  til],  the  reflections  may  be  less  important,  because  they  are 
added  to  a  reduced  stress  from  the  direct  sigLal.  However,  it  must  be 
emphasised  that  layering  always  increases  the  vertical  stress es  if  the 
velocity  of  sound  increases  with  depth.  Reduction  of  pressure,  indicating 
a  protection,  occurs  only  in  the  unusual  case  when  the  sound  velocity  de¬ 
creases  with  depth.  The  reverse  applies  for  vertical  displacements  and 
accelerations;  they  decrease  due  to  layering  if  the  sound  velocity  increases 
with  depth.  The  approach  used  in  Reference  [11]  can  be  extended  to  any 
number  of  layers,  as  long  as  the  velocity  V  is  larger  than  the  sound  veloci¬ 
ties  in  any  layer. 

A  problem  involving  a  lsyer  of  fluid  over  an  elastic  half-space  has 
been  treated  theoretically  and  experimentally  In  Reference  [12],  where  the 


recponse  froa  a  point  explosion  in  tho  lnyor  ha*  boon  obtained,  Pig.  (^3). 
The  caphotis  in  thie  paper  vno  on  vaveo  propagated  on  tbo  interface.  It 
la  experimentally  poosiblo  to  ncaoure  preeaurei  in  a  flnid  reliably,  and 
good  agreement  betveen  theory  and  experiment  vn»  obtained. 
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5.  ONE-DIMENSIONAL  WAVE  PROPAGATION  IN  A  NON-LINEAR  ELASTIC  MEDIUM* 

In  Sections (2)  and  (3) it  could  be  demonstrated  that  for  the  cases  of 
linear  behavior  of  fluids  and  elastic  solids  the  assumption  of  one -dimensional 
vave  propagation  gives  good  results  for  the  early  part  of  the  response  of 
pointB  belov  ground  zero,  r  ~  0.  This  conclusion  can  be  drawn  qualitatively 
for  any  other  medium  if  the  velocity  of  wave  propagation  is  much  smaller 
than  the  velocity  of  the  blast  wave  on  the  surface. 

The  early  pressures  near  r  2*  0  in  actual  situations  are  eo  high  that 
non-linear  effects  must  occur  in  any  conceivable  medium.  It  is  therefore 
very  important  that  the  effect  of  pressure  distributions  of  the  type  given 
in  Fig.  (l-l)  can  be  determined  no  a  problem  in  one -dimensional  wave  propaga¬ 
tion,  because  such  a  problem- -even  if  non-linear— can  be  handled  by  numerical 
methods.  As  a  matter  or  fact  some  useful  general  conclusions  can  ue  drawn 
for  a  general  class  of  non-linear  elastic  cases. 

Consider  a  rod-like  element  of  the  half-space,  Fig.  (5-1).  The  "rod" 
has  unit  cross-sectional  «rea  and  extends  from  the  surface  z  «  0  to  z  - 
In  the  present  situation  the  horizontal  strains  must  vanish  such  tl  at  the 
area  of  the  rod  does  not  change.  Let  Fig.  (5-2)  be  the  appropriate  stress- 
strain  diagram  of  vertical  stress  a  versus  vertical  strain  e.  For  suffi¬ 
ciently  email  stress,  say  0  <  o^;  the  law  may  be  assumed  to  be  linear. 

Let  a  decaying  pressure  pulse  p(t)  «  a  (0,  t)  be  applied  to  the  rod,  such 
that  the  peak  of  p(t)  is  larger  t  ion  o^. 

The  response  for  the  case  described  above  could  be  computed  numerically 
for  any  given  stress-strain  diagram.  However,  an  equivalent  problem  has 
been  considered  in  1 3^1 #  8ec.  II  for  a  bi-linear  (plastic-elastic)  medium. 

The  medium  has  the  stress-strain  diagram,  shown  in  Fig*  (5“3)>  which  is  a 
special  case  of  the  one  indicated  in  Fig.  (5-2).  Fig. (5-4)  shows  the  stress 
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response  for  a  a  example,  which  can  be  considered  typical.  The  example 
was  computed  for  an  input  pro o sure  at  a  «  0  whicn  h &A  a  peak  of  l8J',000 
pel,  while  at  points  below  the  surface,  t  •  500  feet  and  z  *  700  feet, 
much  smedler  peak  stresses  were  found.  On  the  other  hand,  stresses  at 
t  *  100  ms  01  more  after  arrival  at  the  lower  levels  were  larger  than 
the  corresponding  input  stresses.  (The  reduction  of  peak  stresses  must 
be  compensatev  by  some  increase  of  stresses  at  later  time  to  satisfy  the 
lav  of  conservation  of  momentum.) 

From  the  above  example  one  concludes  that  in  the  non-linear  case  a 
reduction  in  peak  pressure  with  depth  (attenuation)  is  to  be  expected. 

This  attenuation  will  however  only  continue  until  the  peak  pressure  has 
dropped  to  the  value  or,  defining  linear  behavior  in  the  stress-strain 
diagram  fig.  (5-2). 

As  long  as  the  applied  load  ia  such  that  the  assumption  of  one- 

dimensional  wave  propagation  is  reasonable  one  must  expect  peak  pressures 

of  the  magnitude  of  the  stress  defining  linear  behavior.  Some  idea 

of  the  level  of  0^  can  be  gained  from  the  simple  tests  on  granites  described 

in  Appendix  F.  It  is  seen  that  the  stress-strain  relations  are  fairly 

linear,  up  to  fracture  on  the  15,000  to  20,000  lb.  level.  In  the  case  of 

plane  wave  propagation  where  transverse  strains  and  therefore  fractures 

similar  to  the  tests  are  inhibited,  it  is  to  be  expected  that  a  nearly 

,  2 

linear  relation  will  hold  up  to  even  higher  stresses,  oL  >  20,000  lo/in* 
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6.  USE  OF  THEORETICAL  RESULTS  FOR  PREDICTIONS  IN  ROCK. 

When  cocaparing  the  complexity  of  an  actual  cituation  vith  the  relative 
simplicity  of  the  problems  analysed,  it  should  be  clear  that  theoretical 
solutions  can  not  give  a  complete  picture  of  stresses  everywhere,  but  can— 
at  best— permit  a  rough  estimate  of  the  stress  field  in  certain  locations 
at  certain  times.  In  view  of  the  scarcity  of  experimental  information, 
the  theoretical  solutions  must  be  exploited  by  searching  through  qualitative 
reasoning  for  situations  where  they 

Before  discussing  specific  locations,  it  is  convenient  to  list  major 
effects  not  appearing  at  present  in  the  theoretical  solutions.  Concerning 
the  effect  of  the  explosion,  only  the  air  pressure  on  the  surface  hQ3  been 
considered,  ignoring  not  only  heat  and  radiation  effects,  but  also  effects 
due  to  the  direct  impact  of  the  casing  at  the  point  of  explosion.  Concern* 
ing  behavior  of  the  material,  high  pressure  levels  requiring  on  equation  of 
state  are  not  within  the  range  of  this  study,  which  eliminates  the  crater 
rogion  from  further  consideration.  On  the  stress  level  considered  here, 
materials  might  not  behave  elastically  but  in  a  viscous  or  plastic  manner, 
or  exhibit  other  disBipative  mechanisms.  The  fact  that  only  small  strains 
are  considered  should  also  be  mentioned  but  is,  on  this  stress  level,  not 
likely  to  be  serious,  as  the  strain  will  not  exceed  a  few  per  cont. 

a)  Target  locations  below  the  crater. 

Consider  a  target  below  the  crater;  it  may  be  at  twice  the  depth  of  the 
crater,  Fig.  (6-1).  Without  knowing  the  pressure  history  exactly,  we  know 
that  very  high  pressures  are  propagated  downward  over  an  area  exceeding  the 
diameter  of  the  crater. 

While  the  formation  of  a  crater  is  evidence  that  the  situation  is  not 
really  a  "plane  wave"  situation,  the  fact  that  the  crater  is  a  flat  bowl 
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dissipative  effects  are  not  likely  to  be  very  Important  (which  does  not 
mean  that  solutions  Including  such  effects  should  not  be  sought).  Esti- 
metes  of  the  first  signal  based  on  steady  state  solutions,  particularly  In 
the  super seismic  range,  should  therefore  be  of  real  significance.  The 
pronounced  difference  according  to  theory  between  ouperoeismlc  and  sub- 
seiamic  ranges,  the  former  having  no  or  negligible  attenuation  of  peak 
pressure  with  depth,  while  the  latter  has  pronounced  attenuation,  can  be 
expected  to  be  a  real  fact. 

The  situation  concerning  the  second  signal  is  much  more  uncertain, 
because  it  is  due  to  forces  acting  in  the  crater  region  where  the  elastic 
theory  can  not  possibly  apply.  Hie  elastic  theory  predicts  Rayleigh  (surface) 
waves  as  the  major  signal  to  be  received;  now,  the  existence  of  such  waves 

requires  on  interaction  between  direct  and  ehear  stresses,  such  that  they 
do  not  occur  Id  a  fluid.  At  high  pressures  in  the  crater  region  the 
material  is  acting  like  a  fluid,  and  one  is  led  to  conclude  that  the  Incep¬ 
tion  of  surface  waves  will  be  inhibited,  and  that  the  elastic  theory  will 
grossly  exaggerate  surface  waves.  For  this  reason  the  Rayleigh  stresses 
computed  in  Appendix  B  must  be  considered  too  high,  presumably  very  much 
too  high.  } 

There  is  a  second  reaoon  why  the  Rayleigh  stresses  in  actuality  must 
be  less  than  found  from  the  elastic  theory.  If  the  stresses  found  wore 
otherwise  correct,  those  near  the  surface  are  larger  than  rocks  are  likely 
to  withstand.  It  would  not  matter  if  such  high  stresses  were  nearly  hydro¬ 
static,  but  the  state  of  stress  in  surface  waves  ie  inherently  deviatory 

(*)  When  arguing  that  the  actual  stresses  due  to  surface  waves  are  smaller 
than  given  by  elastic  theory,  one  should  not  he  misled  into  the  con¬ 
clusion  that  they  can  always  and  everywhere  be  ignored*  Zt  should  be 
remembered  that  at  very  large  distances  seismologies!  experience  proves 
beyond  doubt  that  Rayleigh  waves  ore  the  dominant  effect.  There  must 
therefore  be  a  transition  range  beyond  which  surface  waves  can  not  ba 
ignored. 
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(i.e.  having  considerable  shear  deformation).  The  region  near  the  sur¬ 
face  ought  therefore  to  be  analyzed  as  an  elastic- plant ic  medium  which 

(*) 

would  lead  to  lover  stresses. 

Pending  further  study  the  Rayleigh  stresses  obtained  in  Appendix  B 
can  be  considered  only  as  an  upper  bound.  Such  a  bound  io  useful  because 
whenever  comparison  of  the  firot  part  of  tho  response  indicat  ‘  at  it 
exceeds  the  Rayleigh  stresses  in  importance!  the  latter  con  be  safely 
ignored.  E.g.,  in  the  Buporseismic  range  this  permit-  determination  of  a 
depth  bolcv  which  Rayleigh  stresses  are  cortain  to  bo  small  compared  to 
the  early  reoponoe. 


(*)  Any  other  dissipative  effect  would  also  lead  to  lower  stresses  due  to 
surface  waves.  Such  waves  in  a  visco-elastic  medium  have  recently 
been  considered  (15). 


CENTER  OF  BLAST 


O  TARGET 


•61. 


7.  WAVE  PROPAGATION  IN  GRANULAR  MEDIA. 

The  behavior  of  certain  solids  under  very  high  intensity,  confined 
or  hydrostatic  pressure  is  characterized  by  a  non-linear  stress-strain 
relation  of  positive  curvature  which,  upon  unloading,  shows  large  perman¬ 
ent  strains.  This  response  is  manifested  at  dirferont  pressure  levels  in 
various  materials:  at  low  pressures  in  foam-rubber,  in  the  100  to  16,000 
psi  range  in  granular  soils  and  presumably  at  much  higher  pressures  in 
porous  or  solid  rocks. 

A  typical  stress-strain  diagram  of  this  type  is  shown  on  Figure  (7  l); 
this  diagram  can  be  idealized  by  various  linear  diagrams  shown  in  Figures 
(7-2),  (7  -3)  and  (7-*0. 

Figures  (7  2)  (a),(b)  and  (c)  show  a  class  of  materials  known  as  lock¬ 
ing  media  (a),  or  materials  of  limited  compressibility,  (b)  and  (c). 

Ihese  materials  have  the  property  of  becoming  eo  sent  lolly  incompressible 

or  rigid  upon  compression  when  a  critical  strain  c  is  reached*  At  that 

c 

instant  a  sudden  locking  or  compaction  of  the  medium  occurs,  and  upon 
unloading  the  permanent  strain  remains. 

If  a  certain  amount  of  residual  elasticity  exists  after  compaction 
or  locking,  the  material  can  be  represented  by  the  bi-linear  diagram  shown 
on  Figure  (7-3)*  If  the  irreversible  compaction  proceeds  gradually,  tho 
stress-strain  relation  is  that  of  ~lgure  (7-*0*  This  material  acts  elasti¬ 
cally  on  loading  and  unloads  with  a  permanent  set  as  shown  in  the  figure. 

Rather  extensive  literature  is  available  on  the  static  behavior  of 
locking  media,  [2b]  and  [25]  and  plane  wave  propagation  in  both  locking 
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media,  {Figure  7-2  (a),  (b)  and  (c)],  and  in  the  type  of  material  shown  on 
Figure  (7-*0  which  hao  also  been  investigated  ([26]  and  [27))*  I*Sore  recently, 
problems  of  vnvc  propagation  in  bi-lincar  materials  of  the  type  shown  in 
Figure  (7-3)  have  been  solved.  Spherical  waves  in  locking  media  are  ex¬ 
plored  in  References  [293,  (30)  and  (3lh 

tost  of  these  investigations  find  application  in  researches  connected 
with  the  prediction  of  ground  shock  effects  due  to  thermonuclear  blasts  y 
[32]  and  t33l;or  are  concerned  with  the  detection  of  underground  nuclear 
tests,  (30]  and  [ 3ll *  An  application  to  the  design  of  invulnerable,  deep 
underground  shelters  is  considered  in  References  [18]  and  [I9l>  and  a 
general  discussion  of  these  topics  is  found  in  Reference  [3*0* 

Some  of  the  relevant  results,  such  as  shock* velocity  and  the  attenua¬ 
tion  of  peak  pressure  intensities  with  depth  may  be  summarized  a a  follows: 

If  the  surface  of  a  semi-infinite  locking  medium  of  the  type  shown 
in  Figure  (7-2a)  is  subjected  to  a  decaying  pressure  p(t)  the  shock 
velocity  i  of  the  front  located  at  a  depth  z  is  given  in  [26]  and  [27]  by: 

z  -  Ii2p«c/I{i)dt]-1''2  (7-1) 

where  p  is  the  density  of  the  medium,  and  I,  the  impulse  of  the  applied 
pressure,  is 

I  ■  /  p(t)dt  (7-2) 

The  amplitude  of  the  peak  stress  0  is  related  to  the  depth  z  by: 

z 


(7-3) 
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For  a  medium  vhich  ohovo  a  gradually  irreversible  compaction  (Figure  7-^)> 
tho  above  quoted- reference  give#  the  propagation  velocity  as: 

c0  *  (E/p)"^2  ■  constant  (7-*0 

vheia  E  io  tho  elastic  modulus  of  the  medium.  The  attenuation  of  the  peak 
stress  vith  depth  follows  the  lav: 

°z  "  ^o2'1  (7-5) 

Numerical  evaluation  of  the  otreoces  and  particle  velocities  in  a  locking 
nwdiun,  (Figure  7-2a)#  subjected  to  high  intensity  blast  pressures  dis¬ 
closes  the  following  significant  rosults  [l8]. 

(a)  The  intensity  of  peak  otrcoooo  (and  of  tho  particle  velocities) 
in  lockirg  nedia  attenuates  very  rapidly  with  depth.  (See  Figure  7-5). 

(b)  The  attenuation  of  these  quantities  it  to  $,  surprising  degree 
intuitive  to  variations  of  the  parameters  of  the  medium.  For  example, 
at  a  high  yield  surface  burst  near  ground  zero,  the  intensity  of  stresses 
beyond  a  depth  of  500  feet  differs  only  slightly  for  materials  vhich  have 

a  vide  range  of  critical  strains  Uc  *  0.02  to  0.10  as  shown  on  Figure  7-5)* 

(c)  The  existence  of  residual  elasticity  has  an  equally  slight 
influence:  A  locking  medium  with  a  critical  strain  of  cq  *  0.02  iE  compared 
vith  a  bi-linear  medium  of  identical  critical  strain  but  vith  a  residual 
elastic  modulus  of  E  -  10^  psi  in  Figure  (7-6)  and  ve  note  that  the  attenua¬ 
tion  of  peak  stresses  beyond  about  *^00  feet  of  depth  is  not  significantly 
affected.  These  phenomena  have  been  found  to  exist  for  a  rather  vide  range 
of  variations  for  both  parameters  (cc  and  E)  and  are  the  subject  of  detailed 
investigations  in  Reference  [28). 

These  conclusions  are  of  interest  because  they  show  that  the  essential 
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ond  significant  attenuation  of  bloat  pressures  can  be  anticipated  in  a 
large  variety  of  dry  granular  coils.  Consequently,  the  free  field  para¬ 
meters  should  be  predictable  if  a  limited  amount  of  information  is  avail¬ 
able  regarding  subsurface  conditions. 

As  mentioned  previously  the  theory  of  spherical  vave  phenomena  in 
locking  media  has  been  applied  to  study  the  effects  of  underground  explo¬ 
sions.  In  [29)  the  effect  of  an  exploding  point  maos  is  considered  with 
a  viev  tovardo  exploring  the  effect  of  the  so-called  ‘’direct  ground- shock”, 
(i.e.  the  exploding  casing  of  the  weapon)  in  a  granular  soil. 

It  is  found  that  the  spherical  shock  at  a  radius  R  propagates  vith 
the  velocity: 

r  -  r0(i  +  (~)3ru  (7-6) 

m 

where  the  kinetic  energy  of  the  mass  M  of  the  casing  is: 

1  *2 

and  the  quantities: 


of  a  surface  burst,  as  an  early  time  approximation  of  the  direct  ground 
shock  effect* 


i 
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I 


DxmAcrnoN  effects 

8.  imODUCTION 

Thermonuclear  explosions  on  the  surface  of  an  elastic  half- space  pro¬ 
duce  a  complex  stress  pattern  in  the  medium.  In  the  moot  general  ca3e, 
the  circularly- oymmctric  surface  prcoouroo  from  the  explosion,  vhich  expand 
radially  outward  from  Ground  Zero  in  space  and  decay  in  time,  produco  a 
stress  pattern  at  pointB  in  the  medium  in  vhich  the  two  principal  stresses 
change  in  magnitude  and  in  direction  with  time.  The  diffraction  of  the 
stress  field  by  a  cavity  containing  a  hardened  underground  installation 
must  be  considered  with  a  view  towards  obtaining  the  following  information: 

1)  Stroooos  producod  by  the  dynamic  loading  of  the  cavity  with  a  view 
towards  determining  the  strength  of  the  cavity  ao  a  whole.  In  addition, 
the  velocities,  displacements  and  accelerations  of  points  on  and  near  the 
cavity  boundaries  are  required  for  the  determination  and  evaluation  of  shock 
effects  and  for  the  establishment  of  failure  criteria  for  the  cavity  and 
its  contents. 

2)  Shock  spectra  for  a)  the  total  accelerations  imparted  to  the  con¬ 
tents  of  the  cavity,  and  b)  the  relative  displacements  Of  the  contents  of 
the  cavity  (relative  to  the  cavity  boundaries)  when  the  cavity  is  onvolopod 
by  the  stress  waves  produced  by  the  explosion.  The  contents  of  the  cavity 
may  be  shock  mounted  and  the  shock  opectra  are  required  for  an  optimisation 
of  the  design  of  the  installation. 

Ifce  solution  to  the  general  problem  of  the  diffraction  of  the  rather 
cocaplicated  stress  field  by  the  cavity  in  an  elastic  medium  con  be  constructed 
for  caeep  of  interest  by  the  superposition  of  the  resultB  obtained  from  some 
basic' less  complicated  problomo  which  are  considered  in  Part  II*  specifically, 
the  diffraction  Of  the  stress  field  produced  by  P  (dllatationol)  and  8  (shear) 
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vavao  respectively  vith  plane  wave  fronts  ere  considered.  In  both  cones 
the  general  colution  io  obtained  for  a  step  pressure  in  time  which  can 
be  used  od  influence  coefficients  to  give  the  results  for  waves  vith 
arbitrarily  time  varying  pressures  by  means  of  Duhamel  integrals. 

The  solution  of  these  problems  (stresses,  velocities  and  displace¬ 
ments)  can  be  used  to  construct  solutions  for  the  general  case  in  which 
ta©  incoming  stress  history  has  a  more  complicated  character. 

Moreover,  they  arc  of  direct  use  and  importance  in  the  superooiomic 
range  (V>c^)  in  which  the  loading  on  the  cavity  has  been  shown  to  consist 
of  P  and  S  waves  with  plane  fronts  carrying  pressure  components  which  decay 
in  time  (see  Appendix  A).  For  this  range,  the  solutions  represent  an 
approximate  answer  to  the  actual  physical  .problem  and  can  be  used  directly. 

The  above  theory  has  been  determined  with  the  assumption  that  the 
cavity  containing  the  installation  io  not  lined.  Practically  speaking,  the 
problem  of  the  cavity  lining— structural  or  anti- spalling- -is  quite  important 
in  the  design  of  these  installations  and  must  be  studied.  A  method  of  attack 
in  which  the  results  obtained  for  unlined  cavities  are  u3cd  as  influence 
coefficients  in  integral  equations  for  the  corresponding  solutions  to  cavi¬ 
ties  vith  linings  is  also  developed  in  this  report  [Section  111. 

Section  (9)  gives  the  results  for  stresses,  velocities  and  displacements 
produced  at  the  boundary  of  the  cavity  ao  it  is  enveloped  by  piano  P  and  B 
wave e.  Curves  are  presented  for  the  case  in  which  the  pressure  inputs  in 
the  waves  are  step  functions  in  time,  as  well  as  for  waves  vith  particular 
time  decaying  pressures  produced  by  a  hypothetical  20  MT  surface  burst.  In 
addition,  the  superposition  of  the  incoming  P  and  S  waves  is  shown  for  a 
particular  tunnel  location. 
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Section  (10)  contain®  the  ohock  spectra  curvea  for  the  total  accelera¬ 
tion  and  relative  displacement  of  installations  which  are  ohock  mounted  in 
the  cavity.  It  is  assumed  that  the  equipment  can  be  mounted  in  several  dif¬ 
ferent  vivye  and  shock  spectra  are  presented  for  several  casco  of  interest* 
Heretofore,  shock  epectra  used  in  this  type  of  design  were  computed  from 
free  field  input  pressures,  and  the  effects  of  the  diffraction  of  the  waves 
by  the  cavity  were  neglected.  The  shock  opectra  of  Section  (10)  include  the 
diffraction  effects  and  are  consequently  more  appropriate;  they  will  now 
supersede  the  free  field  spectra  and  should  be  of  considerable  use  in 
obtaining  more  accurate  design  data  for  the  ohock  mounting  of  the  cavity 
contents. 

Section  (11)  presents  an  outline  of  the  theory  by  which  the  various 
results  obtained  for  an  unlineu  cavity  con  be  used  to  determine  the  corres¬ 
ponding  reoultB  for  a  cavity  with  a  given  elastic  lining.  Numerical  solu¬ 
tions  of  the  resulting  integral  equations  are  not  yet  available,  but  should 
be  completed  in  the  near  future  for  several  cases  of  interest. 

Section  (12)  presents  bu-uc  gciigiul  conclusions  relative  to  the  "c suite 
obtained  in  Sections  (9)-(ll). 

Following  the  format  of  Fart  I,  the  theoretical  formulations  of  the 
work  are  presented  in  Appendices  C,  D,  E,  while  the  results  are  reported  in 
Sections  (9),  (10),  (ll),  (12).  The  results  of  some  tests  to  determine  the 
elastic  properties  or  granite  und*  static  loading  are  given  in  Appendix  (F). 
The  following  recurring  symbols  are  adhered  to  in  Part  II  of  this  report. 

Cp  Velocity  of  pressure  (p)  waves  in  a  linear  elastic  solid. 

C  Velocity  of  shear  (S)  waves  in  a  linear  elastic  solid. 

8 

The  velocity  c  appearing  in  the  abscissa  of  the  various  curves  in  Part  II 
refers  to  the  velocity  of  the  incoming  wave  as  shown  in  the  figure  or  noted. 
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9.  STRESSES,  VELOCITIES,  AHP  DISPLACEMENTS  PRODUCED  IN  AM  ELASTIC  MEDIUM 
BY  Tig  DIFFRACTION  OF  PLANE  "P"  AND  MSM  WAVES  BY  A  CYLINDRICAL  CAVITY 

An  infinitely  long  cylindrical  cavity  in  on  infinite  elastic  homogeneous 
and  ic.  tropic  medium  is  acted  on  by  a  plane  check  vnvo  whose  vave  front  io 
parallel  to  the  axis  of  the  cavity.  The  shock  vave  propagates  through  tho 
medium  vith  a  constant  velocity  cp  ("Pn  vave)  and  envelops  the  cavity 
(Fig.  9.1).  For  generality,  it  io  assumed  that  the  direct  stress  components 
carried  by  the  vave  are  eU(t)  and  ccrU(t)  ,  which  are  respectively 
parallel  and  perpendicular  to  the  direction  of  vave  propagation. 

The  solution  of  this  problom  for  tho  inputs  aU(t),  («  »  0)  or  «oU(t) 
only,  may  be  used  to  construct  the  solutions  to  problems  in  which  the  free 
field  has  a  more  general  nature.  For  tho  auperoeiomic  range,  in  which  a 
component  of  the  input  io  actually  a  plane  wave,  the  value  of  «  for  a  plane 
wave  front  roust  be  token  as  c  »  -l/3* 

Preliminary  results  for  the  stress  field  produced  at  the  boundary  of 
the  cavity  by  the  incoming  MP*'  wave  havo  boon  roportod  in  (17),  [l8],  tl9l, 
while  the  complete  plane  strain  solution  to  the  problem  h*"  been  given  in 
l 20 1 ♦  The  preueut  papsr  extends  these  results  to  includo  both  P  and  S  waves* 
In  addition,  results  arc  prooontod  for  the  volocity  and  displacement  fields 
produced  in  the  medium  at  and  near  the  cavity  boundary.  Tho  theoretical 
development  for  the  velocities  and  displacements  are  presented  in  Appendix  D. 

Although  the  true  physical  problem  under  consideration  require s  a  semi- 
infinite  medium  vith  a  roughly  p*anc  boundary  at  z  *  0,  it  io  oaoily  shown 
(Ref,  18,  Pg.  36]  that  in  tho  range  of  practical  interest,  the  plane  boundary 
essentially  has  no  effect  during  times  of  interest  and  the  reoulto  obtained 
by  considering  tho  raodium  to  bo  infinita  will  bo  satisfactory  for  tunnels 
where  depth  "D"  is  greater  than  4  to  5  times  the  radius  "a".  For  such 
installations,  the  major  effects  such  as  maximum  hoop  stresses,  occur  at 
the  boundary  of  the  cavity  at  times  considerably  shorter  than  the  arrival 
tiiaa  of  the  relief  vave  from  tho  piano  surface  at  z  ■  0. 
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a)  Hoop  stress  produced  at  the  boundary  of  a  cavity  [20 ]* 

Figs.  (9-2)- (9-8)  shew  the  hoop  stresses  o00  produced  at  points  on 
the  cavity  boundary  by  various  stress  waves  uhich  envelop  tho  cavity.  The 
stresses  at  tvo  pointo  are  presented:  0  -  0°,  the  point  at  vhich  the  shock 
front  first  hits  the  tunnel;  and  0  «  90°,  the  point  at  vhich  the  maximum 
stress  concentration  is  expocted.  In  determining  0^  by  the  analytical  pro¬ 
cedure  given  in  [20],  only  the  terms  n  -  0,  1,  2  plus  tho  free  field  stress 
components  were  included;  tho  analytical  Justification  for  terminating  the 
aeries  after  the  n  -  2  term  is  presented  in  tlQ],  Pg.  Ul  and  [20 ). 

Fig*  (9-2)  shows  the  hoop  stress  at  r  »  a,  vhich  is  produced  by  an 
incoming  plane  shock  wave  vith  a  step  pressure  distribution  in  time.  Ike 
pressure  component  cU(t)  at  right  angles  to  the  direction  of  propagation  of 
the  wave  is  taken  equal  to  zero.  The  results  obtained  at  long  times  must 
approach  the  veil  known  static  solution  for  a  cylindrical  hole  in  a  uniaxial 
pressure  field,  i.e.,  a  stress  amplification  at  0  »  90°  of  3  (compression), 
vhile  the  stress  at  0  -  0°  is  a  tension  equal  in  magnitude  to  unity.  It  is 
seen  that  the  stress  concentration  is  amplified  by  the  dynamic  loading  in 
the  ratio  of  3*28  to  3  at  ^  transit  times,  at  9  *  9°°  and  l,l6  to  1.00  at 
0  »  0°.  It  may  be  noted  that  for  a  step  pressure  vave,  the  problem  is  essen¬ 
tially  of  a  quasi- static  nature;  since  the  dynamic  amplifications  are  small 
percentagewise.  It  is  also  of  considerable  interest  to  note  that  maximum 
stresses  are  produced  at  about  four  to  five  transit  times  of  the  shock  vave 
across  the  cavity. 

Fig.  (9-3)  shows  the  hoop  stress  at  r  *  a  vhich  is  produced  by  an 
incoming  plane  shock  vave  vhich  has  no  pressure  in  the  direction  of  propaga¬ 
tion  of  the  vave  and  a  pressure  <U(t)  in  the  direction  perpendicular  to  the 
propagation  direction  of  the  vave. 
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Tbc  results  of  Ti g.  (9-2)  and  (9-3)  may  be  used  to  construct,  by  suit¬ 
able  integrations,  solutions  for  cases  in  vhich  tho  free  field  pressure 
is  of  a  more  general  nature. 

Pig.  (9-*0  shows  the  stress  o0g  at  0  «  0°  and  0  •>  90°  for  a  plane 
step  shock  wave  with  a  plane  front  progressing  with  a  velocity  cp,  i.e.  tho 
standard  P  wave  of  linear  elasticity  theory.  The  requirement  of  a  plane 
front  necessitates  a  tranoveroo  pressure  component  with  c  -  -l/3  t*o.#  -  1/3  U(t) 
(Fig*  9-1)  as  veil  as  a  pressure  U(t)  in  tho  direction  of  wave  propagation.  * 

For  the  dynamic  loading,  tho  otreoc  amplifications  are  from  2.667  to  2*92  at 
6  0  90°  (compression)  and  from  0.00  to +  0.11  at  0  «  0°  (tension).  These 
results  may  be  used  directly  to  obtain  meaningful  results  in  the  super- 
Geiomic  range  in  which  the  vaveo  vhich  envelop  the  cavity  are  essentially 
plane  waveB. 

Tho  hoop  stress  o^g  produced  by  a  step  shock  wave  may  be  used  as  an 
s 

influence  function  to  determine  the  corresponding  stress  produced  by  a  wave 
with  a  time  varying  pressure,  P(t),  by  the  Duhamel  integral'  '(Fig.  9-5). 


°ee  3  p0°ee 


+  °ge  (*-»)  dT 


(9-1) 


dT  ”0 

To  illustrate  the  above  procedure,  Figs.  (9-6)- (9-7 )  present  the  hoop 
stresses  produced  by  plane  waves  with  decaying  pressure  time  histories. 

The  pressure  time  histories  of  the  P  waves  which  envelop  the  cavity  are 


token  as  those  produced  on  the  s  jrface  by  a  20  MT  surface  burst  at  the 
65OO  psi  (2l*00  ft.  from  G.£.)  contour  and  the  2000  poi  (3200  ft.  from  G.Z.) 


(*)  Because,  only  a  limited  number  of  modes  have  been  used  in  determining 
the  stress  due  to  the  plane  Bhock  wave,  the  Duhamel  integral  should  not  be 
applied  when  the  pressure-time  history,  p(t),  contains  significant  high  fre¬ 
quency  components  with  time  constants  of  less  than  l/2  of  one  transit 
time. 
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contour.  These  pressure  tine  curves  have  been  given  by  Brodo  [l]»  The 
cavity  is  considered  to  have  a  radius  of  17*5  ft*  in  an  elastic  medium  in 
vhich  cp  -  17,300  ft/sec  and  cq  «  10,000  ft/scc. 

Fig.  (9-6)  gives  the  stress  at  0  «  0°  and  0  **  90°  for  a  plane 

?  wave  with  a  Brode  pressure  input  vith  PQ  =  2000  psi.  Due  to  the  decay 
of  the  pressure  vith  time,  the  Duh&mel  integral  chows  a  maximum  compres¬ 
sive  stress  of  4500  psi  at  0  «  90°  while  the  tension  at  6  -  0°  is  290  poi. 
Pig.  (9-7)  shows  the  corresponding  curves  for  a  pressure  input  vith 
PQ  -  65OO  psi.  The  maximum  compressive  stress  at  90°  is  12,300  psi,  while 
the  tension  at  0  «  0°  is  930  psi.  It  is  of  interest  to  note  that  the 
amplification  of  the  compressive  stress  over  PQ  at  $  ■  90°  was  2.25  for 
the  2000  psi  loading  and  I.89  for  the  6500  psi  loading,  as  compared  to 
2.92  for  a  step  pressure  loading.  Since  the  pressure-time  decay  for  the 
2000  psi  wave  was  considerably  slower  than  that  for  the  65OO  psi  wave,  the 
hoop  stresses  are  closer  to  those  far  a  step  wave. 

Pig.  (9-8)  shows  the  hoop  stress  at  5  -  o°  and  cF  *  45°  for  &  plane 
step  shear  wave  with  &  plane  front  and  a  constant  velocity  c,  i.e.  the 
standard  8  wave  of  linear  elasticity  theory.  At  long  times,  the  etross 
approaches  the  static  solution  for  a  cylindrical  hole  in  a  bi-axial  pressure 
field  vhich  produces  the  shear  stress  distribution  vhich  is  carried  by  the 
wave;  i.e.  a  stress  amplification  at  0  -  45°  of  4  (compression).  It  is  seen 
that  the  stress  is  amplified  by  the  dynamic  loading  in  the  ratio  of  4.37  to 
4  at  45°  and  t  «  4  transit  times,  and  .112  to  0  at  9  *  0°.  It  may  be  noted 
that  as  in  the  case  of  the  P  wave  enveloping  the  cavity,  the  problem  fez-  a 
stop  pressure  input  is  essentially  quasi-static  in  nature  (since  the 
dynamic  amplifications  are  small  percentagewise)  and  that  maximum  hoop 
stresses  are  produced  about  two  to- three  transit  times  of  the  S  wave  across 
the  cavity.  The  theoretical  results  for  the  diffraction  of  the  S  wave  by 
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tfeo  cavity  can  to  derived  fraa  the  theory  for  the  P  wave;  this  dcvelopcacnt 
la  given  in  Appendix  (c). 

In  the  superseismic  range  (V  >  cp),  the  pressure  waves  that  propagate 
through  the  medium  ore  standard  P  and  S  vavoo  vrith  plane  front o  and  propa¬ 
gation  velocities  and  c^roopoctivaly.  Those  vavoo  can  bo  conflidorod 
to  have  the  same  prooourc -tine  hiotory  ao  the  surface  prooourc,  and  ampli¬ 
tude  factors  vhich  are  given  by  Fig.  (A-2),  Appendix  A.  The  P  and  S  waves 

are  inclined  to  the  ourfaco  by  the  angles  a  »  oin~^  V  and  3  *»  sin**'  V  , 

Z  cfl 

P  fl 

respectively.  If  the  cavity  is  located  at  a  sufficient  depth  [see  Sec.  6] 
so  that  surface  vave  effects  are  negligible,  one  can  obtain  a  realistic  pic¬ 
ture  of  tho  stress  build-up  at  tho  cavity  boundary  by  cui>orimpooing  the 
incoming  P  and  8  vavoo  vith  thoir  proper  amplitude  fact  or  o  and  timo  dolayo 
(due  to  different  angles  of  inclination  and  different  velocities  of  propa¬ 
gation). 

As  an  example  of  tho  ouperposition  procedure  described  above,  consider 

a  cavity  of  radius  a  »  17.5  ft*  in  an  elastic  medium  in  vhich  cp  n  17# 300 

ft/eec  and  cq»  10,000  ft/oec.  Hie  cavity  io  located  at  a  depth  of  £00  ft* 

at  the  65OO  psi  surface  pressure  contour  from  a  20  MT  surface  burst,  i.o., 

at  a  distance  of  about  2000  ft.  from  Ground  Zero.  Tho  air  shock  velocity 

at  this  distance  io  about  22,000  ft/oec j  this  correspond^  to  a  Mack  number 
V 

M  «  ~  »  I.27.  From  Appendix  A,  the  P  vave  has  on  amplitude  factor  of 
cp  0 

O.97  and  is  inclined  at  an  angle  ,  «  £1  -5U*  to  the  horizontal  surface, 
while  the  S  vave  has  an  amplitude  factor  of  0.53  an d  an  angle  of  inclination 
of  3  !*  26°-48'.  [See  Fig49->*)]  Fig.  (9-10)  shows  the  hoop  stress 
vhich  is  produced  at  various  points  on  the  cavity  boundary  by  the  super¬ 
position  of  the  P  and  S  waves  with  their  appropriate  time  factors.  Due  to 
the  difference  in  propagation  velocities  of  the  two  waves,  tho  delay  time 
between  the  arrival  of  the  P  and  8  waves  is  considerably  in  excess  of  5  transit 


times. 
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Consequently,  by  the  time  the  shear  vavo  arrives  at  the  cavity,  the 
hoop  stresses  produced  by  the  P  vave  have  decayed  sufficiently  bo  that  the 
peak  of  the  sum  of  the  otreosoo  produced  by  the  tvo  waves  (  0  ^20°)  la 
at  111  lover  than  the  peak  hoop  atreoo  at  ©  «  90°  produced  by  the  P  vavo  at 
early  ti x&a,  before  the  t-rrivul  of  the  S  vavo.  It  nay  bo  noted  that  the 
maximum  peak  stress  of  31,930  psi  (compression)  i6  in  excess  of  the  allow¬ 
able  stress  vhich  a  granite  rock  could  reasonably  be  expected  to  carry; 
therefore,  one  ouat  consider  lover  pressure  contours  for  the  case  of  unlined 
cavities . 

Fig.  (9-H)  shows  the  stress  onQ  at  0  -  0°  and  0  «  90°  for  a  plane 
P  wave  viUk  a  constant  velocity  of  propagation,  cp  »  10,393  ft /sec# 
tc  «  6000  ft/oec.]  and  a  pressure  tine  history  corresponding  to  the  surface 
pressure  at  the  3200  ft  [2000  psi]  contour.  The  air  blast  at  this  con¬ 
tour  has  a  velocity  of  12,100  ft  •*  3ec. ,  .thus  making  the  caoe  superseismic 
with  a  Mach  number  of  1.17  and  consequently  an  amplitude  factor  of  1.0  in 
the  p  vave.  '  The  cavity  has  a  radius,  a  =  17*5  ft*  The  maximum  hoop 
stress  is  produced  at  #  •*  90°  and  ban  a  magnitude  of  4500  psi.  This  is  a 
reasonable  stress  for  a  good  granite  to  carry  and  one  may  conclude  that  a 
reasonable  pressure  limit  for  feasible  underground  installations  will  bo 
in  the  neignborhood  of  the  2000  pei  contour  for  the  hypothetical  20  MP 
surface  burst  under  consideration. 

Comparison  with  the  results  of  Ref.  [£!]. 

The  analysis  presented  in  Sec.  (9)  is  based  on  the  mode  approach  re¬ 
quiring  an  expansion  of  the  stresses  as  a  function  of  0  in  a  Fourier  series. 
In  the  response  to  a  step  wave,  the  higher  terms  of  the  series  could  be 
ignored,  and  only  the  terms  up  to  and  including  n  a  2  were  used.  It  is 


(*)  See  Appendix  A,  Fig.  (A-2).  As  pointed  out,  the  steady-state  solution 
for  a  Mach  number  so  close  to  unity  may  not  he  very  accurate  and  the 
actual  b,ppp  stresses  at  the  points  in  question  will  be  somewhat  lover* 


interesting  to  note  that  moults  aro  available  from  an  alternative  analysis, 
vhich  permit  checking  of  tho  present  results  over  part  of  the  circumference 
of  the  tunnel. 

She  alternative  analysis  [21  ]  io  also  approximate;  it  applie*  only 
for  limited  valuao  of  tho  angle  |e|<  60°.  Ihc  hoop  stress  OqQ  taken  froa 
Pig.  (7)  Of  Ref*  [21]  is  plotted  in  Pie*  (9-12).  This  figure  also  shews 

Qqq 

the  peak  valuoo  of  thiG  otreoo  -j-  found  by  the  present  analysis'  9  for 
$  “  22.5°;  4$°,  60°  and  90°*  It  io  seen  that  the  values  for  22.5°  and  45° 
agree  quite  veil;  for  60°  the  present  analysis  gives  nearly  10^  more  otreoo 
than  Ref.  (2l).  As  tho  latter  becomes  more  and  more  approximate  as  Q 
increases,  the  present  theory  io  believed  to  be  the  better  value.  Finally, 
for  0  o  90°,  Ref.  [2l]  gives  no  result  at  all,  but  states  an  extrapolation 
that  the  maximum  value  of  is  about  2.  The  present  analysis  indicates 
that  this  extrapolation  vus  not  Justified,  the  actual  value  being  about  3* 

0 QQ 

It  is  easily  seen  that  the  maximum  stress  -j-  for  a  step  vave  can  not 
possibly  bo  2,  but  must  bo  larger  than  2.  Ibis  can  bo  concluded  from  the 
fact  that  the  dynamic  peak  stress  mot  be  larger,  or  at  least  equal  to  the 

0 QQ 

static  peak  stress.  Tho  latter  i°  "  2*52  for  v  **  0.33  and  2.67  for 

V  »  0.25* 

To  find  the  maximum  amplification  of  stress  for  a  step  wave,  or  for 

a  decoying  vave  vith  a  decay  constant  of  more  than  one  half  of  a  transit 

time,  it  io  therefore  appropriate  to  use  the  present  analysis  rather  than 

the  extrapolation  in  [21].  To  find  the  response  due  to  pro  a  sure  waves  which 

have  very  fast  decays,  the  present  mode  approach  is  not  suitable,  but  the 

approach  of  [21]  could  be  used  instead,  numerical  verk  for  such  a  purpose, 

comiiderably  beyond  that  presented  in  [ell,  uould  be  required/**^  Moreover, 
(*)  Note  that  Ref.  [21]  and  the  present  analysis  use  a  dif ferent  volue  of 
Poisson's  ratio,  i*e,  v=  0,33  and  v  *  0.25,  respectively. 

(**)An  estimate  of  the  pressure  at  early  times,  given  in  [l 0],  pp.  42-10, 
is  also  available  and  may  somet lines  suffice. 
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it  should  he  noted  that  situatloiiS  in  vhich  the  pressure  decay  io  sufficiently 
font  such  that  the  present  analgia  becomes  unsuitable,  will  occur  only 
very  clooo  t<?  Ground  Zero,  i.o.,  in  a  region  of  extremely  high  preosuro. 

b)  Velocities  and  displacements  produced  at  points  in  the  clastic  medium. 

The  displacement  and  velocity  components  v,  v,  w  and  v  vhich  arc  pro¬ 
duced  at  points  in  the  elastic  medium  by  a  plane  step  shock  wave  are  given 

by  Eq.  (34)- (37) i  Appendix  D.  Numerical  results  arc  presented  for  those 
quantities  at  points  on  tho  cavity  boundary,  r  a  a.  As  in  the  case  of  the 
stresses,  tho  velocities  and  displacements  produced  by  a  step  shock  wave 
may  be  used  as  influence  functions  to  determine  tho  corresponding  quanti¬ 
ties  produced  by  a  vave  vith  a  time  varying  pressure  P(t),  by  a  Duhorael 
integral  similar  to  Eq.  (9-l).  In  turn,  these  reoultc  can  be  used  as 
input  functions  for  tho  determination  of  shock  effects  vhich  are  imparted 
to  shock  mounted  installations  vlthin  the  cavity*  They  win  be  used  in 
the  determination  of  the  acceleration  and  displacement  shock  spectra  vhich 
are  presented  in  Secxion  (10). 

Fig*  (9-13;-(9-l4)  show  tho  dioplacomonto  and  volocitios  produced  by 
a  plane  step  shock  wave  at  the  boundary  points,  e  «  0°  and  0  »  900.  hx, 
these  points,  the  ooticn  is  purely  radial,  i.e,  v  and  v  =»  0.  More  compre¬ 
hensive  results  arc  given  in  Fig.  (7)- (8),  Appendix  1),  in  vhich  the  veloci¬ 
ties  £r  and  v  are  given  for  other  points  on  the  cavity  boundary. 

The  rigid  body  translational  motion  of  the  cavity  boundary  in  vhich 
the  cavity  maintains  its  cylindrical  shape  and  translates  in  the  direction  of 
the  incoming  step  shock  vave,  Is  extracted  from  the  total  motion  in  Sec.  (6), 

Appendix  D.  Fig.  (9-15)  shows  tho  rigid  body  displacement,  velocity  and 

(*)  See  Appendix  (E),  Part  lll-a,-  for  the  formal  dcrivation  of  expres-  * 
sions  for  tho  velocities  and  displacements  produced  at  points  on  the 
chvlty  boundary  by  waves  vith  Brode  pressure  inputs. 


acceleration  of  the  cavity  boundary,  under  the  step  pressure  shock  vavt 
loading*  These  results  may  also  je  used  as  influence  coefficients  in 
IXihscel  integrals  to  obtain  co  >  n  ng  umtitiea  produced  by  shock 
vaves  vith  time  varying  pressure  histories.  Results  of  this  type  for 
vaves  vith  Erode  pressure  inputs >  will  be  used  as  input  functions  for 
the  determination  of  acceleration  and  displacement  shock  spectra  in 
Section  (lO). 
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10*  SHOCK  SPECTRA  FOR  INSOLATIONS  IN  CYLINDRICAL  CAVITIES  IN  EU5TIC  MEDIA  . 

In  many  casco,  the  variouo  installations  vhich  arc  placed  in  under¬ 
ground  cavities  will  contain  components  vhich  arc  quite  shock  sensitive. 
Consequently ,  they  vill  require  special  mountings  vhich  are  capable  of 
absorbing  the  shock  effects  produced  by  the  pressure  loading  on  the  cavity. 
Shock  spectra  for  accelerations  and  displacements  are  frequently  utilized 
in  the  design  of  shock  mounted  equipment.  The  present  section  presents 
typical  shock  spectra  for  a)  the  motion  of  various  points  on  the  cavity 
boundary,  and  b)  the  rigid  body  (mean)  motion  of  the  cavity  as  a  vhole. 

The  analytical  formulation  for  the  development  of  shock  spectra  has 
been  presented  in  detail  in  Appendix  (E).  Two  particularly  usefu*  types 
of  spectra  arc  considered:  l)  spectra  for  the  peak  absolute  acceleration 
vhich  is  imparted  to  a  shock  mounted  installation  in  the  cavity;  2)  spectra 
for  the  peak  relative  displacement  of  the  installation  in  the  cavity,  with 
respect  to  the  motion  of  points  on  the  cavity  boundary  to  vhich  it  is 
attached.  The  former  spectra  give  the  acceleration  design  requirements 
for  a  given  shock  mounted  inczallatic  the  latter  give  clearance  require¬ 
ments  for  the  mounting  of  the  installation  vithin  the  cavity. 

The  shock  spectra  developed  in  the  present  Section  and  in  Appendix 
(E)  include  the  effects  due  to  the  diffraction  of  the  shock  wave  by  the 
cavity.  Heretofore,  shock  spectra  vhich  wore  used  in  the  Assign  of  under¬ 
ground  installations  vere  computed  from  free  field  input  pressures  only, 
that  is  from  the  readings  of  a  pressure  gage  at  a  point  in  a  medium  vith 
no  tunnel,  and  the  effects  of  the  diffraction  of  the  shock  vave  by  the 
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be  determined  from  the  analysis  in  Appendix  (D);  the  procedure  for  obtain¬ 
ing  the  total  displacement  under  the  combined  ?  and  S  wave  loading,  is 
identical  to  that  described  in  Section  (9)  for  the  determination  of  the 
hoop  stress. 

(a)  Shoes  Spectra  for  the  Rigid  Body  Motion  of  the  Cavity* 

Consider  a  mass  M  which  is  mounted  to  the  boundary  of  a  cavity  by 
either  of  the  two  methods  shown  in  Figure  (10-1).  In  the  first  case,  the 
mass  is  connected  to  the  boundary  at  many  points;  each  mounting  is  charac¬ 
terized  as  being  a  linear  spring.  For  this  case,  the  relative  motion  of 
any  two  points  on  the  cavity  boundary  will  be  averaged  out  by  the  springs 
which  are  connected  to  the  other  support  points.  Consequently,  the  mass 
M  will  react  to  the  average  motion  of  the  cavity  which  is  a  rigid  body 
translation  of  the  cavity  in  the  direction  of  the  shock  wave  propagation. 

In  the  second  case,  the  masB  is  connected  as  rigidly  as  possible  to  a  stiff 
structural  lining.  Again,  the  mass  will  react  to  the  rigid  body  transla¬ 
tion  of  tne  cavity  as  a  whole. 

Figure  (10-2)  shows  acceleration  and  relative  displacement  spectra 
for  the  rigid  body  motion  of  a  cavity  of  radiuo  r  -  17*5  feet  in  an  elas¬ 
tic  medium  where  cp  »  17,300  feet/oec*  The  cavity  ie  subjected  to  a  P 
wave  carrying  a  Brode  pressure  input  with  ?o  *  6500  poi  for  a  20  MT  surface 
burst.  For  comparison,  the  corresponding  acceleration  spectrum  from  the 
free  field  pressures  only  is  shown  as  a  dotted  line*  It  is  noted  that 
for  higher  frequency  components,  the  free  field  accelerations  are  much 
larger  than  those  in  the  spectrum  from  the  analysis  of  Appendix  (E  )which 
includes  the  diffraction  effects;  the  difference  increases  rapidly  as  the 
frequency  increases*  The  difference  io  duo  to  two  cau&oss  l)  the 


bo  determined  from  the  analysis  in  Appendix  (D);  the  procedure  for  obtain¬ 
ing  the  total  displacement  under  the  combined  P  and  5  vave  loading,  is 
identical  to  that  described  in  Section  (9)  for  the  determination  of  the 
hoop  stress* 

(a)  Shook  Spectra  for  the  Rigid  Body  Motion  of  the  Cavity  * 

Consider  a  mass  M  vhich  is  mounted  to  the  boundary  of  a  cavity  by 
either  of  the  two  methods  shown  in  Figure  (10-1).  In  the  first  case,  the 
mass  is  connected  to  the  boundary  at  many  points;  each  mounting  is  charac¬ 
terized  as  being  a  linear  spring.  For  this  case,  the  relative  motion  of 
any  two  points  on  the  cavity  boundary  will  be  averaged  out  by  the  springe 
vhich  are  connected  to  the  other  support  points*  Consequently,  the  mass 
M  will  react  to  the  average  motion  of  the  cavity  vhich  is  a  rigid  body 
translation  of  the  cavity  in  the  direction  of  the  shock  vave  propagation. 

In  the  second  case,  the  mass  is  connected  as  rigidly  as  possible  to  a  stiff 
structural  lining.  Again,  the  mass  will  react  to  the  rigid  body  transla¬ 
tion  of  the  cavity  as  a  whole. 

Figure  (10-2)  shows  acceleration  and  relative  displacement  spectra 
for  the  rigid  body  motion  of  a  cavity  of  radius  r  **  17*5  £ect  ln  clas¬ 
tic  medium  vnere  c  17,300  fect/ccc.  The  cavity  is  subjected  to  a  V 
vave  carrying  a  Brode  pressure  input  with  PQ  =  6500  poi  for  a  20  MT  surface 
burst.  For  comparison,  viie  corresponding  acceleration  spectrum  from  the 
free  field  pressures  only  is  shown  as  a  dotted  line.  It  is  noted  that 
for  higher  frequency  components,  the  free  field  accelerations  are  much 
larger  than  those  in  the  spectrum  from  the  analysis  of  Appendix (E  )which 
includes  the  diffraction  effects;  the  difference  increases  rapidly  as  the 
frequency  increases.  The  difference  is  due  to  two  causes t  l)  tho 
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diffraction  effect  vhich  is  maximum  at  relatively  early  times,  and  thuo 
affects  the  high  frequency  portion  of  the  shock  spectrum,  and  2)  the  fact 
that  the  rigid  body  notion  of  the  cavity  produces  the  shock  spectrum  for 
the  configurations  in  Figure  (10-1)  rather  than  the  notion  of  a  particle 
subjected  to  the  full  free  field  pressure.  Uha  spoctrua  for  the  relative 
displacement  of  the  shock  mounted  mass  M  vith  respect  to  the  notion  of 
the  cavity  boundary  is  also  shown  in  the  figure. 

Figure  (10-3)  oho vo  similar  results  for  the  sane  cavity  in  the  elastic 
medium.  In  this  case,  the  cavity  is  subjected  to  a  P  wave  vith  a  Brode 
pressure  input  vith  PQ  *  2000  psi.  Again,  a  20  WT  surface  buret  is  consi¬ 
dered.  This  case  is  felt  to  be  more  realistic  vith  respect  to  an  actual 
installation,  since  it  has  bean  shown  in  Section  (9)  that  shelters  in 
rock  at  the  2000  psi  contour  night  be  theoretically  feasible. 

(b )  Shock  Spectra  for  the  Motion  of  Individual  Points  on  the  Cavity 
Boundary . 

Consider  a  maos  Mvhichis  mounted  by  a  linear  spring  to  a  point  on 
the  boundary  of  the  cavity.  Figure  (10-4)  shows  the  shock  spectra  for 
the  motion  at  the  boundary  point  0  *=  0°.  Again,  the  input  wave  is  pro¬ 
duced  by  a  20  MT  surface  buret,  and  the  peak  pressure  Po  is  6500  psi.  It 
is  seen  that  the  free  field  accelerations  in  thic  case  are  considerably 
lover  than  the  accelerations  from  the  analysis  of  Appendix  (E),  particularly 
at  hign  frequencies.  This  is  to  be  expected,  oince  the  presence  of  the 
free  cavity  boundary  gives  rise  initially  to  a  doubling  of  the  particle 
velocity  at  the  boundary  point  0  =  0°  and  consequently,  to  an  increase 
in  the  displacement  at  the  boundary  points  at  early  times.  Hence,  the 
high  frequency  accelerations  at  the  point  0  =  0°  are  considerably  higher 
than  those  obtained  from  the  free  field  analysis  and  the  more  accurate 
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anolysis  which  includes  the  diffraction  effects  must  he  used.  The  spectrum 
for  the  relative  displacement  of  the  shock  mounted  mass  M  with  respect 
to  the  motion  of  the  cavity  boundary  at  0  »  0°  is  also  shown  in  the  figure. 

Figure  (10-5)  shows  the  shock  spectra  for  the  motion  of  the  boundary 
point  0  ■  180°  undU-r  the  action  of  a  P  wave  from  a  20  KT  surface  burst 
with  a  peak  pressure,  Po  *  6500  psi.  In  this  case,  the  free  field  accelera¬ 
tion  spectrum  is  considerably  higher  than  the  spectrum  including  the 
diffraction  effects;  the  diffraction  actually  decreases  the  displacement 
input  and  consequently  the  high  frequency  accelerations. 

The  motion  of  the  cavity  bounda;*y  at  tho  points  1  *»  0°  and  0  ■<  l80° 
that  have  been  considered  thus  far  is  purely  radial;  the  tangential  conpon- 
ent  v  is  equal  to  zero  at  all  tines'  .  Hence,  the  linear  oscillator 
concept  can  be  applied  in  determining  shock  spectra  for  equipment  in  the 
cavity  which  is  connected  to  thccc  points.  However,  other  points  on  the 
cavity  boundary  will  undergo  both  radial  and  tangential  displacements 
and  shock  spectra  for  both  directions  are  required.  Consequently,  the 
shock  mounting  of  equipment  to  such  points  must  be  capable  of  absorbing 
accelerations  in  both  the  radial  and  tangential  directions,  i.e.  ar 
essentially  "two-way"  shock  mounting  is  required. 

(c)  Conclusions 

It  is  obvious  that  the  most  favorable  conditions  for  shock  effects 
will  be  encountered  If  the  equipment  Is  shock  mounted  in  the  cavity  such 

that  it  will  react  to  the  average  motion  of  the  cavity  (rigid  body  motion) 
rather  than  to  the  motion  of  individual  points  on  the  boundary  of  the 
tunnel*  Such  mountings  may  be  difficult  to  obtain  practically.  However 


(*)  It  should  be  noted  that  this  is  only  true  if  the  S  wave  effects  can 
be  neglected. 
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if  ouch  a  mounting  con  ho  achieved,  the  equipment  will  he  subjected  to 
ohock  accelerations  vhich  are  considerably  smaller  than  those  given  by 
the  free  field  shock  cpectra. 

In  the  case  of  equipment  vhich  in  attached  to  a  point  on  the  tunnel 
circumference,  the  spectra  allowing  for  diffraction  differs  considerably 
from  the  free  field  spectra.  The  peak  accelerations  for  a  point  9  «  0° 
are  considerably  increased  by  the  diffraction.  However,  the  opposite 
situation  prevails  at  9  *  180°,  in  vhich  the  shock  accelerations  includ¬ 
ing  diffraction  effects  sire  lower  than  the  "free  field'Values. 

In  each  case,  the  design  accelerations  for  high  frequency  components 
dro  greatly  influenced  by  the  diffraction  effects  of  the  wave  by  the 
cavity.  Consequently,  it  is  felt  that  the  more  accurate  theory  of 
Appendix  (E),  vhich  includes  these  effects,  should  be  used  in  the  develop¬ 
ment  of  shock  spectra  and  design  criteria. 

It  should  be  noted  that  the  shock  spectra  which  have  been  presented 
in  this  Section  as  examples,  art  based  on  an  input  of  a  sharp  fronted  pressure 
(P)  vave  with  a  aero  rise  tin*.  For  the  caae  of  a  wave  with  a  finite  ria# 
time,  the  intensity  of  the  high  frequency  portions  of  the  ehock  spectra  would 


be  substantially  decreased. 
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11.  Diffraction  of  Preopuro  Waves  by  a  Cylindrical  Elastically  Lined 

Cavity  in  an  Elastic  Medium. 

The  stresses  and  motions  produced  by  the  diffraction  of  pressure 
waves  by  an  unlined  cylindrical  cavity  in  an  elastic  medium  have  been 
ctudicd  in  Reference  [2o]  and  Appendiceo  (C)  and  (d).  It  is  of  interest 
to  evaluate  the  effect  of  the  use  of  a  structural  elastic  lining  for  the 
cavity  boundary.  Particularly,  the  effect  of  such  a  lining  on  the  strength 
of  the  cavity  as  a  whole  vnd  on  the  boundary  displacements  which  are  the 
input  functions  for  the  determination  of  shock  spectra  for  installations 
within  the  cavity,  must  be  studied. 

The  problem  considered  is  that  of  an  elastic  shell  in  an  infinite 
elastic  medium  under  the  action  of  a  plane  step  pressure  wave  which  travels 
through  the  medium  and  envelops  the  cavity.  A  method  of  solution  which 
utilises  the  solution  for  the  corresponding  unlined  cavity  problem  as 
influence  coefficients  can  be  used  to  obtain  the  corresponding  results 
for  the  elastically  lined  cavity.  Such  a  method  has  already  been  applied 
successfully  for  a  radially  symmetric  problem  in  an  acoustic  medium  [22]. 

The  generalization  of  this  method  for  the  present  problem  is  considerably 
more  conplex  and  represents  a  major  computational  effort.  This  effort  is 
currently  under  way;  consequently,  only  an  outline  of  the  method  of  solution 
is  presented  in  this  Section. 

An  infinitely  long  elastically  lined  cylindrical  cavity  in  an  infinite 
elastic  homogeneous  and  isotropic  medium  is  acted  on  by  a  plane  pressure 
wave  whose  front  is  parallel  to  the  axis  of  the  cavity*  The  shock  wave 
propagates  through  *he  medium  with  a  constant  velocity  and  envelope 
the  cavity  [Fig#  (11-1)]. 
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a)  Determination  of  the  equations  of  Motion  for  the  Shell. 


The  response  of  the  elastic  lining*  ^  to  the  transverse  shock  wave 
is  studied  by  considering  the  shell  in  vacuo  as  a  separate  eiructuro 
responding  to  the  dynamic  forces  exerted  by  the  surrounding  elastic 
medium  and  by  the  incoming  shock  wave.  Using  the  modes  or  free  vibration 
of  the  shell  in  vacuo  as  generalined  coordinate*,  its  response  can  be 
expanded  in  terms  of  the  infinite  number  of  these  modes*  The  mode  shapes 
and  the  corresponding  frequencies  that  are  required  have  been  developed 
in  Reference  [23] •  Let  the  displacement  of  the  shell  be  written  in  terms 
of  the  generalised  coordinates  q^#  and  q^^  (n  f  0) : 


v(e,t)  »  q0(t)  +y[qn(t)  +  ^(t)] 


cos  n 9 


r*l 


(li.i) 


sin  n$ 


(U.2) 


The  quantities  q^^  arc  the  coordinates  for  the  primarily  inex^nsional 
motions  of  the  shell;  the  quantities  *qo  and  ^  aro  the  coordinates  for  the 

t  * 

primarily  extensionel  motions. 


The  coefficient  is  given  by  the  relation 


?or  properties  of  practical  interest  and  for  the  lower  modes  (n  <  5)# 

there  is  little  coupling  between  the  bending  and  extenslonal  effects 
o 

because  a  A  »  L 

n  "  . — 

7  The  lining  ia  now  essentially  an  elastic  shell  in  an  infinite  elastic 
atdium* 


(U-3) 


*\  o  o 

For  such  cases,  *  n  and  the  frequencies"'  ^  and  a£  are  given  by 


—  p  2  p 

Mg .  mi £ljL 

n  ma4(n2+l) 

(11-4) 

“n  “  \  (n2+1) 

(H-5) 

where  E  -  - ^  ,  I  ana  A  are  Tihe  moment  of  inertia  of  the  cross  section 

1-V 

of  the  shell  with  respect  to  a  principal  axis  at  right  angles  to  the  pl*ne 
of  the  ring  and  the  croso  sectional  area  of  the  ring  respectively,  and  m 
is  the  maos  per  unit  area  of  the  elastic  lining. 


The  equations  of  motion  for  the  shell  may  be  written  in  terms  of  the 
generalised  coordinates  q^  and  q^  j 

n  -  1,  2,  3,  ...  (11*6) 


2  «n 


%  +  “n*n  "  5“ 
n 


90  o  Q— 

\+*XmT 


n  m  0,  1,  2,  3# 


(11-7) 


where  and  are  generalised  forces  and  mn  and  mQ  are  the  generalized 

masses :  r 


mn-“(1  +  i) 


(U-8) 


\  “  n(l  +  <**)  1 


(H-9) 


The  generalized  forces  Qfl(t)  and^(t)  will  be  evaluated  later  in  this  Section. 


*) 


For  cases  in  which  the  coupling  between  i-he  bending  and  the  extensional 
motions  is  significant,  the  frequencies  must  be  determined  from  the 
egressions  given  by  Eq.  (A-8)  of  Reference  [23] . 
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b)  The  Response  of  the  Boundary  of  an  Unllned  Cylindrical  Cavity  ** 
Determination  of  Influence  Coefficients . 

The  displacements  at  the  boundary  of  an  unlined  cavity  to  a  etep 

pressure  wave  have  been  evaluated  in  Appendix  (p)  of  this  paper. 

Essentially,  the  results  were  obtained  by  the  addition  of  an  effect 

produced  by  the  free  field  component  of  the  incoming  shock  wave  plus 

the  effect*  produced  by  the  application  of  corrective  boundary  traction* 

which  are  required  to  produce  a  traction  free  surface  at  the  cavity 

boundary.  For  the  present  purpose,  the  free  field  components  must  be 

expanded  into  a  Fourier  series  in  0,  and  added  to  the  corresponding 

components  from  the  corrective  tractions.  Finally,  the  total  motion  of 

points  on  the  boundary  of  the  unlined  cavity  is  obtained  in  the  form 
00 

v  cos  nO  (11-10) 

n*Q 

00 

v  -^0n(t)  sin  n0  (ll-ll) 

n»l 

Two  additional  coefficients  are  required  from  the  unlined  cavity 

problem.  They  are  respectively!  l)  the  boundary  displacements  vA(t)  coo  n9 

and  vA(t)  sin  n 9  produced  by  the  applied  boundary  tractions  »  U(t)  cot  a© 
a  rrn 

°r9  *  0  (11-2) 1#  2)  the  boundary  displacements  w0(t)  cos  n0  and 

n 

Vrs('C')  »ia  uv  produced  by  the  applied  boundary  tractions  *  0; 

*  .  *rn 

a^Q  •  U(t)  sin  n0  [Fig.  (ll-3)]*  Thess  quantities  can  be  evaluated  from 
n 

the  results  of  Appendix  (p)  end  will  be  used  as  influence  coefficients  in 
Duhsmcl  Integrals  for  the  analytical  determination  of  the  generalized 
coordinates  ^(t)  and  q^t)  • 
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c)  Derivation  of  Integral  Equations  for  thu  Evaluation  of  the  Generalized 
Coordinates  q^(t)  and*5n(t). 

Consider  an  elastically  lined  cavity  in  the  medium  under  the  action 
of  the  step  shock  wave.  The  total  radial  and  tangential  displacements  of 
the  shell  boundary  con  be  derived  from  the  supcrpooitl.cn  of  the  corresponding 
displacements  of  the  unlined  boundary  plus  the  displacement  produced  by 
unknown  radial  and  tangential  forces  X(t)  and  Y(t)  respectively  which  are 
required  to  force  the  unlined  cavity  into  a  compatible  displacement  with 
the  elastically  lined  cavity.  The  forces  X(t)  and  Y(t)  [Fig,  (11-4)]  are 
expanded  into  a  Fourier  series  in  0: 

00 

X(t)  »^Xn(t)  cos  (11-12) 

n^O 

00 

f  (t)  t.  £  Yn(t)  sin  n0  •  (11-13) 

D*1 

The  compatibility  equation  on  the  radial  and  tangential  displacements  of 
the  shell  can  be  written  in  terms  of  the  generalized  coordinates  q^t)  and 
qV  (t) ,  and  the  I'orces  .  (t)  and  Y  (t).  The  sum  of  the  displacements  due  to 
1)  the  shock  wave  on  the  unlined  cavity;  2)  the  applied  boundary  tractions 
°rr  “  X  (t)  cob  n 0,  ar0  »  0;  3)  the  applied  boundary  tractions  -  0, 

°r0  0  ®i&  mu°t  be  set  equal  to  the  actual  displacements  of  the 

shell  itself  [Eq.  (U-l)-(2)J, 

The  radial  and  tangential  displacements  produced  by  the  traction 
arr  0  °r&  ~  G/  applied  to  the  cavity  boundary,  can  be  evaluated 

in  terms  of  the  known  displacement  coefficients  y^(t)  and  v^(t)  from  the 
*) 

Uihamel  integrals  7 : 

'  The  derivation  of  formulas  of  the  type  given  by  Eq.  (11-14)-(11-15) 
is  shewn  in  detail  in  Reference  [22  f. 
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r  «>) 

W'J  “dr-vA(t't)dt 

O 

r  dX  (t) 

V“J 

o 


(11-lU) 


(11*15) 


In  a  similar  manner,  the  radial  and  tangential  displacements  due  to  the 
boundary  tractions  ofr  -  0*  V,  ■  V  Q(t)  sin  nO,  are  evaluated  in  terms 
of  the  known  displacement  coefficients  v^(t)  and  v^(t)  from  the  IXihamel 
integrals  j 


f  dY (t) 

V-J  ___vB(W)dT 
o 

r  dY„(t ) 

VBJ  irV^4’ 

o 


(11-16) 


(11-17) 


The  compatibility  equations  for  the  shell  displacements  thus  become; 
f  dx  (t)  y?  dY  (t) 

F(t)  +  J  “dr —  wA^t"T^dT  +  J  at . “  <in(t)  +  %(*)  (11-18) 


/dX  (?)  r  U1  \1  /  U\l>/ 

-3T-  vA(t-T)dx  ♦  J  -a-.  vB(t-T)dT  -  +  dn^(t)  (11*19) 


dY  (t) 


^(t) 


Eq.  ( 11-18) •(11*19)  are  0  pair  of  coupled  integral  equations  on  the 

unknown  elastic  forces  X  (t)  and  Y  (t). 

n'  '  n 


The  generalized  forces  (^(t)  and  Qn(t)  which  appear  in  E$*  (ll-6)-(ll-7) 

are  now  determined  as  functions  of  the  unknown  forceB  X  (t)  cos  nO  and 

n' 

Yn(t)  sin  n9.  A  set  of  two  simultaneous  linear  equations  in  X^t)  ^4  Yfi(t) 
are  obtained.  These  equations  are  then  solved  simultaneously  for  XQ(t)  and  ‘ 
YQ(t)  as  functions  of  the  generalized  coordinates  and  : 
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y t) .  j  (V  v  «a,  \) 

(11-20) 

-  g  (qn,  qa,  ©a,  \) 

(11-21) 

Subatitutlng  Eq.  (11-20) -(11-21)  into  Eq.  (U-l8)-(U-19),  a  »et  of  coupled 
integral  equations  on  the  generalized  coordinates  q^  and  ^  are  obtained. 

These  equations  can  be  solved  by  numerical  methods,  using  finite 
difference  theory*  Essentially,  the  integrals  in  Eq.  (ll-l8)-( 11-19)  are 
replaced  by  finite  difference  summations  from  which  recurrence  formulas 
for  the  coordinates  a^t+k)  and'q^t+k)  are  obtained  in  terms  of  their 
known  values  at  previous  time  steps* 

Once  the  generalised  coordinates  q^  and  are  evaluated,  the  displace¬ 
ments  of  the  shell  are  computed  from  Bq.  (ll-l)-(ll-2) *  Similarly,  the 
shell  velocities  and  accelerations  are  also  computed* 

The  direct  stress  (hoop  stress)  in  the  shell  is  evaluated  from  the 
relations 

o  -  ——  t  (11-22) 

(1-V2) 

«  -  |  (v  +  vfl).  (11*23)  * 

Substituting  Bq.  (ll-l)-(ll-2)  Into  tb«  eh  or*  equation*,  the  n^  component 
of  the  direct  hoop  stress  beeches 

ffan  - - y—  (»  +l)q„  co»  nC 

Mn  (l^)a 


(11-24) ' 
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The  n  component  of  the  flexural  stress  on  in  the  ehell  is  given  by 

D 

°B  "  +  ^n)  cob  n0  (U-25) 

n 

vhorc  d  is  the  distance  from  the  neutral  axis  of  the  shell  to  its  extreme 
fiber.  The  total  hoop  and  bending  stresses  are  evaluated  by  summation  of 
tiie  component  a  and  cr^  respectively. 


FIG.II-2 


FIG.  11-3 
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12.  CONCLUSION  AND  SUMMARY  —  DIFFRACTION  EFFECTS 

Sections  (9)»(ll)  have  considered  three  haaic  theoretical  problems 
vhich  nay  ho  of  ubq  in  the  analysis  and  design  of  underground  structure* 
and  their  contents.  These  problems  are  respectively;  l)  the  stress, 
velocity  and  displacement  fields  produced  in  on  elastic  medium  by  the 
diffraction  of  plane  P  and  S  vavco  by  an  unlined  cylindrical  cavity; 

2)  the  determination  of  shock  spectra  for  installations  that  are  mounted 
in  cylindrical  cavities  in  elastic  media;  and  3)  the  stresses,  velocities 
and  displacements  for  the  cose  of  the  diffraction  of  plane  P  ana  S  vavee 
by  a  cylindrical  elastically  lined  cavity,  i.e.  a  structural  shell  vhich 
would  house  an  installation. 

The  shock  waves  have  been  assumed  to  have  plane  fronts  vhich  are 
parallel  to  the  axis  of  the  cavity.  The  plane  wave  solutions  can  be 
applied  directly  to  the  superseiBinic  range  (V  >  cp).  Moreover,  they 
can  be  used  to  construct  by  superposition,  the  solution  to  more  general 
diffraction  problems  in  vhich  the  stress  field  produced  by  the  surface 
explosion  does  not  consist  of  plane  waves  but  is  considerably  more  com¬ 
plicated.  The  problem  of  the  envelopment  of  the  cavity  by  a  shock  wave 
with  a  wave  front  perpendicular  to  the  axis  of  the  cavity,  has  not  a3 
yet  been  treated,  and  remains  to  be  analyzed. 

It  should  be  noted  that  the  limitations  on  the  use  of  an  elastic 
theory  vhich  are  noted  in  Part  (i)  also  hold  for  the  present  cases.  The 
application  of  the  theoretical  results  obtained  in  Sections  (9)~(lt)  is 
valid  only  in  those  portions  of  the  medium  in  vhich  the  assumption  of 
linear  elasticity  is  tenable.  It  may  be  inferred  frem  Appendix  (F) 
that  thi*  is  not  an  unreasonable  assumption  outside  of  the  crater  region 
in  a  good  granite  rock. 
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Specific  conclusions  on  the  uoo  and  applicability  of*  tho  various  theo¬ 
retical  results  arc  given  in  each  Section  and  will  not  be  repeated  horc. 

It  nay  be  noted  hovevor  that  for  a  pressure  vavc  with  a  step  pressure 
input,  the  hoop  stresses  produced  on  the  boundary  of  an  unlincd  cavity  are 
essentially  quasi- static  and  reach  their  maximum  values  after  3-?*  cnvelop- 
nent  tines,  i.e.  the  dynamic  amplifications  are  snail  percentage  vise# 

Results  for  vaves  vith  decaying  pres  sure -time  hiotories  arc  easily  derived 
from  the  corresponding  step  pressure  results  by  neons  of  IXihsmel  integral*. 

*1110  velocities  and  displacements  of  the  cavity  boundary  are  used  as 
input  functions  for  determining  acceleration  and  relative  displacement 
frequency  spectra  for  shock  mounted  installations  in  the  cavity.  The 
shock  spectra  vhich  include  the  diffraction  effects  and  are  consequently 
more  appropriate,  differ  considerably  from  the  free  field  spectra,  particularly 
in  the  high  frequency  ranges.  In  addition,  it  appears  that  the  most 
favorable  conditions  for  shock  effects  will  be  encountered  if  the  equipment 
is  shock  mounted  so  that  it  will  react  to  the  average  (rigid  body)  motion 
of  the  cavity. 

Results  for  an  elastically  lined  cavity  in  on  elastic  medium  cun  be 
obtained  from  the  corresponding  solutions  for  the  unlincd  cavity.  The 
theoretical  method  leads  to  a  set  of  simultaneous  linear  integral  equations 
vhich  may  be  solved  numerically.  This  large  effort  is  currently  under 
way,  and  it  is  hoped  that  the  results  will  be  available  in  the  near  future. 

It  is  felt  that  the  theoretical  problems  considered  here  represent 
a  basic  first  theoretical  development  leading  to  a  rational  design  procedure 
for  underground  structures.  A  great  deal  of  further  theoretical  study  in¬ 
volving  diffraction  effects  in  dissipative  and  granular  (compressible)  media 
will  also  be  required  for  locations  of  practical  importance  and  interest 
in  vhich  the  elastic  theory  will  not  be  applicable. 
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Appendix  A 


STEADY  STATE  RESPONSE  OF  AN  ELASTIC  HALF* SPACE  DUE  TO  A  SURFACE  PRESSURE 
MOVING  WITH  SUPERSEISMIC  SPEED. 


This  Appendix  considers  the  plane  strain  problem  of  finding  displace¬ 
ments  and  stresses  produced  by  the  uniform  motion  of  a  distributed  normal 
load  p(t  -  x/V)  on  the  surface  of  an  elastic  half-spar.e  (Figure  A-l). 

The  speed  V  of  the  loading  is  greater  than  the  propagation  velocities 
c  and  Cp  in  the  medium*  It  is  assumed  that  the  load  has  been  acting  for 
a  sufficiently  long  time  such  that  a  steady  state  situation  is  established. 

Figure  (A-l)  depicts  the  geometry.  The  coordinate  system  x,  z  is  fixed 
in  space,  while  x,  z  is  the  coordinate  system  with  origin  at  the  front  of 
the  load  and  moving  with  the  load  with  constant  velocity  V.  (if  the  two 
coordinate  systems  coincide  at  t  ■  0,  then  x-Vt-x)  z  «  z“. )  The  dis¬ 
placements  and  stresses  are  obtained  by  superposition  from  the  solution 
for  the  uniformly  moving  concentrated  line  load,  Reference  [8].  The 
results  listed  hereafter  apply  for  the  case  where  the  two  Lame  constants 
are  equal,  \  *  \i.  The  horizontal  and  vertical  displacements,  u  and  w, 
respectively  are 

“  ■ ;  [  v  “  -  —  *  k2  «**  p  <«  -  — -p^)] 

».i[  -Kjcota  P  (»  -  ^!S)  +  Kg  p  (t  -  i±f2ffi)] 

4 

The  stresses  are: 


(cot2fl 


1)  Kjp't  - 


X  +  ZCOttt' 


)  -  ScotflKgpft 


X  +  zc< 


2*fc) 
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172- 


*73- 


A  pl(St  of  the  quantities  o/p  and  t/p  is  shown  in  Figure  (2)  •  It 
is  seen  tb*%  except  near  V/cp  «  1,  the  ratio  |  a/p|io  nearly  equal  to 
unity;  j  f/p{ varies  appreciably,  it*  maximum  being  about  0.60.  The 
consequent}**  of  the  *en*itivity  of  the  solution  near  V/cp  *  1  aro  dis¬ 
cussed  in  the  body  of  the  report* 
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Appendix  B 

SURFACE  WAVES  IN  AN  ELASTIC  HALF  SPACE 


At  largo  distances  from  a  disturbance,  the  major  effects  near  the 
surface  are  due  to  Rayleigh  Waves.  This  has  been  discovered  theoretically 
by  Rayleigh  and  ia  well  confirmed  by  seiemological  experience.  It  is 
the  purpose  of  the  present  paper  to  give  suitable  eaipreosions  for  the 
determination  of  Rayleigh  effects  due  to  transient  normal  pressures  on 
the  surface.  Specifically,  closed  form  solutions  are  presented  for  the 
effect  of  a  concentrated  load,  suddenly  applied  and  maintained  thereafter. 
Any  general  pressure  distribution  may  then  be  treated  by  integration  in 
space  and  time. 

The  complete  effect  of  a  suddenly  applied  concentrated  load  has  been 
treated  by  Pekeris  [JL,  a)  ,  who  also  obtained  response  curves  for  certain 
displacements  requiring  lengthy  numerical  integrations#  Ao  noted  by  Pekeris 
(2)i  and  long  ago  found  by  Sommerfeld  [3],  the  surface  effects  can  be 
obtained  by  the  transform  approach  as  contributions  of  certain  poles, 
ignoring  branch  integrals  which  occur  in  a  complete  solution.  This  ap¬ 
proach  is  used  hereafter.  Introducing  the  clearly  appropriate  oimplifi- 
cation  of  considering  the  depth  small  versus  the  radial  distance  from  the 
force,  relatively  simple  closed  form  expressions  will  bo  obtained  for 
stresses  and  displacements. 

By  superposition  this  concentrated  force  solution  can  be  used  to  ob¬ 
tain  approximate  solutions  for  various  typos  of  distributed  pressure  load¬ 
ings  on  the  surface  of  the  half-space.  As  an  example,  expressions  are 
given  for  the  displacement  and  stress  components  caused  by  a  surface  pres¬ 
sure  (with  on  intensity  varying  in  time)  uniformly  distributed  over  a 
circular  disk  of  increasing  radius. 

'  References  ore  listed  at  the  end  of  this  Appendix. 


Analysis 

The  geometry  of  the  problem  is  depicted  in  Figure (l) where  a  con¬ 
centrated  force  P(poeitive  downward,  and  varying  with  time  ao  the  step 
function)  acts  on  the  surface  of  an  elastic  half-space.  Reference  QTj 
(Equations  (3),  (9),  (l6)-(19))  gives  the  Hankel -Laplace  transforms  of 
the  potentials  from  which  the  transforms  of  the  vertical  and  radial  dis¬ 
placements,  u  and  u  .  respectively  may  be  derived: 
z*  r 


»  a+i« 
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where 


hS  -  (#)2  k2  -  (|)2  c2  „  B  ,  #m  .  3c2 

p  C  P  P  P 

ka  -  (I2  +  h2)1/2  k0  a  (£2  +  k2)1/2  i  «.\PT 
M(t)  -  (2t2  ♦  k2)2  -  I*k2i2a0 


(la) 


For  simplicity,  similar  to  ^lj  ,  the  two  Lame  constants  have  been  set 
equal,  \  m  In  order  to  insure  decay  of  the  displacements  at  infinity, 

the  branch  of  the  square  root  must  be  chosen  to  yield  a  positive  real  part. 

The  approximate  solution  for  the  Rayleigh  phase  is  obtained  from 

Equations  (l)  by  utilising  only  the  residues  of  the  integrals  at  the 
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(Rayleigh)  poles,  p  m  t  #  With  7  -  (|)(3  +n/T)1/2  (  reference  DO). 
The  contribution  from  the  poles  is: 

L  P(7)  J  l(\?  -  1)1/2  (u)2-l)1/2J 


U  a  •  ■ 

z  n^ir 
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r  rt*ir 
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U) 

(to  -  1) 

where,  again,  that  branch  of  the  square  root  must  be  taken  which  yields 
a  positive  real  part,  and  whore 
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0)  -  [j  ♦  U  (y2  -  l)1/2]  jy 
V-  [,^v-i>vg/r 

T-  ,  Hr)  .  •  7 
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Since  only  the  Rayleigh  phase  is  being  considered,  these  expressions  axe 
only  applied  to  shallow  depths  (  /«  1)  and  to  values  of  the  time  near  the 
arrival  time  of  the  Rayleigh  wave  (ti  7).  Using  these  restrictions,  Equa- 
tionb  (2)  may  be  approximately  simplified  to  the  final  form: 
us  -  ^(y2  -  i)1/2  Im  [(1  -  2y2)  Z'1/2  +  2y2  Z^2] 

ur  -  Kj_  7  Re  ["(1  -  272)  +  2 (y2  -  1)1/2  {y2  -  i)1/2  Z‘1/2j 

where 
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From  these  cxpresbiono  for  the  displacements,  expressions  for  the  stress 
components  may  bo  derived  by  routine  differentiation  where  !L  «  1, 
t  «  7  is  used  again: 

°rs  “  h  to*  *  1)2  R°  [.2,i3/2  *  Z2V3 

°rr  “  h  (-2?Hr2  -  i)l/a  In  |^3/a  -  <*]  (4) 

°ee  “  h  )  Ro  [ZI3/2] 

arr  -  ^  a  -  2’,S)  H°  jj2y2  +  3)  zl3/2  +  a  -  2r2)  z23^2 1 


where 


16 


(4a) 


The  numerical  evaluation  of  Equations  (3);  (4)  io  simple.  Not  only 
ore  the  expressions  in  closed  form,  but  each  of  the  responses,  for  a 
given  medium, is  solely  a  function  of  a  single  parameter,  the  non-dimen¬ 
sional  time  t.  Figure (2) shows  as  example  the  vortical  stress  a 


(note 


) 


Equations  (3)  and  (4)  may  be  utilizod  to  obtain  approximate  solutions 
for  various  types  of  pressure  loadings  on  the  surface  of  the  elastic  half¬ 
space  t  Consider  as  an  exomple  the  loading  depicted  in  Figure  (3),  A  nor¬ 
mal  load  of  intensity  p  (varying  in  time)  is  distributed  uniformly  over 
a  circular  disk  of  increasing  radius  R. .  Assuming  further  that  R  «  r, 


the  solution  may  bo  written 
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where 
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Equations  (3)  and  (t)  break  down  for  z  «  0,  and  the  integration  procedure 
(5);  (6)  can  therefore  not  be  used  for  the  surface  z  «  0  either*  The 
determination  of  responses  from  transient  distributed  loads  for  a  *  0 
requires  a  separate  approach  starting  from  the  transform  expressions  for 
thfeee  loads,  or  equivalent  procedures  of  going  to  the  limit  z — >-0. 


,  0 


1\.»  J.  ,  _rf  1  /-*<  » cures.  it  U'  n<_  i  t  hoc* .  ccn  ■:  f 

. . -v  .,  t-w  2\.  e.'Jjuio,*  v*  U>.  uCiC^c  *  ,tc£Tul 

r  >,  (  ).  '***  pr>*  ciu1  u  r  ijor  cc  ^utatioxvvd  tifeit  for 

■  "  '  r  *o  .  ■:  J  Interest.  I  i.  vo  *.  don<  nu;  criC'uly , 

i  -iiijorjnx  ip'  y^'clc  r.  for~ulur.  ( SI -noon*  t>  Rule)  for  both 

'  .ire  r'kicri*  .  Inc  ''uriuble  upp_r  ar*u  lr»*er  limits 

,  *•  m;1  i  'ivrj.jr.u.  It  r.  in  adued  i  ^  AriU.)..  to  tm  nur^rl  cn.1  cc"  pi<  ,*J  t;» 


”  '  tcpral  cv<  r  p *ac <■  (in  toe  /triable  \)  vnicn  utili.-co  1 ne 

,  coefficients  c  *'as  replaced  by  u  finite  difference  re.’fc  lor 

VSi 

application  to  Sir.pro r/a  rule  for  integration*  7r*e  rapidly  ooci listing 
ciiuracteriotic  of  the  influence  coefficients  o  ncceocitatea  the  u&e 

&H 

:f  «jui  otre^elv  fine  in  these  calculation;  in  many  caeca  ceverul 

hundred  points  ver*  required  for  each  Integration  corresponding  to  an 
uoper  ijnd  lover  limit  of  R(f).  Curves  for  tneoe  coefficients  o'  are 
given  in  Flg<  (*)-(7)  tend  illustrate  the  necessity  for  the  fine  tteth 
required  iq  the  space  integration. 

Tne  preai>ure-time  cur/y,  p(t),  and  the  range -time  cur.’e,  R(t), 
correspond* hg  to  a  hypothetical  £0  M«T.  surface  explosion  are  cnovn  in 
Fi (o) -(>»)*  As  vc*  are  primarily  interested  in  the  Rayleigh  effect ^ 
generated  oy  the  nigh  intensity  surface  pressures  in  the  neighborhood 
of  Ground  Zvro,  the  pre  score  and  range  curves  vere  cut  o ti  after  a 
certain  time  t^  m  the  compute tioaad  program*  The  integration  in  time, 
using  these  functions  a,*  inputs  vae  alco  trade  using  SJrpeon's  rule  fo, 
cne  integration  procedure. 
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The  computations  vcro  made,  by  necessity  using  on  IBM  7^4  coesputer. 

The  ccuin\«te  numerical  analysis  as  veil  as  the  detailed  computer  program 
vill  be  presented  in  detail  in  a  forthcoming  paper  1 4]  • 

Results  vere  obtained  at  various  depths  below  throe  surface  pressure 
contours  for  the  20  M.T.  surface  burst*  First,  an  elastic  medium  in  which 
e  m  10,000  ft/sec*,  v  -  1/4  was  considered.  The  pressure  contour*  chosen 
were  at  R  »  2000  ft.  [10,000  psi  contour],  R  «  3200  ft.  (2,000  psi  contour] 
and  R  «  500  ft.  [700  psi  contour]. 

rigs.  (ic)-(12)  show  the  attenuation  of  the  peak  stresses  orr ,  Qq0 
and  0  with  depth  at  the  2000  ft.,  3200  ft.  and  5000  ft.  distanced  from 
Ground  Zero  for  the  elastic  medium  with  c  ■  10,000  ft/eoc.  and  a  weight 
of  167  lb/ft^.  These  curves  are  useful  in  determining  the  depth  at  each 
location  below  which  Rayleigh  wave  effects  are  small  and  say  be  neglected* 

The  reader  is  referred  to  the  discussion  of  the  determination  of  Rayleigh 
wave  effects  in  Sec.  (3)  of  the  report  and  to  the  results  presented  in  . 

Fig.  (3-10). (3-13)* 

The  stresses  produced  by  a  20  M.T.  surface  explosion  on  a  olower 
elastic  medium  in  which  c  -  6000  ft/oec.  but  whoso  weight  was  the  oamc, 
i.o.  167  lb/ft?,  voro  also  obtained  for  R  »  3200  ft.  [2000  psi  contour] 
at  various  depths.  Fig*(l3)  show  the  attenuation  of  the  peak  stresses 
°VT9  641(1  axi  wltb  dspth  at  this  location. 

Some  typical  stress  versus  time  histories  for  Rayleigh  effects  are 
shown  for  orr,  Ogg  and  c%%  at  the  indicated  range  and  depth  in  Figs.  [(l4)-(l6)], 
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CE-(?)  INFLUENCE  COEFFICIENT 


°ee[c]  INFLUENCE  COEFFICIENT 


PRESSURE  (psi) 
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STRESS  DUE  TO  RAYLEIGH  WAVE  VS  DEPTH 


o 
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APPENDIX  C 

DIFFRACTION  OF  A  SHEAR  WAVE  (S  WAVE)  BY  A 
CYLINDRICAL  CAVITY  IN  AN  ELASTIC  MEDIUM 


I  INTRODUCTION 

An  infinitely  long  cylindrical  cavity  in  an  infinite  olaotic  homogen¬ 
eous  and  isotropic  medium  iB  acted  on  by  a  plane  shear  \,ave  whose  front  is 
parallel  to  the  axis  of  the  cavity.  The  shock  wave  propagates  through  the 
medium  with  a  constant  velocity  cq  and  carries  the  shear  stresses  tU(t) 
face  Fig.(l)]. 

By  means  of  a  suitable  transformation  of  coordinates,  the  problem  can 

be  solved  by  a  procedure  which  utilizes  the  results  obtained  in  the  solution 

of  the  problem  of  the  diffraction  of  a  P  wave  by  a  cylindrical  cavity  in  an 

elastic  medium.  An  integral  transform  technique  is  used  to  determine  the 

stress  field  produced  in  the  medium  by  the  diffraction  of  the  incoming 

shock  wave  by  the  cavity.  Expressions  for  the  hoop  stress  Oqq  ,  the  radio.'' 

stress  cr  and  the  shear  otross  cr  n  are  derivod  and  numerical  results  ore 
rr  re 

presented  for  the  hoop  stress  at  the  cavity  boundary.  Although  the 

problem  is  considered  for  S  waves  with  a  step  distribution  in  time,  the 
results  obtained  for  this  case  may  be  used  as  influence  coefficients  to  deter¬ 
mine,  by  means  of  Duhamel  integrals,  the  stress  field  produced  by  waves  with 
time  varying  pressures. 

II  CENERAL  PROCEDURE 

Consider  an  infinite  elastic  medium  which  does  not  contain  a  cavity. 

ine  stresses  a  ,  o  and  produced  by  the  incoming  shock  wave  at 

rr  to  uo 

points  lying  on  a  circle  of  radius  "a"  [Pig^2)l  are  given  by  the  expressions 


*ai4- 


■  T  sin  2  5 
rr 

(i) 

ar0  “  *  co0  2  ® 

(2) 

Oqq  ■  -  T  sin  2  ? 

(3) 

The  superposition  of  the  tractions 

0  m  -  t  sin  2  5 

rr 

(4) 

°r0  "  *  T  008  2  ^ 

(5) 

vhich  are  equal  and  opposite  to  those  given  by  Eq.  (l)  and  (2)  on  the 

surface  r  ■  a,  produces  a  traction  free  surface  which  can  bo  considered 

to  be  the  boundary  of  a  cavity  of  radius  "a”  in  the  medium*  The  total 

stress  field  produced  by  the  incoming  S  wave  is  obtained  by  superimposing 

the  free  field  stresses  [Eq.  (l)-(3)l  and  the  stresses  produced  by  the 

application  of  the  surface  tractions  a  and  a  (Eq,  (4)-(5)]  to  the 

rr  rtf 

boundary  of  the  cavity* 

A  transformation  of  the  coordinate  S  can  now  be  made.  This  transform¬ 
ation  brings  the  applied  tractions  of  Eq.  (U)-(5)  into  coincidence  with 
those  of  Eq.  (4)-(5)  of  [20  ]*  Consequently,  certain  of  the  formulae  and 
influence  coefficients  developed  in  [20 ]  for  the  case  of  on. incoming  P  wave, 
may  be  used  to  obtain  similar  results  for  the  case  of  the  S  wave.  Letting, 
[Fig.  (3ajl,  3-6-74  (6) 

and  substituting  Eq.  (6)  In  Eq.  (4)«(5),  the  tractions  become 

orr  «  •  t  sin  2(e  »  */4)«  y  cos  2  6  (?) 

a[6  m  •  t  cos  2(6  -  n/h)m  •  y  sin  2  6 


(8) 
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The  boundary  tractiena  orr  and  [Eq.  (7)-(8)]  ore  e^anded  into  a 
Fourier  Series  in  0  [Fig^3b)] : 


where 


and 


°rr  000  n0 

D>1 

(9) 

°rfl  -K8(t)  sin  n9 
n*l 

(10) 

Jt/^  +  a(t) 

%(*)  “  *7®rr^  cos  n0  de 

n/h  -  a(t) 

(11) 

+  a(t) 

bn8(t)  -  i /or<j(0>  oin  n0  dfl 
*/4  -  a(t) 

(12) 

a(t)  -  cob"1(1  -  f£) 

(13) 

),  (8)  into  Eq.  (11),  (12),  the  expansion  coefficients 

become 

During  envelopment  t  < 


After  envelopment  t  > 


2a 


»2g(t)  -  -  b2fl(t)  «  T 


sin  (2+n)  a]  n* 

— &r-\ 8in  -5 

(1*0 

^sln^Iij8ln^ 

(15) 

n  /  2 

(16) 

(17) 
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It  ahould  be  noted  that  there  is  no  contribution  from  a  series  term  for 

n  «  0,  and  that  since  only  the  series  term  for  n  «  2  has  a  non  zero 

coefficient  after  envelopment,  only  the  n  ■  2  mode  will  contribute  to  the 

long  time  (static)  solution  for  the  stress  field  in  the  medium. 

The  plane  strain  problem  for  the  stress  field  around  a  cavity  of 

radius  MaH  which  has  the  surface  tractions  orr  and  or0  applied  to  the 

n  n 

boundary  r  -  a,  is  now  considered.  The  surface  tractionB  arr  and  arQ 

r.  n 

are  given  by 

a  *  a  (t)  cos  n9  (lB) 

rr  ns 


ari  “  bna^  8ln  a$  (*9) 

n 

As  in  fee],  to  make  the  problem  more  tractable,  the  tractions  arr  and 

n 

arc  first  applied  to  tho  cavity  boundary  as  stop  functions  in  time  r 
r  n 

°rr  "  :U^t)  coe  (20) 

n 


=  fcrU(t)  sin  xfi 


(21) 


The  stress  field  components  produced  by  the  pressure  inputs  E<^.  (20)-(2l) 
are  evaluated  In  Section  III  of  Reference  [20]  and  are  used  as  influence 
coefficients  in  IXihamel  integrals  for  tho  determination  of  tho  corresponding 
components  produced  by  the  true  boundary  tractions  of  Eq*  (18)-(19).  For 
convenience,  the  auxiliary  problem  is  solved  for  the  following  two  sets 
of  applied  boundary  tractions, 


rr 


U(t)  cos  n0 


°r0  -  0 


(22) 


(23) 


As  an  illustration  of  this  procedure,  consider  the  hoop  stress  at 

#  ## 

the  cavity  boundary.  Lat  the  quantities  and  Oqq  be  tho  time 

dependent  portions  of  the  hoop  stress  Qqq  at  r  «  a  produced  by  the 
applied  tractions  given  by  Eq.  (22)  and  (23)  respectively.  The  correspond- 
ing  stress  components  due  to  the  boundary  tractions  of  Eq.  (18)-(19)  are 
given  by  the  IXihamel  integrals! 


The  total  hoop  stress  aQd  produced  at  the  boundary  of  the  cavity  by 
the  incoming  3  wave  is  given  by  the  superposition  of  the  free  field  stress, 
Eq.  (3)i  and  the  stresses  from  Eq.  (24)-(25)i 


°30 


(t)  m  -  T  COS  2 9  + 


oo 

'The  <*>♦»*  <o] 
n  n 


cob  nd 


(26) 


Tho  method  of  superposition  can  also  be  used  in  a  similar  manner  to  evaluate 
the  stress  field,  arr  ,  Oqq  ,  o^q  at  points  in  the  medium  beyond  tho 
boundary  of  the  cavity,  i.e.,  r  >  a. 


Ill  AUXILIARY  PROBLEM  -  STRESS  FIELD  PRODUCED  BY  THE  BOUNDARY  TRACTIONS 

°rr  “  U(t)  cob  n0  and  ar0  *  kU(t)  sin  n0  [Fig.  (4)1* 
n _ n 

The  stress  field  produced  by  the  boundary  tractions  arr  «  U(t)  cos  n0 

n 

and  a  q  *  kU(t)  sin  n9  has  been  determined  in  [20] ,  Section  III,  and  the 
n 

analytical  results  are  given  in  detail  in  that  reference.  The  hoop  stress 
component,  a ^  ,  can  be  evaluated  from  an  Inversion  Integral  of  the  type 
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oeg  (r,e,t) 

n _ 

cos  na 


c  t 
it 


cp-i7 

_L  f  I(r»S)e  1  dt 

2x1  J  i  [FB  +  DEJ  1 
»oo— i  7 


(27) 


vhere 


t»  +  k£]  (2tf  -  {(2n2  +  2n)2g  +  t2}  H^2)(i|)]  + 


♦  [kF  -  D)  [an(n*l)%  H^-^—)  - 


<U*  2o„na  (2)  At* 


r2  ‘*n  N  o^a  '  cr  *Hn-l  ^  cjk  ^  ^ 


(28) 


and 


F  -  [2n(n+l)  -  3t2]  hJ2)(0  -  2t  hJ2|(0 
2 

B  ■  [-2n(n+l)  +  t2]  H<2)(^  l)  +  :>  ^  t  J^c  t) 


(29) 


(30) 


D  -  2n(n+l)  h£2)(1)  -  Sa  t  H^(t) 


(31) 


E  -  2n(n+l)  l)  -  2n  ^  S  H^|(^  0  (32) 

ffrrn  ffr$n 

Similar  inversion  integrals  for  and  aiiTnff  ***  Proeented  in  (20 ]. 

The  integral  of  Eq.  (27)  has  been  evaluated  in  Section  (XV)  of  [20]  for 


the  hoop  stress  at  the  cavity  boundary,  r  ■  a.  The  hoop  stress  Oqq 

v  n  n 

given  by  the  relation 

aQ9 

m  [•  1  +  R  +  I  ]  cos  n0 

vhere  R  and  I  are  given  by  Eq.  (75)  and  (71)  of  p#]  respectively. 


is 


(33) 
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IV  NUHSHIC/lL  RESULTS  AND  CONCLUSIONS 

Numerical  results  are  presented  for  the  hoop  stress  at  r  «  a,  vhich 
is  produced  by  on  incoming  plane  shear  (S)  vave  with  a  step  pressure  dis¬ 
tribution  in  time.  As  in  the  case  of  the  P  wave,  it  was  found  that  the 
maximum  vuluos  of  the  stress  ore  not  materially  affected  by  the  early-time 

streoo  contributions  which  come  from  tho  coefficients  0^  and  , 

n  n 

n  larger  than  2,  and  consequently,  tho  series  in  Eq.  (26)  may  be  terminated 
after  the  n  «=  2  term  for  computational  purposes*  Moreover,  only  the 
contribution  from  the  n  =  2  component  combines  with  the  free  field  stress 
in  Eq.  (26)  to  give  the  asymptotic  long-time  valuo  of  Oqq  . 

Fig.  (9-0)  shows  the  hoop  otreBo  a  at  the  cavity  boundary  for  the 
locations  5  •  0°  end  8  «*  U50. 

The  hoop  stress  produced  by  tho  step  shear  wave  may  be  used  as  an 
influence  function  to  determine  the  hoop  stress  produced  by  a  vave  with 
varying  pressure-time  history  by  means  of  the  Duhmoel  Integral  of  Eq.  (9-l)* 
[a**  Pi 5. 
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KO&gKCUglgg 

The  following  nocaeaclature  is  used  in  this  paper. 


e,  r 


polar  coordinates,  nee  Figure  (l). 


u(r,  t) 


displacement  vector  of  a  point  in  the  medium. 

radial  and  tangential  displacements  of  a  point 
in  tho  medium,  (hoto  that  a  pooitivo  displace¬ 
ment  v  is  outward). 


a 


radiUG  of  cylindrical  cavity. 


an(t),  ^n(t) 


Cp'  Ce 


* T  .c  ,  f  c 
dn'  V  hn'  fln 


Fourier  series  coefficients  for  expansion  of 
stresses  o  and  o  - 
rr 

Propagation  velocity  of  dilatational  and  shear 
vaves  respectively  in  medium. 

Fourier  series  coefficients,  See  Equations ( 102) and  (103) 


d(t) 


r(r,  e,  t) 
f(r,  9,  a) 


} 


Time  dependent  coefficients  for  rigid  body  motions. 

function  of  r,  0 ,  t  and  corresponding  transformed 
function  with  respect  to  tine. 


F,  B,  D,  E 


functions  (See  Equation  (59)- (62)) 


Honhe!  function'?  of  r\nd  second  kind  of 

order  n. 


r),  Q(i,  r)  =  Functions  appearing  in  transformed  displacements 

and  velocities. 


in(y )>  \(y) 


modified  Beooel  functions  of  the  iirst  and  second 
kind  respectively,  of  order  n. 


V  ke 


unit  vectors  in  r  and  0  directions  respectively. 


k 


coefficient. 


t 

* 

t 

p(t) 

u(t) 


=»  time 

=  delay  time  [Equation  (8)]. 

“  variable  pres  sure -time  history  of  incoming 

pressure  vuve,  sec  Section  (VII). 

«  unit  step  function. 
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a(t) 

n 

< 

n 

v 

t 


0 


n 


angle  of  envelopment* 

number  of  circumferential  vavea,  integer. 

ratio  of  otrcoo  conponcnt*  in  incoming  vavc, 
ooe  Fifture(l). 

Load  con ot ant. 

shear  modulus  of  medium. 

Poioeoa*o  Ratio 
n % 

CP 

case  denoity  of  medium. 

radial  and  shear  stresses  respectively  at  a 
point  in  the  medium. 

stress  intensity  of  incoming  wave. 

potential  functions  of  reflected  and  radiatod 
vaves  in  the  medium. 

transform  variable. 


Additional  symbols  are  defined  as  they  occur  in  the  text.  Subscripts  and 

dots  used  vith  displacement  and  potential  functions  indicate  differentiation; 

d2u  ..  d2u  . 
e.g*  u  -  —x  ,  u  ■  — x  ,  etc. 

**  dr*  dt* 


-226 


I.  INTRODUCTION 

An  Infinitely  Iona  cylindrical  cavity  in  on  infinite  elastic  homogene¬ 
ous  and  isotropic  medium  tu  acted  upon  by  a  plane  shock  wave  vhoee  front 
is  parallel  to  the  axis  of  the  cavity*  Tho  shock  vavo  propagates  through 
the  medium  with  a  constant  velocity  and  envelops  tho  cavity  [Fig.  (l)] . 

The  dlreot  stress  components  aU(t)  and  coU(t)  which  aro  respectively  parallel 
and  perpendicular  to  the  direction  of  wave  propagation,  arc  carried  by  the 
shock  vavo. 

An  integral  transform  technique  is  used  to  determine  the  displacements 
and  velocities  produced  at  various  points  in  the  medium  by  the  diffraction 
of  the  incoming  shock  wave  by  the  cavity*  Expressions  for  radial  and 
tangential  components  of  the  displacement  and  velocity  ore  derived  and 
numerical  results  ore  presented  for  these  quantities  at  points  on  the 
boundary  of  the  cavity.  Although  tho  problem  is  conoidorcd  for  pros sure 
vsvee  with  -  step  distribution  in  time;  the  results  obtained  for  thio  case 
may  be  usod  as  influence  coefficients  to  determine,  by  meono  of  Duhamel 
integrals,  the  displacements  and  velocities  produced  by  waves  with  time- 
dependent  pressures. 

II.  GENERAL  PROCEDURE 

The  stress  field  produced  by  the  diffraction  of  a  plane  step  shock 
wave  by  a  cylindrical,  cavity  in  an  elastic  medium  was  evaluated  in  Ref.  [20]. 
The  present  paper  utility*  a  similar  approach  in  determining  the  displacements 
and  velocities  of  mass  points  In  the  elastic  medium.  Essentially,  the 
displacement  and  velocity  consonants  at  any  point  in  the  medium  are  obtained  by 
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superimposing  the  free-Cield  velocity  and  displacement  component*,  i*e 
thoee  produced  by  the  pressure  wave  in  the  medium  vith  no  cavity,  and 
the  velocity  and  displacement  components  produced  by  the  application  of 
corrective  tractions  at  the  boundary  of  the  cavity  in  order  to  make  the 
boundary  turf ace  traction  free. 


a)  Frec-Fleld  Velocities  and  Displacements. 


Consider  an  infinite  elastic  medium  vhich  dcea  not  contain  a  cavity* 

Th*  streese?  ?rr  or£  produced  by  the  incoming  shock  vave  at  points  lying 
on  a  circle  of  radius  ’’a'*,  [Figure  (2aJ  are  given  by  the  expressions 


r  m  -o  cos^O  -  (  iin^0 

(i) 

rr  L  J 

w  -  °  [¥]  °in  20 

(2) 

The  radial  and  tangential  components  of  the  particle  velocity  of 
points  in  the  medium  behind  the  step  vave  front  are  given  by  the  relations 

V -  «  -  —  COB  0  (3) 


vf  *  ~  Bin  6 
f  ocp 


(4) 


The  displacement  cf  a  point  in  the  medium  contains  tvo  components 
-f  “  vf(r»e>t)kr  +  vf(r,Ott)k0,  (5) 

where,  hy  integrating  Eq.  (3) -(4),  the  radial  and  tangential  components 
of  the  dieplacement  are  given  by 


v.(r,0,t)  •  •  rr*tcos  flj(t-t  ) 
r 


(6) 
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vf(i',0,t)  -  ^f-Uin  A]  (t-t*) 
P 


(7) 


whoro  the  time 
point  r  -  a,  0 


t  io  moaoured  from  tho  arrival  of  the 
-  0°  [Fig,  (3)1  and  the  delay  time  t* 

a  -  r  cos  0 

*  — — 


shock  wave  at  the 
is  given  by 

(3) 


b)  Velocities  and  Displacements  produced  by  corrective  boundary  tractions 
applied  to  the  surface  of  the  cavity,  r  ■  a* 

The  stresses  produced  in  an  infinite  elastic  medium  with  no  cavity, 
at  points  lying  on  a  circle  of  radius  "a",  arc  given  by  Eq.  (l)  -  (2), 

The  superposition  of  the  tractionB 

m  cr  £coo^0  -  «  oin^J  (9) 

°re  -  *°  r ’]  8ln  20  (10) 

which  ore  equal  and  opposite  to  those  given  by  Eq.  (i)  and  (2)  on  the 
surface  r  ■  a>  produces  a  traction-froc  surface  which  can  then  be  considered 
to  be  the  boundary  of  a  cavity  of  radius  "a"  in  the  medium.  The  total  dis¬ 
placement  and  velocity  field  produced  by  the  incoming  pressure  wave  io 
obtained  by  superimposing  the  free  field  quantities  [Eq.  (3)-(4),  (6)-(7)l 
and  the  velocities  and  displacements  produced  by  the  application  of  the 
surface  tractions  of  Eq,  (9)*(10)  to  the  boundary  of  the  cavity. 
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Proceeding  as  in  Ref.  [1],  the  boundary  tractions  orr  and  ar5 
[Eq,  (9)  •  (10)]  are  expanded  into  a  Fourier  Series  (Fig,  (2b)  ]t 


a  (t) 

orr  -  — +  )  aa(t)  cos  n0 


n»l 


(n) 


'r0  =  £bn^ 


sin  nf? 


n-1 


where 


«(t) 


■„<*>  ■  !/ »rr(»)  cos  nO  dO 
a(t) 

bQ(t)  •  ar&  (e)  ein  n6  d0 


(12) 

(13) 

(14) 


<  2ft 

For  tines  which  are. less  than  one  full  envelopment,  i.a.,  t  -  —  > 

CP 

the  angle  a(t)  is  given  by  the  expression 


a(t)  »  cos  d-^) 


c_t 


(15) 


For  t  >  —  ,  the  angle  a(t)  «  tt*  Substituting  Eq*  (9)  •  (10)  into 


Eq«  (13)  •  (14),  the  expansion  coefficients  become: 

ta 

i 

p_ 


^  2ft 

During  envelopment  t  -  ~ 


!2(t)  .  £  +  (i££)  sin  2a]  (l6) 

aQ(t)  »  £,ln  wx(co»^a  •  saln^a)  +  (j^J)  8in  (2-n)a  *  (^jj*-)  8ln  ***] 

ja  ^  oin(2+n)aj 


V*>  •  •  h  <*">  Nr-  -  -as 


(17) 


(16) 
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Aftor  envelopment  t  >  ~ 


•»w .  „,1S) 

2  2  ' 

(19) 

an(t)  -  bn(t)  m  0  n  /  0  and  2 

(20) 

«fe(t)  -  -  b8(t)  .  o$£) 

(21) 

For  a  wave  with  a  plane  wave  front,  the  parameter  c  becomes 


€  ■  • 


1-V 


(22) 


The  plane  • train  problem  for  the  displacement  and  velocity  fields 

around  a  cavity  of  radius  "a"  which  has  the  surface  tractions  o  ^  and 

rr 

n 

ar0  to  the  boundary  r  «  a,  i»  now  considered.  Ihe  tractions 

n 

orr  and  or0  are  given  by 
n  n 


"rr  -  “nW  00•  n0 
n 


are  "  *in  n0  J 

n 


2d 


rr  O 
0 


a_(t)  ;  arg  -  0 


(23) 

(2U) 


To  make  the  problem  more  tractable,  the  tractions  tj  and  c  are  first 

rr  rt> 

n  n 


applied  to  the  cavity  boundary  as  step  functions  in  time: 

a  U(t)  cos  n0 


rr 


or^  «k0  U(t)  sin  n  9 


Orr  -  0  U(t) 

O 


(25) 

(26) 
(27) 


vhera  k  1.  an  arbitrary  number  which  is  used  to  Identify  that  part  of  the 


solution  that  cornea  from  the  o, 


r6. 


traction 
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Ths  displacement  and  velocity  components  produced  by  the  traction  inputs 
of  Eq.  (25)  -  (27)  arc  evaluated  in  Section  (U 3) and  arc  used  as  influence 
coefficients  in  Duhaael  integrals  for  the  determination  of  the  corresponding 
components  produced  by  the  boundary  tractions  of  Eq.  (23)  -  (24).  For 
convenience,  the  auxiliary  problem  Is  solved  for  the  following  three  sets 
of  applied  boundary  tractions 


0 rr  -  U(t)  cos  n 0 
n 


(28) 


rr 


•  0 


<*r0  "  8in  n0 
r  n 


(2?) 


arr  -  U(t)  ;  arQ  m  0  (30) 

0  0 

To  illustrate  the  procedure,  let  the  quantities  vn* ,  v  **  and  vQ* 

represent  the  r  and  t  dependent  portions  of  the  radial  displacement  at 
any  point  r,Q,  produced  by  tho  appliod  tractions  of  Eq.  (25)  •  (27) 
repectively.  The  corresponding  displacement  components  due  to  the 
boundary  tractions  of  Eq.  (23)  -  (24)  are  given  by  the  IXihsmel  integrals 


wn(r*t)  ■ J  [ftn  (t)  vn#  (r,t-T)]  dr 

n  >  1 

(31) 

^n(r,t)  « J  (t)  vQ**  (r,t>T)J  dt 

n  >  1 

02) 

r  V  n 

w0(r/t)  -  /i[»0  (t)  w0#  (r»t»t)J  dt 

n  ■  0 

(33) 
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The  total  radial  displacement  w(r,0,t)  of  a  point  in  the  medium, 
which  ie  behind  the  incoming  shock  wave  front  is  obtained  by  superimposing 
the  components  given  by  Eq*  (6)  and  Eq.  (31)  -  (33) J 

v(r,0,t)  -  -  p2-[co»  0](t-t*)  +  vQ(r,t)  +  ^£vn(r,t)  +  *n(r,t)J  co*  nO  (3k) 

In  a  similar  manner,  expressions  for  v(r,0,t),  w(r,0,t)  and  v(r,0,t) 
can  be  obtained  by  superposition  of  the  free  field  quantities  and  those 
produced  by  the  application  of  the  corrective  boundary  tractions  to  the 
surface  of  the  cavity,  r  -  a  i 


v(r,0,t) 

■  sf-t,in  -n 

vD(r,t)  ♦  $n 

(r,t)"|  sin  nd 

(35) 

V 

n*l 

w(r,9,t) 

0 

m  -  — -  COB 

pC„ 

o  +  w0(r,t) 

+  ^[wn(r,t) 

+  vn(r,t)j  cos  nO 

(36) 

P 

n*l 

v(r,e,t) 

■  -2-  .in  6 

ft 

+  £[$n<r,t)  ♦  vn(r,t)] 

oin  n 0 

(37) 

D»1 


r  (3k)  •  (37)  are  given  in  a  form  suitable  for  those  points  1a  the 
medium  over  which  the  shock  wave  has  already  passed  and  which  have 
already  received  signals  from  the  reflected  and  the  radiated  waves 
originating  at  the  cavity  boundary.  For  other  situations,  the  free  field 
and  corrective  traction  components  must  be  superimposed  with  appropriate 
delay  times,  relating  to  the  arrival  of  the  various  waves  at  the  point 


under  consideration. 
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III. 


AUXILIARY  PROBLEM  -  Displacements  and  Velocities  Produced  by  the  Boundary 

Tractions  *  tfU(t)  coa  n 0  and  “  koU(t)  sin  n0  [Figure  (4)]. 

n  n  _  _ 


The  equation  of  motion  for  the  linearly  elastic  medium  is  given  by: 

PV*u  +  (\  +  p)vV'H  *  PJi  (38) 

vhere  for  the  plane  strain  problem,  the  displacement  u  contains  two  components 
u  -  v(r,e,t)kr  +  v(r,a,t)ke  (39) 

Defining  two  potential  functions,  $f(r,0,t)  and  ♦(r,0,t)  such  that 

v(r,e,t)  -  0r  +  p  *e  (**0) 

v(r,e,t)  -  J  -  *r  W 

and  substituting  Equations  (40)  -  (41)  into  Equation  (30);  the  functions  0 
and  *  satisfy  the  following  wave  equations: 

Cp  v2<f  •  (j  (1*2) 

c2V%  -  V  0*3) 

where  _  i 

s -Ih*  ’ -  /T  ""*> 

are  the  velocities  of  propagation  of  dilatational  and  shear  waves  respec¬ 
tively  in  the  infinite  elastic  medium* 


A  transform  with  respect  to  time  is  applied* 


?(r,M)  -  2«  ^'f(r,e,t)*’lftt  at 

w-17 

t(r,6,  t)  »  J dft 


m 

(U6) 


Writing  the  transformed  potentials  0  and  ^  as: 

—  M 

0(r,  8,  0)  «  0(r#  n)coon0 


^(r,  Q,  0)  -  $(r,  fl)oinnO 


the  transformed  vave  equations  becomo 

^rr  +  r  ^r  +  ^  n  0  (** 

%  r 

♦„  +  r  *r  +  ^  ~  \  >♦  "  0  (5 

cs  r 

and  tho  solutions  for  divergent  vavoo  aro  givon  by  the  expressions 

(k  (r*  o  n^  -  A  u(2)  tto\  ('rtnwfl  /  r 


J*n(r,  »,  «)  -  (^)coonO 

P 

T„(r,  0,  n)  -  B  Hi2)  (^)oinnO 


The  coefficient o  and  B  are  obtained  from  the  boundary  conditions  at 

n  n 

the  cavity  boundary  r  -  a: 


au(t)cosn0 


a q  |  ■  kcM;(t)oinn0  (5*0 

-Wa 

The  stresses  ^  and  o^  are  first  expressed  in  terras  of  the  poten¬ 
tial  functions  0  and  ^ 

«„  -  Ww--|  +  (55) 

r 

V  "  “I  ^r0  ^  +  ?  *r  '  *rr  +  “I  *00  1  (56> 
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Transfomlng  Equations  (55)-(56)  and  substituting  Equations  (5l)-(52) 

ints 


Into  the  transformed  equations,  the  constants  Afl  and  Bq  nay  be  computed. 


The  transformed  potential  functions  become' 


7  ca2(B  +  kE) 

“  2X1H  jjiTfB  +  t)Ej 


j(2)  (to) 


H'1-'  Gr)  CO.  n9 
n  °P 


T  oa2(-P  +  k?)  n(2)  /(te\  . 

vn  "  2«.n  pira  +  dej  nn  'cJ  Bln  n 


(57) 

(58) 


where 


F  -  [2n(n  +  l)  -  3&2]  H*2)(0  -  2*H ^  (0  (59) 

2 

B  -  t-2n(n  +  l)  +  \  t2]  H*2)  (^2  5)  +  2^5  H (-2  5)  (60) 

(61) 

(62) 


D  -  2n(n  +  1)  hJ2)(|)  -  2nl  hJ2J  (&) 


E  ■  2a(n  +  1)  H<2)  (f  t)  .  2af  t  H<2>  t ) 


and 


*  ■  ? 


The  transformed  displacements  are  obtained  by  substituting  Equations 
(57)- (50)  into  the  transformed  Equations  (40)-(4l). 


oa  Pn(i»  r)c0B  n0 


wn(r,  o,  n)  »  2xlflJ1  (fb  +  ub7 

„  Q  (E,  r)  sin  n0 
vQ(r,  6,  ft)  -  (FB  +  bB) 


(63) 

(64) 


vhere 


Pn(t,  r)  -  (B  +  kE]{-f  H<2)(f§)  +  &H<2J(£j[))  "  (»  -  kFlt2?^-2-)] 


(65) 


(* )  She  following  results  are  given  for  the  specific  case  of  \  -  |i,  i.e. 
for  a  value  of  Poisson's  ratio  v  »  1/4. 
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The  corresponding  inversion  intogrolo  for  end  may  be  derived  by  re¬ 
placing  the  quantity  ?nU>  r)  cob  n0  with  Qn(t,  r)  sin  nO  in  Equations 
(69)-(70)  respectively. 

The  inversion  integrals  for  tho  velocity  components  and  at  the 
bountUry  of  tho  cavity,  r  «  a,  will  be  evaluated  in  Section  (IV)* 


IV.  INVERSION  OF  EQUATIONS  FOR  THE  RADIAL  AMD  TANGENTIAL  VELOCITY  COMPONENTS, 


v  AND  v  AT  THF  CAVITY  BOUNDARY,  r  -  a. 
n -  n - ' - 


The  radial  velocity  vn  at  r  •  a  Is  evaluated  by  Inverting  the  Integral 


c  t 

u4- 


vn(a,e,t)u  x  «■'  lPnU»5)e 


oc  cos  n(?  *  2*1 
P 


r 

J  (FB  +  DS]  d* 

®-l7 


(71) 


where 


Pn(a,0  -  {[3  +  HE)  (-nH^(t)  +  (0)  -  (D  -  kF]  tnH^2)(^)l)  (72) 


All  computations  were  made  using  a  value  of  Poisson's  Ratio  v  *  1/4;  this 
correopondo  physically  to  a  granite  rock  medium.  For  v  »  1/4,  the  relation 

Cp  *=  /I  cB  (73) 

exists  between  the  diiat&tional  and  shear  wave  velocities  and  this  will 
be  introduced  into  Equation  (72)  in  the  following  transform  inversion. 

The  singularities  of  the  integrand  of  Equation  (7i)  arc  a  branch 
point  at  the  origin  t  *  0,  and  simple  poles  defined  by  the  roots  of  the 
equation  F3  +  DE  m  0: 

4/T  t2(n2  -  1)h£2|  ({)H^2|  (JT  0  +  +  1)  -  3t2]H^(0  *)  + 

(2) 

.+  [-6fc3  +  4nl(l  -  n2)  ]H^(t)H^(/3  0  +  [-6/3  63  +  4/F  t  n(l-n?)  )Hn(t)H^(/3t) 

m 

Consider  the  contour  ADCDHF  shown  in  Figure  (  5  )  where  the  arcs  AB  and 
EF  are  of  infinite  radius  and  the  lines  BC  and  DE  are  branch  cuts.  By  the 
residue  theorem  and  Jordan’s  lemma,  the  integral  of  Equation  (71)  may  be 


evaluated  as 
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Pn(i,  a)o1{V/a 

(FB+  DE) 


BCDE 


P_(5»  a)«llV/a 

'  "057^1 - At  +2Rifj 


(75) 


vhero  J  is  the  sun  of  the  residues  of  the  poles  defined  by  Equation  (7^) 
of  the  pertinent  branch  of  the  integrand. 

The  contribution  to  the  velocity  coenponcnt  vq  of  tho  integration  over 
the  path  CD  is  equal  to  zero.  This  value  represents  the  asymptotic  value 
of  vn  for  very  long  times. 

Tho  contribution  of  the  integration  over  the  branch  cuts  BC  and  DE 
leads  to  a%  branch  integral  which  must  be  evaluated  numerically.  Letting 


<«  ■  ^’/a  <  V.  ■  ^ 

and  using  tho  recurrence  ro lations/  ' 

“  2Jn(0  -  H M(l) 

U^hl***)  m  .2JnU)  -  H*X)(0 
the  branch  integral  becomes: 


(76) 


(77) 


(70) 


where 


< 


12a2(n  +  l)2y2  [^(y)  +  ^(y)]  + 

n  n 

+  3y2t*m2(n  +  l)2  +  12n(n  +  l)y2  +  9y **llA£(/3  v)  +  ^/Ty))  + 

+  x&»Vl>ftJ.1(jr)  +  ^(y))  +  36yU(^in.i(/3  y)  +  ^(*3  y))  + 

-  al*n2(n  +  l)y3i>^In.1(y)In(y)  -  ^(yj^y))  - 

-  12/Ty3[2n(n  +  i)  +  y)iJUTy)  *■  y)l 


> 


(79) 


(#)  "A  Treatise  on  the  Theory  of  Bessel  Functions",  by  o.N.  Watson,  Cambridge 
University  Press,  1952>  Pg.  75* 
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&»(n  ♦  l)y2(2n(n  +  l)  +  3>'2H«2I^(y)  +  l^(y)] 

+  6n(n  +  l)y2  (2n(n  ♦  l)  «■  3y2}t"2I2( /3  y)  +  Kjj(/3  y)J 
<  ■*  12nyU  (A^y)  ♦  J^_x(y)l 

+  36nyU  +  K^/fy)] 

•  6ny3t3y2  ♦  2(n  +  l)2H"2Xn.1(y)ln(y)  -  K^yJt^Cy)] 

.  -  6/3  ay3(3y2  +  2(n  +  l)2)!*!2!^^ y)ln(/3  y)  -  \ml(  13 y\(./J y)) 


and  where  Dn  Is  given  by  eq.  (8l),  Pg.  242. 

Equation  (78)  can  be  evaluated  mxnoricaHy  for*  value 0  of  the  tijao 
fit 

paraewter  -2-  , 
a 

The  contributions  from  the  poles  ^  which  arc  rooto  of  Equation 
(74)  may  be  evaluated  \ry  the  Residue  Theorem: 


(80) 


R  »  2niJ 


where 


2^ 


2/J  (1  -  nk)^  H<f>(6j)  B &  {/3  5j) 
t2/3n(nk  -  l)4j  H*2^)  I^2J  (/3  5j) 

+  [2n(n3t  -  l)  +  3£2Hj  (l^)  !^2)(/3  6j) 

;  3«(1  +  Ij) 


ic  t 


.12/3  t]  a£>  (tj)ag>  (JS  tj) 

+2l3{3(3n  -  2)  +  Ml  -  /3  Sj) 

+3£2(4n  -  3I2)  a£l  (5j)H^2)(/3  £j) 

[8i2(n2  -  1)  +  3£2(l*n  -  3*2)]  (/3  tj) 


(82) 
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'Ric  velocity  ccraponent  «n  then  bcccr\co : 


V 

oc 


°  (I  +  U]  cos 


(03) 


In  a  similar  manner,  the  tangential  velocity  at  r  ■  ^  ia  evaluated 


by  inverting  the  integral  t 

i  *y  «„(», 

2*1  ./  r«  5 


vn(a,  0,  t)n 
oc  sin  a  2*1 


-—17 


[PB  +  DE] 


(84) 


where 


^(a,  0  -  [tB  +  kEll-nH^U))  +  (D  -  kF][-nH^(^fJ  +  i)]}(&5) 


Again/  using  the  relation  v  «  l/4,  i.e.,  *»  /3  cc,  the  singularities 

of  the  integrand  of  Equation  (84)  arc  a  branch  point  at  the  origin  1  =  0, 
and  simple  poles  defined  by  the  roots  of  the  equation  FB  +  DE  *  0, 

Equation  (74).  The  integral  io  inverted  over  the  contour  shown  in  Fig¬ 
ure  (5).  Proceeding  as  in  the  inversion  for  vq,  the  tangential  velocity 
component  Y^(a,  Q,  t)  can  be  evaluated  from  the  relation 


v_(a,  <?,  t)n 

-  -  -  (X  +  R)  sin  a0 


ffC 


The  branch  integral  I  io  given  by 

n  r  Hh  -  MU  • 

f-kr*- * 


yy 

a 


dy 


(86) 


(87) 


whew 


-6n(n  +  l)y2[2n(n  +  1)  +  3y2ll»&2(y)  +  K^(y)]  - 

-6n(n  +  l)y2[2n(n  +  l)  +  y)  +  K^/Jy))  -  12ay4(«2I2<>1(y)  +  ^(y)] 

-36nyU(»i2I2_1(/3  y)  +  ^(/J y)) 

+  6ny3[2(n  +  l)2  +  3y2H’&0_1(y)ln(y)  .  K^yjKjy)] 

+  6/3  ny3t2(n  +  l)2  +  3y2lU2i  t(^3  y)1^^ y)  *  Kh.i^’^"y^Kn^y^ 

n**i 


(88) 
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H 


2*iP. 


-3y2[Un2(n  +  1  f  +  12n(n  +  l)y2  +  9y4H«2I^(y)  +  K^y)] 
-12n2(n  +  l)2y2  Is&^YF  y)  +  ^(TJy)] 

-W*lAJ.1(y)  +  ^(y)] 

-36n2y4[  y)  +  ^(/Jy)] 

+  12y3[2n(n  +  l)  +3y2K«2II>1(y)ll|(y)  -  Kfi_.(y)Ku(y)] 

.+24^  n2(n  +  DyVi^/J  y)  Ia(/3y)  -  K^/3  y)^(/Ty)] 


(89) 


and  D  is  given  by  Equation  (8l).  Equation  (87)  may  be  evaluated  numerically 
n  c  t 

for  values  of  the  time  parameter  -jp  * 

The  contributions  from  the  poles  ^  vhich  ore  roots  of  Equation  (7*0 
ray  be  evaluated  by  the  Residue  theorem: 


E  =  2«5  Ej 


vhere 


2/3  (-n  +  k)^2  (5j)H^(/3  ?j) 

V3  [2n(n  -  k)  +  3k£2)^  )H^^3  5,) 

'»“<»  -  k>!j  «£}  <v  K?>('®  V 

-3n(l  .  i>|*  4\»<*>IJ3 

k.  > 

ie  t 

_ E_ 

>  ie  * 

«S«j) 

1 

+  2(343(3n  -  2)  +  4n(l  -  n2)^] 

+  («»  *  362)  H<£>  (Sj)Hi2)(/3  5j) 

.+/3  f8^n2  *  1)  +  -  3I2)]  Ilf >(0  £j) 
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V.  DISPLACEMENTS  AND  VELOCITIES  PRODUCED  BY  THE  BOUNDARY  TRACTIONS  Orj.  «  U(t) 
_ 0 

AND  n  -  0. 

Tv 

0 


The  displacement  u  for  the  case  n  a  0  has  a  radial  component  only 

u  -  v  (r,t)  k  (91a) 


The  coefficient  wQ(r,t)  is  expressed  in  terms  of  a  single  potential  function 
0(r,t): 


wQ(r,t)  =  0r(r,t) 


(91b) 


Proceeding  as  in  Section (IV),  the  radial  displacement  wq  and  the  radial 

velocity  vQ  can  be  obtained  from  an  inversion  of  the  following  integrals: 

c  t 


and 


t 

®-i  7 

f 

h{2)(s|)  e  df, 

2xip  J 
-00. 

■iy 

S(-2h[2)(0  +  35Hj2)(03 

c  t 

0 

o-iy 

f 

n[2)(i|)  0  4  a  dt 

J 

f 

(92) 


(93) 


-iy 


The  quantities  wq  and  wo  at  the  cavity  boundary  are  evaluated  by  integrating 
Eq#  (92)  -  (93)  with  r  ■  a  over  the  contour  ABCDEF,  Figure (5)*  Proceeding 
as  in  Section  (IV),  the  velocity  wq  is  computed  from  the  relation 


w  (a,t)(i 

— ; - ®  Oc  +  R 

oc 


(94) 


where  the  branch  integred  I  and  tho  oum  of  tho  rooiduoo  R  aro  givon  bolov: 
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c  t 
O 


I  * 


Js. 


J2L 


4  [2K,(y)  +  3yK0(y))2  +  «2(2i1(y)  -  3yi0(y)]a 


(95) 


it 


1Hi2)(t.i)e  * 

yL  «tJ'-rT»F(9  -“Tv 


. .  c  t 

it.i  JL 


(96) 


The  summation  over  j  in  Eq.  (96)  is  carried  out  for  all  values  of  $j 
which  lie  *  the  first  and  fourth  quadrants  and  are  roots  of  the  equation 


3iH^2)(t)  -  2Ilj2)(0  -  0. 


(97) 
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VI.  DiJTERMIKATIOH  OCT  THE  RIOID  BODY  MOTION  OF  THE  CAVITE  BOUNDARY  [Fig.  6]  . 

A  portion  of  the  total  ruction  of  the  boundary  of  the  cavity  consists 
of  a  rigid  body  translation  in  which  the  cavity  maintains  its  shape  and 
translates  in  the  direction  of  propagation  of  the  incoming  wave.  This 
rigid  body  motion  can  be  extracted  from  the  total  motion  of  points  on 
the  cavity  boundary*  Essential ly,  one  must  superimpose  a  portion  of  both 
the  n  -  1  component  of  a  Fourier  expansion  of  the  free  field  motion  and 
the  n  »  1  component  of  the  motion  produced  by  the  corrective  boundary 
tractions,  orr  »•  a^(t)  C0B  <*rQ  m  b^(t)  oin  0  [Eq.(23)l* 

Consider,  for  example,  the  determination  of  the  velocity  components, 
w  anti  v,  for  the  rigid  body  translation  of  the  cavity*  The  free  field 
velocity  components,  wf  and  vf,  Eq.  (3)-(^)  respectively,  which  are 
produced  at  points  on  the  cavity  boundary,  r  3  a,  by  the  incoming  plane 
step  shock  wave  are  expanded  into  a  Fourier  series  in  Q: 


where 


d£(t)  r-' 

vf(a,e,t)  =  — —  +2^(t)  cos 
n«l 

00 

vf(a ,0,t)  =  ^Th*(t)  sin  n0 
n=l 

a(t) 

d£(t)  -  |^f(a,0)  cob  nO  d0 
a(t) 

h^(t)  *  oin  n0 


n0 


(9S) 


(99) 


(100) 


(101) 


and  wf  and  are  given  by  Eq.  (3/-(*0* 


The  velocity  components  v«  and  v-  ,  corresponding  to  n  ■  1,  can 

*1  rl 

he  written  in  the  following  manner: 


.  fdf(t)  ♦  hf(t)] 

\  *  L - 2 - J 

.  [•<£(*)  +  hf(t)-| 

\  '  L  s-2  J 


rdf(t)  -  hf(t)n 

COS  0  +  - £ -  COS  $ 


raf{t)  - 

ein  6  -  — - s - sin  9 


(102) 


(103) 


The  c-refficlent, 


A  +  hi' 


represents  a  deformational  component  of  the 


I  ^1  *  nll 

velocities  while  the  coefficient  - « -  represents  the  rigid  body 


velocity  components,  ror  the  step  shock  wave  under  consideration 


i[(t) 


•  —  i  [a  *  ein  a  coo  a) 
pc 


t  <  — 
-  c 


a 

pc 


t  > 


2a 


(io*0 


h£(t)  «J 


2—  fa  -  oin  a  coo  a] 

pc 


t<*> 

-  c 


t  >  ~ 
c 


f  ,  (*) 

In  a  similar  manner,  the  velocities,  v  and  v  ,  '  produced  by 

C1  C1 

the  corrective  boundary  tractionB  orr  «  a^(t)  cos  0  and  org  *»  b^(t)  sin  0 
can  be  vr.ttten  in  the  form; 


(105) 


(*) 


The  quantities  v  (0,t)  and  v  (0,t)  can  be  evaluated  from  the  formulae 
°i  i 


in  Section  (IV). 
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rd'(t)  +  h’(t)-i  rd-(t)  -  h-(tn 

Vc  (6,t)  ■  d£(t)  COS  8  m  - g - J  008  6  +  - g - - - J  COB  0 

fd?(t)  +  h®(t)*i  rd°(t)  -  h?(t)n 

v  (fi.t)  -  h*(t)  oln  e  a  |^— - g — — — J  oin  8  -  [— - g— = J 


(106) 


Bln  e  (107) 


where  the  terns  containing  the  quantities 
body  components  of  the  velocities. 


-dfo)  -  h°(t?j 


are  the  rigid 


The  total  rigid  body  velocity  of  points  on  the  cavity  boundary  is 
obtained  by  superposition  of  the  rigid  body  components  in  Eq.  (102) -(103) 
and  Eq-  (106)-(107): 

\/(0,t)  m  d(t)  C08  0 
v(0,t)  *•  -  d(t)  sin  0 

where 

rd*(t)  +  d*(t)  -  h£(t)  -  h°(t)n 
d(t)  b  I - g - J 


(108) 

(109) 


(IK) 


In  a  similar  manner,  results  for  the  rigid  body  displacements  and 

accelerations  can  also  be  determined.  The  corresponding  free  field 

f  f 

expansion  coefficients,  d^  and  for  the  displacements  and  accelerations 
are  given  below: 

Displacement  coefficients 

v£(t)  »  d£(t)  COfl  0  >  v£(t)  »  -  h£(t)  Bin  0 
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where 


■ 


ga  1 
2  * 


[• 

rc  t  i 

-f--1 


cob  a  -  sin  a  - 


sin' 


,] 


sin^  a 


h*(t) 


[  ^2  x  ["  a  000  a  +  oin  a  '  ”3 
rc  t 

l+* 


oa 
oc_ 


1 


•sf 


t  > 


2a 


t<sa 

—  c 


t  > 


2a 


Acceleration  coefficients: 


v^(t)  -  d^(t)  cob  6  ;  vx(t)  -  ta^(t)  oin  0 


where 


<^(t) 


hj(t) 


-  O 

“  -  cot  a  cob  a 
pa  a 

0 


—  -  ein  a 
pa  7t 


c  -  C2 


*>? 


•it2 


t  > 


2a 


(111) 


(112) 


(113) 


(114) 
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vn.  numerical  results  and  conclusions  . 

Numorical  results  arc  presented  for  the  displacement  and  velocity 
components  at  various  points  on  the  cavity  boundary  r  ■  a;  these  quanti¬ 
ties  are  produced  by  a  plane  step  shock  wave. 

As  in  tho  cane  of  the  otreoo  [19] ,  it  was  found  that  the  maximum 
value 0  of  tho  displacement  and  velocity  arc  not  materially  affected  by 
the  early-time  contributions  which  come  from  the  terms  for  which  n  is  larger 
than  2,  and  consequently,  the  series  in  Equations  (3*0-(37)  W  generally 
be  terminated  after  the  n  *  £terro  for  computational  purposes. 

Figures  [ (9-M)-(9-l*0]  show  the  displacement  and  velocity  components, 
w  and  w  for  tho  locations  0-0°  and  0  «  90°  on  tho  bound  ry  of  tho  cavity. 
Figurco  [(7  )-(<$)]  of  the  Appondix  show  tho  radial  and  tangential  velocity 
components  w  and  v  at  the  locations  0  -  0°,  22.5°,  45°,  67 *5°/  9^°/  135° 
and  l80°  on  the  cavity  boundary. 

The  rigid  body  motions  of  the  cavity  under  a  step  wave  loading  are 
computed  from  flection  (Vi)  and  the  components  from  the  free  field  and  correc¬ 
tive  traction  (n  -  l)  effects  ore  shown  in  Figures  [(9)-(l0)]  respectively* 
The  total  rigid  body  displacements,  velocities  and  accelerations  are  shown 
in  Figure  (9-lp)  of  Section  {9}. 

The  motions  of  points  in  the  media  which  are  produced  by  the  step  shock 
wave  may  be  used  as  influence  coefficients  to  determine  the  corresponding 
motions  produced  by  a  wave  with  a  varying  pressure-time  history#  For  ox- 
ample,  denoting  the  radial  component  of  the  velocity  produced  at  the  bound¬ 
ary  by  a  step  wave  as  vg(0,  t),  the  rudial  velocity  due  to  a  wave  with 

the  pressure-time  history  P(t)  is  easily  computed  from  the  Duhamel  integral 

t 

W  -  POWB(0,  t)  ♦ J  2|ill  w8{t  •  t)dt 
o 

by  numerical  integration  [  See  Figure  (9-5)] • 


-253- 


\  e.j 

r-  -j 


F  I  G .  I  GEOMETRY  OF  PROBLEM 


FIG.  2  BOUNDARY  TRACTIONS 


-?55- 


■ra 


-257- 


F  I G  . 5  -  INTEGRATION  CONTOUR 


NOTE  l  THE  QUANTITY  d(t)  IS  INHERENTLY  NEGATIVE. 
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FI 6.7  RADIAL  VELOCITY  w(e,t) 


FIG. 8  TANGENTIAL  VELOCITY  V(0,t) 
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appendix  E 

DETERMINATION  OP  SHOCK  SPECTRA  FOR  INSTALLATIONS  IN  CYLINDRICAL  CAVITIES 
IN  ELASTIC  MEDIA# 

I.  Introduction# 

The  motion  of  points  on  the  boundary  of  a  cavity  subjected  to  shock 
vaves,  impart  accelerations  to  installations  which  are  located  within  the 
tunnel  and  are  attached  to  these  points.  In  many  cases,  these  installa¬ 
tions  vill  be  quite  shock  sensitive  and  consequently,  they  may  require 
special  mountings  to  absorb  the  shock  effects  produced  by  the  pressure  waves. 

In  order  to  optimalize  the  design  of  shock  mounted  equipment,  chock 
spectra  are  frequently  utilized.  Two  types  of  spectra  are  particularly 
useful  in  the  present  problem;  l)  Spectra  for  the  peak  relative  displace¬ 
ment  of  the  installation  in  the  cavity,  with  respect  to  the  motion  of  the 
points  on  the  cavity  boundary  to  which  it  is  attached;  2)  Spectra  for 
the  peak  absolute  acceleration  which  is  imparted  to  the  shock  mounted  in¬ 
stallation  in  the  cavity.  The  latter  cpectra  give  data  which  is  required 
for  the  actual  design  of  the  installation  to  be  shock  mounted,  nnmoly  the 
peak  accelerations  to  which  it  vill  be  subjected;  the  former  spectra  give 
clearance  requirements  for  the  mounting  of  the  installation  within  the 
cavity. 

The  theoretical  formulation  of  the  problem  of  the  determination  of 
shock  spectra  of  the  type  described  will  be  presented  in  this  Appendix  for 
waves  which  carry  decoying  pressure-time  histories  as  given  by  Brode  (35l; 


(i) 
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Tho  quantity  PQ  la  the  peak  overpressure  in  the  shock  wave,  and  the 

*•  mm  m. 0 

constants  Dp,  A,  a,  B  and  0  are  defined  in  (35) *  The  displacements  of 
points  on  the  cavity  boundary  for  vaveo  with  the  decaying  pressure-time 
history  are  obtained  from  the  corresponding  displacements  due  to  vaveo 
with  step  pressures  by  IXihamel  integrals  [Section  (9)].  The  spectra 
eurvoG  for  typical  cases  are  nhovn  in  fioetion  (10)  of  tho  main  roport. 

It  is  of  interest  to  note  that  heretofore,  shock  spectra  which  were 
used  in  the  design  of  underground  inotallationo  wore  computed  from  free 
field  input  prooouroo,  and  tho  offocte  of  tho  diffraction  of  tho  ohock 
wave  by  tho  cavity  vero  neglected.  The  opcctra  developed  in  thio  Appendix 
include  the  diffraction  effects;  consequently  they  are  more  appropriate 
and  should  supersede  the  free  field  spectra.  It  will  be  noted  that  for 
shock  mounted  equipment  with  high  frequency  components,  the  free  field 
spectra  are  quite  erroneous  and  the  more  accurate  spectra  presented  in 
Section  (10)  should  be  used. 

II.  General  Procedure. 


Consider  an  installation  which  is  shock  mounted  to  th»  vallr.  of  a 
cylindrical  cavity  in  an  elastic  medium.  Tho  problem  can  be  idealized 
by  considering  a  linear  oscillator  consisting  of  a  concontrated  maos  M 
on  a  linear  spring  of  constant  K.  The  oscillator  is  attached  to  a  support 
which  is  subjected  to  a  motion  equal  in  magnitude  to  the  displacement  of 
the  point  or  points  on  the  cavity  boundary  at  which  the  masG  is  attached 
[Figure  (l)].  Let  the  displacement  of  the  support  and  of  the  mass,  both 
relative  to  a  fixed  datum,  be  U(t)  and  Y(t)  respectively;  the  relative 
motion  of  the  mass  M  with  respect  to  the  support  is  Y(t)  -  U(t*b  The 
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equation  of  motion  of  the  oscillator  is: 

MY  +  K  (Y  -  Uj  »  0  (2) 

Defining  the  frequency  of  the  oscillator, 

®  =  Jl 

Equation  (l)  is  written  as 

Y  +  (o2Y  -  u>2U  (4) 

For  initial  rest  conditions,  i.e.  Y(o)  *  0  and  Y(o)  =  0,  the  solution 

of  Equation  (4)  is  given  by  the  integral. 

t 

Y(t)  =  /uU(t)  sinu)(t  -  t)dT  (5) 

J 

o 

Shock  spectra  showing  the  peak  absolute  acceleration  and  the  peak  relative 
displacement  (relative  to  the  moving  support)  as  a  function  of  the  frequency 
a)  of  the  oscillator  are  required.  For  a  given  value  of  o>,  the  peak 
acceleration  of  the  mass 

|m2(Y  -  U)  I  J 

and  the  peak  relative  displacement 
^  Y  “  U  ^  maximum 

are  evaluated. 

For  this  study,  Brode  pressure  curves  produced  at  different  locations 
by  a  hypothetical  20  M.T.  surface  burst  were  considered,  and  the  result¬ 
ing  displacements,  U(t)  at  various  points,  0  s  ,  on  the  cavity  boundary 
were  determined.  These  were  introduced  as  forcing  functions  into  the  right 
hand  side  of  Equation  (4),  and  for  various  values  of  o>,  the  acceleration 
and  displacement  responses,  given  by  Equation  (6)  and  (7)  were  computed 
and  plotted  as  a  function  of  o>.  These  shock  spectra  are  presented  in 


(6) 

(7) 
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Figure  (10-2  )-(lO-J>  ). 

The  computations  were  made  using  an  IBM  704  computer.  The  analytical 
formulation  of  the  problem  in  detail  is  given  in  Section  III  of  this 
Appendix. 

III.  Analytical  Formulation  of  tho  Problem. 

The  determination  of  acceleration  and  displacement  shock  spectra  may 
bo  formulated  analytically  in  tvo  steps.  First,  the  displacement  U(t) 
producod  at  a  point  0  »  0^  on  tho  cavity  boundary  by  a  shock  vavc  with  a 
Brode  pro o sure  decoy  must  be  evaluated.  Secondly,  tho  values  of  u(t) 
are  substituted  as  input  functions  into  the  right  hand  oide  of  Equation 
(4);  tho  oquatioxi  is  then  integrated  and  tho  oxprcooiono  for  the  maxi¬ 
mum  acceleration,  Equation  (6),  and  relative  displacement,  Equation  (7) 
are  evaluated  for  each  value  of  u>.  The  analytical  formulation  Of  the 
problem  follows: 

a)  Evaluation  of  U(t),  produced  at  a  point,  0  =  0^,  of  the  boundary 

of  the  cavity  by  a  shock  wave  with  a  Brodo  pressure  decay. 

Expressions  for  the  displacement,  U(t)  produced  at  a  point,  0  «  0^t 

on  the  cavity  boundary  by  a  stop  shock  wave  (P  wave)  are  given  in  Appendix 

— 

D#  *  The  corresponding  displacement  curve,  U(t),  produced  by  a  wave 
carrying  a  pressure  with  a  Brode  decay,  may  be  evaluated  by  means  of  the 
IXihomel  Integral 

t 

U(t)  «  I'  P(t)  ---(ly-'t  dT  (8) 

'  0 

* 


S«e  Appendix  (D),  Eq*  (34)-(35}« 
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vhere  ?(?)  is  given  by  Equation  (l),  To  perform  this  integration,  ninth 
order  polynomial  curves  are  fitted  to  the  function  U(t).  Due  to  the 
nature  of  the  U(t)  curve,  a  single  polynomial  fit  is  only  valid  over  a 
certain  region  of  the  curve,  and  hence  several  polynomials  arc  used,  each 
being  valid  in  a  certain  region.  Using  the  notation  below,  the  displacement 


(ri>  / 
u  1  ( 


- 1  * 


c  t 

T  «* 

a 


is  valid  in  the  region  R  which  is  bounded  by  and  r^+^, 

(rj 

Region  R:  <  T  <  ri+1  ,  U  1  (t).  (ll) 

i  -  0,  1,  2,  ...  . 

Jr.) 

Substituting  Equation  (9)  into  Equation  (8),  the  displacement  U  (T) 
vhich  is  valid  in  the  region  R,  r^  <  T  <  r^+^,  is  obtained  by  integra¬ 
tion,  This  quantity  is  written  in  the  form 

(r, )  JA,  r. )  JB,  r, )  (r. ) 

U  1  (?)  =  U  1  (?)  +  U  1  (?)  +  l  1  (?).  (12) 


The  expression  U 


U  1  (?) 
_(A,  O 

T,  A 


(t)  is  given  by  the  expression 


U^'  ^  (T)  ■»  2 


(A,  rj 

U  1  (T) 


U(A*  ^(T)  -  A |  V  [  (a^J  1)  -  \  ^V)]  +  a^rl^A)(T)  (14) 

U»l 

*»'  M  .  |  U)-1  [.  .  „<„  .  T))  ««(,,)  *  ^  J  .-(*  ‘  *J>  <»> 

F^(T)  -  (-l)n'^In  -  <x(y  -  T))V^^(T)  + 


(16) 


(17) 
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« 

V<A>(T)  -  |  £  (-l)Dgj.  lPaB-n 


n«n 


and  the  constants  A,  a,  and  n  are  given  In  term  of  tho  poranatora  of  the 

Broda  curves  (See  Equation  (1)], 


Aa  aa 

A  a  ~  p  ;  a  a  =- 

C  D  o'  C  D 

P  P  P  P 


c  D 

,  JL2 

a 


(10) 


•Die  expression  for  U^D'  ri^(T)  is  obtained  by  substituting  B  for  A 

and  £  for-  a  in  Equation  (l3)-(l8),  where 


B  a 


5a 

CD 

P  P 


tJSL 

CP°P 


(19) 


\ 

The  term  Z'ri'(T)  represents  the  displacement  duo  to  the  atmospheric 
pressure  torn  in  Equation  (l)  and  is  given  by  tho  relation 

ztri)  (T)  =  14.7  £ 


’  a(fl)  t" 
0  L  n  T 
n«0 


(20) 


b)  Integration  of  Equation  (4)  and  tho  evaluation  of  the  maximum  accelera¬ 
tion  [Equation  (6)]  and  the  maximum  relative  displacement  [Equation  (7)]* 

Once  the  displacement  U(T)  is  evaluated,  the  displacement  Y(n?)  is  ob¬ 
tained  by  integrating  Equation  (U)  for  a  particular  value  of  the  frequency 
to*  In  term*  of  tho  non-dimensional  variable  T,  Equation  (4)  bocomco 


^2  +  5^  =  S?U{T) 

vf 


(21) 


where 


m 


-  -  <r  <«> 

The  computed  values  of  Y(T)  are  then  used  to  determine  the  maximum  accelera¬ 
tion,  Equation  (6)  and  tho  maximum  relative  displacement,  Equation  (7). 
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Generally,  Equation  (2l)  can  easily  be  integrated  numerically  by  a 
forward  step  integration  in  time.  Letting  ,fkM  be  the  interval  of  tho  uon- 
dimensional  time  steps,  and  using  the  Itouaerov  procedure *  9  for  second 
order  differential  equations,  a  recurrence  formula  for  Y(T  +  k)  is  obtained: 

2—2  2-2 
Y(T  +  k)  »  -  V  [U(T  -  k)  +  10U(T)  +  U(T  +  k)]  +  [3ft  '  ]Y(t)  -  Y(T-k) 


12  +  k  o 


Equation  (23)  allows  the  determination  of  the  displacement  Y(T*  +  k)  in 
terns  of  the  previously  computed  values  of  Y  at  the  tvo  previous  time 
steps,  T  and  T  -  k.  When  the  integration  is  started  from  the  initial  time, 
T  -  0,  the  starting  formula 


Y(0)  »  0 


12+k  a) 


[U(k)  +  5U(0)] 


must  be  used. 

Stability  considerations^  for  the  numerical  solution  of  Equation  (2l) 
by  tar  Noumerov  recurrence  formula,  Equation  (23),  requires  that  tho  time 
step  k  satisfy  the  condition 

k  <  iiS  (26) 


A  further  requirement  for  accuracy  of  the  solution  results  in  the  use  of  a 

substantially  reduced  value  of  V  from  that  of  Equation  (26).  In  the 

metrical  computations  reported  in  Section  (10),  a  value  of  k  *  was 

0) 

used  in  each  case. 


(*)  See,  Numerical  Methods  in  Engineering,  by  M.  0.  Salvadori  and  M.  L. 
Baron,  Prentice-Hall,  Second  Printing,  1955/  1>a6e  118  ff. 
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For  high  frequency  oscillators  in  which  the  value  of  5  is  largo,  it 
becomes  impractical  to  uqq  a  numerical  procedure  becauoe  of  the  email 
values  of  the  spacing  'V  that  are  required.  Consequently,  for  use  with 
high  frequency  oscillators,  on  analytical  expression  for  Y(T),  valid  iu 
the  range,  0  <  T  <  2,  has  been  developed  fran  the  integral  in  Equation  ($)• 
Consider  the  displacement  U(T)  in  the  region  -  0,  corresponding 
to  the  bounds  0  <  T  <  2.  Pron  the  Equation  (12), 

T5^(T)  *  U(A’  °^(T)  +  U(B»  °^(T)  +  Z^(T)  <27) 

Substituting  Equation  (27)  into  Equation  (5)  and  performing  the  integration,  the 
displacement  Y^(T)  valid  in  the  region  0  <  T  <  2,  is  obtained: 

Y^(T)  *  °^(T)  +  Y^D'  °^(T)  +  Y^1  °^(T)  (28) 

vhere 

^A>  °)(T)  „  V  y£A'  °^(T)  (29) 

v— i  a 

n^l 

y(B.-  °)(t)  .  ^  °^(T)  (30) 

n»l 

9 

Y(Z»  0)(T)  -  1U.7  |  £  0)(T) 

n-0 


(31) 
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Tho  function  'M'  °\l’)  is  given  by 


*£A’  0)(T)  -  a^Vl)0"1^  ni 


(n.c)^M)J£^  [.(.if y(.i)*i§£ 

j-n  ®  L  j  k«J+2  1 

’*  \k  oven] 

,1+1  »  k-1  . 

♦  <-o*  tuf*- 

.  n  |  J+l  m  k  |k  odd/  .  j+2  «  k-1  , 

^  1  <-«*  fu>5  %£ *  (.iF  £  u)T 


J*n 

0) 

L 

[k-J+3  ] 

k-J+3 ' 

k  even  l 

k  odd 

[j  odd  j 

J  even 

i=i£5 


a? 


n  w 


-u>5  £  (-D5  +  (.i)"2"  £  uf^ 


n+1  « 


n-1 


f k^n+2  1 

\k,n  even} 


fknn+2  1 
\k,n  odd) 


(32) 


Tbs  exproooion  for  Y^'  la  obtained  by  aubotituting  B  for  A 

n 

and  0  for  a  in  Equation  (32).  The  coefficient o  D  and  0  in  terms  of  the 
Erode  parameters  are  given  by  Equations  (l9)-(20). 


The  expression  for  Y^Z|  °^(T)  is  given  by  the  following  expressions: 


(33) 


In  performing  the  computations  for  determining  the  shock  spectra  for 
acceleration  and  relative  displacement  shown  in  Section  (10)  of  the  main 
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report,  it  vas  found  advantageous  to  uoe  the  analytical  expression  for 
Y(t)  (Equation  (3  )]  to  determine  Y(t)  in  tho  rungo  0  <  T  <2#  ovon 

for  the  intermediate  range  of  the  frequency  to.  Tho  recurrence  formula* 
Equation  (23)  >  was  then  uced  to  prolong  the  solution  numerically  heyond 
the  range  T  *  2.  In  such  cases,  the  numerical  integration  van  started 
vith  the  formulas 

Y(2)  -  Ygj  Y(2  -  k)  -  Y2.k  (34) 

vhoro  Yg  and  vore  computed  from  Equation  (28). 

The  displacements  U(t)  and  Y(t)  vere  computed  using  an  IBM  704 
computer.  Once  these  displacements  are  known,  the  acceleration  and 
relative  displacement  time  histories  can  be  determined  from  Equations  (6) 
and  (7)  respectively.  In  each  case,  the  maximum  value  of  the  quantity  is 
evaluated.  This  procedure  is  repeated  oyer  the  entire  range  of  the 
oscillator  frequency  co.  The  curves  of  maximum  acceleration  and  maximum 
relative  displacement  versus  the  frequency  cd  are  the  required  shock 
spectra. 

IV.  Determination  of  the  ”Frce  Field1*  Shock  Spectra  for  Comparison  Purposes. 

The  shock  spectra  developed  in  Parts  (I-III)  of  Appendix  (E)  include 
the  effects  due  to  the.  diffraction  of  the  shock  wave  by  the  cavity.  It 

is  of  interest  to  compare  these  spectra  vith  tho.  "free  field"  shock  spectra 
which  heretofore  have  been  used  in  the  design  of  underground  installations. 
These  "free  field”  spectra  are  computed  from  the  effects  of  the  free  field 
input  pressures  only;  all  diffraction  effects  due  to  the  cylindrical  cavity 
are  neglected.  Hence,  essentially  the  "free  field"  spectra  are  developed 
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froca  the  pressures  that  a  gage  placed  at  a  point  in  the  medium  vlth  no 
tunnel;  vould  read.  Analytical  expressions  for  the  determination  of  the 
"free  fiold"  accoloration  and  ,r«lntivo  displacement  oi>ectra  aro  developed 
in  this  Section. 

Let  U(t)  be  the  displacement  of  a  point  in  the  medium  produced  by 
a  plane  step  nhock  vave  with  a  constant  pressure  P  and  a  velocity  of 
propagation  c^.  She  displacement  U(t)  is  measured  in  tho  direction  of 
propagation  of  the  stop  vrave. 

The  particle  velocity  at  the  point,  produced  by  the  step  vave 
becomes: 

u(t)  =  JL  (35) 

The  displacement  curve  U(v)  pivduccu  by  a  wave  carrying  d  pressure 
with  a  Erode  decay,  [Equation  (l))  may  be  evaluated  by  means  of  the 
Duhaael  integral  of  Equation  (6).  Substituting  Equations  (35)  and  Equa¬ 
tion  (l)  into  Equation  (8)  and  integrating  the  dioplac\,aent  U(t)  i0 
obtained: 

■  •  Ao  ’  f  •  i  «  •  » /VH i  m r/P 

u(t)  -  ‘J=y—  (}i  -  l/a)(l  -  o  “*)  +  To  + 

pcpa  L  J 

♦  -$T-  (  n  -  1/P)(x  -  e'^)  +  Te-W  1  (36) 

P°pP  l 

where  the  non-dimensional  time  T  is  given  by 


Th*  coefficients  A,  B,  a,  8  and  \i  are  given  by  Equations  (18)-(19). 


i 
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Ihe  displacement  of  the  oscillator,  Y(T)  io 
Equation  (36)  into  Equation  (5)  and  integrating: 


Y(T)  - 

pcpa 


—2  -aT  —  -2 

(u  -  sill  -  -  SLS — —  — sIcSt  +  -y~ 3  cocar]  + 


+  -jp — 5  [S2^  +  -  -g  1  +  Si-  °a-  "-j-jL  g  JoiiSp  -  — yD‘^3"p  C0Q  Sr  1  + 

a2**2  (c?  +  ET  (a2  +  a?)2  j 

2  ^  2 

+  ---75-  ,  [|&  •  jg  ][L  ~  cocaff  -  pp — ^2  —  oiriSr  +  -?p— coca^T]  + 

pC^P  0+0  {J*~  +  Q  P  +  CD 

♦  G?r  +  -|^g-  J  +  oin  ®  *  ,-g-^-^g  coo  SB?  l 

B3  +  m2  B2  (B2+52)2  (B2  ♦  Sf  )2  J 

(37) 


where  a>  “  —  . 
c 

Using  the  values  of  U(T)  and  Y(T),  the  acceleration  and  relative  dis¬ 
placement  time  histories  con  be  determined  from  Equations  (6)  and  (7)  res¬ 
pectively.  Frequency  spectra  for  the  acceleration  and  the  relative  displace¬ 
ment  of  the  oscillator  mass  may  then  be  constructed  following  the  procedure 
in  Section  (ill)  cf  this  Appendix. 


V.  numeric  .*1  Results. 

Acceleration  and  relative  displacement  ohock  spectra  have  been  obtained 
for  cavities  vhich  arc  subjected  to  ?  vaves  with  Erode  pres  cur'-  decays,  from 
a  20  MT  surface  burst.  Spectra  for  the  rigid  body  motion  (average  motion) 
of  the  cavity  are  presented  in  Secticu  (10)  for  waves  with  Pq  *  6500  psi 
[Figure  (10-2))  and  PQ  -  2000  psi  [Figure  (10-3)1*  In  addition,  spectra 
based  on  the  motion  of  the  cavity  boundary  at  the  points  0-0°  and  G  «  180° 
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are  presented  in  Figure  (10-4)  and  (10-5)  respectively,  for  P0  -  pel. 

For  coaparition  purpotea,  tho  "frae  field"  acceleration  epeotrua  Is 
also  given  in  each  case,  ft ic  reader  is  referred  to  Section  (10)  for  a  dis¬ 
cussion  of  the  results  and  a  comparison  of  the  spectra  developod  in  the 
Appendix(vhich  inalu&i  diffraotion  affactsAvith  the  free  field  ipeotr*. 
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Appendix  F 

ELASTIC  PROPERTIES  OF  GRANITES  UNDER  STATIC  LOADING^ 

I*  FR03SDURE, 

Compression  teats  were  performed  on  two  different  domestic  granites 
of  structural  and  architectural  quality  to  determine  tj-picAl  stress-strain 
diagrams  and  claotic  constants  for  these  materials. 

Specimens  having  the  shape  of  parallelepipeds  of  approximately  three 
inches  height  and  one  square  inch  cross-sectional  area  were  cut  using  a 
diamond  rock  eav.  Electrical  resistance  strain  gages,  [SR-U  AX-51 , 
which  measure  both  longitudinal  and  transverse  strains  were  affixed  to 
the  specimens.  The  specimens  were  loaded  statically  cn  a  Bold"!*  Universal- 
testing  machine. 

Table  (l)  presents  the  description,  composition  find  elastic  proper- 

) 

ties  of  the  materials  used  as  listed  by  the  U.  8.  Bureau  of  Standards.'  ' 
Ta^le  (2)  shovs  the  dimensions,  conditions  of  loading  and  teat  results  for 
the  specimens  tested.  Figure  (l)-(3)  show  the  stress- strain  diagrams 
which  were  obtained  from  the  tests. 

II.  Discussion  of  Test  Results ♦ 

It  appears  from  Tables  1  and  2  that  the  tested  specimens  failed  at 
lower  loads  than  expected.  However,  the  Bureau  of  Standards  does  not  list 


(*)  Ihe  tests  were  performed  by  Dr.  Robert  Heller;  Assistant  Professor  of 
Civil  Engineering,  Columbia  University. 

(**)  See,  D.  W.  Kassler,  E«  Irwiey,  W.  H.  Sligh,  "National  Bureau  of 

Standards  ”,  Paper  No.  R*  P-  X320,  J.  of  Research,  V.  25,  19^0,  Pugo  l6l-S0(> 
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more  than  one  of  tvo  test  results  for  each  materiel  and  consequently  a 
statistical  distribution  of  properties  for  a  particular  typo  of  granite 
cannot  be  deduced  from  the  tabulation*  Tho  test  results  for  each  material 
are  remarkably  veil  duplicated  since  the  specimens  wore  all  cut  from  one 
large  block  and  consequently  have  nearly  uniform  properties.  Hie  ordinary 
grey  granite  has  a  higher  compressive  strength,  but  eomcvhat  lover  averagQ 
modulus  of  elasticity  than  tho  pink  one.  Grain  structure,  in  evidence  in 
tho  grey  gx-anite,  had  no  apparent  effect  though  on  one  specimen  tho  load 
vac  applied  perpendicular  and  on  tvo  others  parallel  to  the  grain,  The 
pink  granite  which  docs  not  exhibit  any  particular  grain  structure  had  a 
high  elastic  modulus  because  of  a  greater  quartz  content. 

All  specimens  failed  in  transverse  tension  forming  columnar  fibers. 
The  shape  of  the  stress- strain  diagrams  is  similar  to  those  obtained  by 
other  investigator/).  9 

Poisson' &  ratio  and  a  shear  modulus  may  bo  computed,  and  Poiccon'o 
ratio  io  found  to  be  variable  for  all  specimens.  A  typical  curve  of 

Poisson' »  ratio  ie  shown  in  Figure  (l).  Same  data  indicating  the  loading 

(*) 

rate  sensitivity  for  these  materials  is  presented  by  Wuerker.'  9  The 
loading  rate  in  the  present  tacts  vac  held  to  approx  lately  5000  lb/min# 
For  higher  rates  of  loading,  the  modulus  of  elasticity  increases  con¬ 
siderably,  but  no  appreciable  increase  in  compressive  strength  is  shown. 
Hydrostatic  prossux’e  combined  with  axial  loading  has  al p»  been  used  by 
some  investigators  where  work  could  be  surveyed  in  the  future  to  determine 
the  necessity  end  f edibility  of  future  tests* 

(*)  See,  R*  0*  Wuerker,  Petroleum  Branch,  AlME,  Doc.  1956  (D552V95.)* 


W.  Kassler,  H.  Insley,  W-  E.  Sllgh,  "national  Bureau  of  Standards, "  Paper  Ko.  R.  P.  1320. 
of  Research,  V.  25,  19^0,  Pgs-  161-20S. 
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